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ABSTRACT 


This thesis presents one method for modelling and analyzing a helicopter design using 
Flightlab. Flightlab is a computer program that provides for engineering design, analysis and 
simulation of aircraft using non-linear dynamic modeling techniques. The procedure to model a 
single main rotor helicopter is outlined using the sample helicopter design in the book "Helicopter 
Performance, Stability, and Control" by Ray Prouty. The analysis procedure contains computer 
program scripts for determining the time response of the helicopter to standard control inputs such 
as a longitudinal impulse, a lateral step, and a pedal doublet. A linear model of the helicopter can 
be extracted from the non-linear model, and a comparison of the time response to the control inputs 
based on these two models is presented. The procedure for conducting frequency sweep testing for 
the linear model is also discussed. This guide to using Flightlab for aircraft modelling and analysis 
is designed to make it easier to use Flightlab for creating additional aircraft models for use in control 
system analysis and additional engineering design. 
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I. INTRODUCTION 


A. BACKGROUND 

The most recent trend in the procurement of helicopters by 
the United States government is to have the helicopter design 
developed and flown using some type of flight simulation 
program prior to development of the first flyable prototype. 
The selection of the Boeing/Sikorsky team to build the new 
Light Helicopter (LH) for the U.S. Army, the Comanche, was the 
result of a competition based solely on the flight simulation 
of an engineering design. No prototypes have yet been built. 

The development of a helicopter flight simulation of this 
type is usually done with a very complex set of computer 
programs, a very powerful and expensive computer system, and 
a full motion based simulator. This type of model also takes 
a great deal of manpower and expertise, is very expensive to 
operate, and it takes a very long period of time to develop. 
This type of flight simulation is not feasible for use within 
the learning environment and resources of a university. 

Most university aircraft design courses are conducted 
without any type of flight simulation model to verify the 
acceptability of the aircraft design being developed by the 
students. Most often the students are required to make 
complex calculations by hand or with the help of many 
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different computer programs to resolve the aircraft's 
stability derivatives. Programs such as MATLAB™ are then 
used to evaluate the matrix representation of the aircraft 
design to determine the eigenvalues and mode shapes of the 
dynamic system response. This type of program is limited to 
analyzing the aircraft as only a six degree of freedom model 
and it can only produce numerical data output which can then 
be plotted graphically. 

The advent of low-cost, powerful engineering workstations 
combined with multi-processing computer systems has led to the 
development of a new computer program called Flightlab. This 
program allows you to develop a helicopter design that 
includes the rotor system degrees of freedom in addition to 
the six degrees of freedom of the body. This allows for a 
more accurate representation of the aircraft and subsequent 
improvement in the analysis of the flight characteristics of 
the aircraft. This model can then be used to create a 
computer flight simulation of the aircraft design. The flight 
simulation model provides real time operation with pilot in 
the loop capability to analyze the model. 

B. OVERVIEW OP FLIGHTLAB 

The Flightlab program is written in a higher level 
language which uses a combination of the C and FORTRAN 
computer languages. Although many of the basic features of 
Flightlab are similar to MATLAB™, Flightlab offers the 
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advantage of being able to do a non-linear representation of 
the aircraft system. The program is written for use on 
computers utilizing a Unix™ operating system. The basic 
computer requirements to run Flightlab are a Silicon Graphics 
International (SGI) workstation with 36 megabytes of random 
access memory (RAM) and 100 megabytes of mass storage space 
for the program. Additionally the program requires that you 
use an Tektronix 4014 terminal emulation window, such as 
xterm, to display any of the plots generated by the 
simulation. 

In order to use Flightlab the user must be familiar with 
the basic features of Unix and the use of Unix editors. There 
are two programs that make up the Flightlab computer flight 
simulation program package. These program differ in the 
method by which users interface with the program to develop 
the aircraft models and simulations. The first is called 
Scope and the second is called Gscope. The only difference 
between the two is the manner in which the user interacts with 
them. In Scope the user interfaces with the program directly 
at the programs command prompt. This requires the user to 
type commands and inputs to develop models and execute 
commands directly within the program. This method requires a 
very high degree of knowledge of command syntax and format. 
In Gscope the interface with the program is through the use of 
graphically displayed windows. The Gscope program window 
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based method is more familiar to most computer users and is 
generally less difficult to use. Gscope offers the additional 
advantage in that an aircraft model can be built by linking 
various graphically displayed model components in a building 
block approach. Data required for each component can be 
entered in the associated object fields and then the Gscope 
program can be used to generate the executable computer code 
that defines the model in the correct syntax and format. 

Once the aircraft model is developed, a program must be 
written to run the flight simulation and develop the desired 
output. These programs are called script files and are 
written in the Scope language. Through the use of these 
scripts and the built in functions, the nonlinear model can be 
solved in real time to determine the forces and moments of the 
aircraft and its components. This solution is used to provide 
the aircraft dynamic response in both the time and frequency 
domains, and you can extract a linear model of the aircraft 
for comparison to other linear system models. 

C. MOTIVATION 

Although the Flightlab program allows for a dynamic real 
time simulation of an aircraft model with the modern 
engineering workstations and computers now found at most 
universities, it still requires a great deal of man hours and 
expertise to develop the model and associated programs. The 
developer of Flightlab provides a user's manual that consists 
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of four sub-manuals: Scope Users Tutorial [Ref. 1] ; 
Component Reference Manual [Ref. 2] ; and a Scope Theory Manual 
[Ref. 3]; and a Scope Command Reference Manual [Ref. 4], 
These manuals are all based on using direct interface with 
Scope and do not discuss how to use Gscope. The Theory Manual 
discusses the aerodynamic or control theory used to develop 
each of the components in FLIGHTLAB and is useful for 
determining what sources are available for deriving data 
needed for each of these components. The Component Reference 
Manual describes the overall modelling and solution procedure 
and each component in detail, including what data is required 
by the program for each component. The Tutorial deals with 
much of the basics of modeling various items, but it does not 
show you how to model a complete aircraft. 

This thesis is an attempt to provide a procedural guide 
for using Gscope to model a helicopter and analyze its flight 
characteristics. This guide can be used as a reference manual 
to reduce the amount of time required to learn the FLIGHTLAB 
program and should make it possible for FLIGHTLAB to be used 
during a one quarter design course to set up a model of a 
helicopter or other aircraft and analyze its flight 
characteristics. This guide can also improve the users 
understanding of the program and therefore allow better 
utilization of its capabilities for developing a model of 
more complicated thesis research topics, such as higher 
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harmonic control for helicopters or Unmanned Air Vehicle (UAV) 
control system analysis . 

The following chapters outline a standardized method of 
modeling and analysis of a helicopter using the example 
helicopter design developed by Ray Prouty. The 
characteristics of this helicopter are presented his book, 
"Helicopter Stability and Control" [Ref. 5:pp.669-682]. 
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II MODEL DEVELOPMENT PROCEDURE 


A. GETTING STARTED WITH GSCOPE 

Gscope is an X windows program that is run in the Unix™ 
operating system environment. This program is installed on 
the Naval Postgraduate School's parallel multi-processor 
computer in the Computer Center Visualization Laboratory 
(alioth.cc) and the Department of Aeronautics SGI Indigo 
computers. The program can be run remotely from any Unix™ 
workstation on campus capable of running X. Table 1 lists a 

Table I UNIX ACCOUNT SETUP 


Running on alioth.cc : 

path must include /usr/local/flightlab and 
/usr/local/flightlab/bin 
setenv FL_DIR /usr/local/flightlab 
setenv GVSDIR /usr/local/flightlab/gvs 
setenv PS_PRINTER 'tiVis' 

add xhost+alioth.cc to local machine for remote access 

Running on Department of Aeronautics Machines : 

path must include /usr/local/flightlab and 

/usr/local/flightlab/bin 

setenv XAPPLRESDIR $HOME/app-defaults 

setenv FL_DIR /usr/local/flightlab 

setenv PS_PRINTER 'hp2p_ps' 

mkdir app-defaults 

Copy /usr/local/flightlab/Scope.ad to 
app-defaults/Scope 









few necessary commands that must be used to specify path and 
environment variables for FLIGHTLAB. 

The Gscope program has five main windows used to develop 
models and programs, the main, model, editor, plot and help 
windows. The program is started by typing the command 
"gscope" with or without an ampersand (&) sign behind it in an 
xterm window. The ampersand allows the program to run in the 
background keeping the xterm window active for other use. 
This is useful for editing script files while viewing the 
model window at the same time. 

The program will start by opening a window that provides 
information about the program and then the main window will 
appear as shown in Figure 1. 



Figure 1 Main Window 

The upper area of this window below the file menu is where 
the commands and output of the Gscope program are shown. The 
box at the bottom part of this window beneath the S> is where 
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Scope commands are entered. The area just above this box 
displays the most recently used commands. Selecting the 
windows menu button allows the user to open the other windows 
that are available for use. 

The model window shown in Figure 2, is used to design the 
model. 



Figure 2 Model Window 


This windows file menu bar allows you to open, edit, and 
save files. The model is built by selecting and placing the 
desired components in the area on the right side of this 
window. The components are then connected and the required 
data fields are entered for each component. The different 
components that can be utilized are shown on the left hand 
side of the window. The components shown in Figure 2 are the 
kinematic components, however, other components can be 
selected by clicking the mouse button over the box labeled 
Kinematic. Other groups include aerodynamic, control, non- 
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linear control, and transducer components, and are shown in 
Figures 22-26 in Appendix A. 

B. MODELLING THE HELICOPTER 
1. Model Hierarchy 

A Flightlab model is usually built from the bottom up, 
i.e, you build subsystems, test them, and then combine them to 
form a complete system [Ref. l:p. 1-2]. This method allows 
you to rapidly reconfigure your model by replacing a subsystem 
within the model with a different variation of that subsystem, 
e.g. a main rotor system with four blades instead of three 
blades. All the components in a subsystem collectively form 
one group within the model hierarchy. This group is given the 
name of the sub-system. Many sub-system groups can exist in 
a higher order group within the model. All of the different 
subsystem groups connected together then form the group named 
"model". 

The remainder of this section describes the model 
hierarchy used for Prouty's sample helicopter (PSH) and can be 
used with minor modifications for most single main rotor 
helicopters. The hierarchy of the model will be presented 
from the top level down to the bottom level. This is the way 
models are presented when first opened in Gscope. For ease of 
understanding, all component names will be italicized, all 
group names will be bold, and all program script file names 
will be underlined . 
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2. The World Group 

The top level in the hierarchy of all models is called 
the world level, and this group for the model of psh is shown 
in Figure 3. 



As shown in Figure 3, the world level consists of the 
model group, and the atmosystem and aerosystem components. 
The model group contains all of the other groups (subsystems) 
that make up this helicopter. The atmosystem and aerosys 
components are system components used to collect information 
from all groups within the model and put it all together in 
one place [Ref. 2:p. 6]. The arrows indicate that the system 
components are connected somewhere in the model group, which 
indicates that data is shared between these groups. 

Each component within a group has a set of mandatory 
data fields which must be assigned a value when building the 
model. These data fields are displayed by double clicking the 
left mouse button on a component symbol in what is called the 
object field window. You may insert the value of the data 
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fields in the space provided or you can use pointer variable 
names to direct Flight lab to look in a data file for the value 
of that variable name. The naming convention for all files 
created for the psh model is to use psh with an appropriate 
extension. Data files are given the .prolog extension, so the 
data file for psh is called psh.prolog . The value of all the 
variables in this file is loaded into the scope program prior 
to execution of the model script file. Using one data file 
for all components in the model allows rapid changes to be 
made, such as changing the location of the center of gravity 
(c.g.), without having to generate a new model script. 

All the data fields required for each component used 
for psh, along with the variable name used to refer to the 
data in the psh.prolog file are listed in the psh.exc file 
provided in Appendix A. This appendix also includes a set of 
five figures that show the symbol for each type of component 
and also the component names used in the following sections. 
Additionally this appendix contains all the program scripts 
that include tables of data values used by aerodynamic 
components. Explanatory remarks about entries in a script 
file are entered by using either a comment marker, //, or by 
using the describe feature of Flightlab. The describe feature 
allows data fields in a program script to have a description 
string added after the input, e.g., stc=zeros(1,6); "stc is a 
1x6 matrix of zeros". The description string will be 
displayed for all variables in a group whenever the describe 
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command is used during execution of scope. This feature is 
highly beneficial for keeping track of the variable names and 
the data to which they refer. Where appropriate, comments and 
descriptions of data fields in the script files for psh have 
been used to make it easier to understand what each line of 
program code means. 

Selecting the connections box in the lower right 
corner of a components object field window will show the other 
components in the model which the selected component is 
connected. The aerosys component is connected to the rotor 
group and the atmoays component is connected to model group. 

3. The Model Group 

The model group is the first level in the model 
hierarchy below the world group. This is where all components 
and subsystem groups that define the helicopter are located. 
Double clicking the left mouse button will take you to this 
level and will display the subgroups within the model group as 
shown in Figure 4. 



Figure 4 Model Group 

4. The Hell Group 

The heli group consists of three sub groups and an 
atmosphere component and inertial component which represent 
the helicopter model as shown in Figure 5. 
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Figure 5 Heli Group 


The atmosphere component models the atmospheric 
conditions and shares this data with other components within 
the group. This component uses the ARDC62 model of the 
standard atmosphere and is listed in atmo. tab . The inertial 
component represents the inertial reference coordinate axis 
system. All velocities and accelerations are measured in this 
frame and are transformed to each component' s frame of 
reference using the appropriate transformation matrix. The 
inertial coordinate system is oriented positive x axis forward 
along the nose of the aircraft, positive y axis to the right 
side of the aircraft, and positive z axis down towards earth, 
a. The Body Group 

The body group is used to model the fuselage and 
tail section of the helicopter as shown in Figure 6. 
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Figure 6 Body Group,PSH 
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The first component used is the dof6 component 
which represents the six degrees of freedom of the fuselage. 
Connected to this is the dmass component which is used to 
model the distributed mass and the effective forces about the 
center of gravity of the helicopter. The data needed for this 
component includes the mass of the helicopter minus the mass 
of the main rotor system, the inertia matrix for the 
helicopter and the value of the gravity vector to be used with 
the model. 

The dmass and dof6 components are connected with a 
translate3 component which is used to locate the dmass 
component at the center of gravity. The location is specified 
by a vector which consists of the fuselage station, buttline 
and waterline station locations of the center of gravity. 

Another translated component connects the dof6 
component with a aero3ds component which represents the three- 
dimensional aerodynamic characteristics of the fuselage. Data 
requirements for the fuselage include the lift, drag, and 
pitching moments as a function of angle of attack and 
sideslip. The values for psh were taken from the charts in 
Prouty's book [Ref. 5:pp. 679-682]. This component requires 
data from +90 to -90 degrees angle of attack and since the 
reference did not provide this data, a user-designed function 
(odli) was used to perform one-dimensional linear 
interpolation of the available data to meet this need. This 
function can be used whenever insufficient data is available, 
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however, you must verify the accuracy of this extrapolated 
data using acceptable aerodynamic theory. 

Two tables of the aerodynamic characteristics of 
the fuselage are required for this component. One is the high 
resolution data for low angles of attack, in this case from - 
25 to 25 degrees using five degree increments. The low 
resolution data table provides the characteristics from -90 to 
90 degrees in ten degreee increments. The aerodynamic 
characteristics must also include cross coupling effects due 
to sideslip from -180 to 180 degrees. The data fields and 
tables for the fuselage component for psh is loaded into the 
scope program by executing the c_wfaero.exe file. The method 
used by Prouty for determining the aerodynamic characteristics 
of his theoretical fuselage shape is presented in the USAF 
Stability and Control Datcom manual [Ref. 6]. 

Similar tables are needed for the aero2d3d 
components used to model the horizontal and vertical tail 
sections two-dimensional aerodynamic characteristics with 
three-dimensional flow. The tables for main rotor blade 
segments aero2d component, include the lift, drag and pitching 
moment characteristics based on angle of attack and mach 
number. The files c_htaill.tab and c_htai!2.tab contain the 
high and low resolution lift and drag tables for the 
horizontal tail, and the vertical tail tables are in the files 
cvtaill.tab and c_vtai!2.tab . The commands to load the 
tables for the horizontal and vertical tail section are 
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horizontal and vertical tail section are contained within the 
appropriate section of the psh.Prolog file. 

Each of the aero2d3d components used to model the 
horizontal and vertical tail sections are connected to the 
dof6 component with translates components and rotate 
components. The translates components provide the location 
of these sections in terms of their fuselage, buttline and 
waterline station. The orientation of the coordinate system 
for aerodynamic components including the horizontal and 
vertical tails and rotor blades is different than the 
inertial coordinate system. .Rotate components are used to 
create the transformation from the inertial frame to the 
component frame of reference. Each rotate component has data 
fields to indicate the axis of rotation and the amount of 
rotation from the previous frame of reference to the next. 
Two rotate components are needed to change the orientation of 
the coordinate axis so that the x axis is to the right along 
the span of the tail section, y is forward along the chord and 
z is up. An extra rotate component is used for the horizontal 
tail to allow adjustment of its angle of incidence. 

The tail rotor is modeled using the Bailey 
component. This component is a simplified model of a tail 
rotor based on the theory presented in NACA Report No. 716 
[Ref. 7]. The orientation of this component's frame of 
reference requires one rotate component to rotate the y-z 
plane so that the z axis points in the direction of the tail 
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rotor thrust. One of the required data field parameters is 
for tail rotor blockage due to the vertical tail. These 
values represents the amount of tail rotor thrust loss due to 
the vertical tail at low speeds. Since there was no data 
available for psh concerning tail blockage, this parameter was 
set to be equal to one, which represents the assumption of no 
losses. 

The final component in the body group is a msensor 
component. This component provides data about the motion of 
the fuselage rigid body to the control system including the 
inertial position, body axis velocities and accelerations, and 
body axis angular velocity and accelerations in the inertial 
frame of reference. The data fields for this component state 
the number of outputs desired and a gain matrix used to select 
which information is provided to the control system. 

Each of the components of the heli group are 
connected to others as indicated by the arrows between them. 
There is a limit to the number of connections that some 
components can have and is listed in the appropriate section 
of the Component Reference Manual [Ref. 2]. Each connection 
between components is defined in terms of the parent node and 
child node for the connection. The correct connection between 
the parent and child frame is important for example in the 
case of a sumj component, which represents a summing junction. 
This component is limited to two connections which are summed 
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together based on an assigned gain value for each, which is 
either plus or minus one. 

Jb. The Rotor Group 

The rotor group represents the helicopter blade 
system, its rotor hub and its swashplate as shown in Figure 7. 



Figure 7 Rotor Group, PSH 


The translate3 component of the rotor group 
provides the same generalized location data in terms of 
stations for the rotor hub as was used for the aerodynamic 
components in the body group. The rotate component changes 
the orientation of the x-z plane such that x is pointed aft 
along blade one at the 0 degree azimuth position and z is 
pointing up. This rotation also includes the tilt of the main 
rotor shaft. Connected to this rotation component is the tpp 
component which is used to compute the tip path plane angles 
for the blades. Also connected to the rotate component is a 
chinge component (chinge) which is used to model a controlled 
hinge which provides the angular velocity input from the drive 
train group to the main rotor blades. This provides the 
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blades with the rotation force to keep them moving at the 
selected main rotor speed (rpm). 

The control hinge is connected to four identical 
blade groups which are in turn connected to the swashplate 
component. The swashplate component is connected to each 
blade at the feathering hinge and is used to provide control 
input for the main rotor system. The data for the control 
inputs is provided to the swashplate by the control system 
without being connected to the swashplate through use of the 
data variable field input. A rigid blade element model was 
used for psh, and a representative blade group is shown in 
Figure 8. 
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Figure 8 Rigid Blade Model,PSH 


The first rotate component represents the 
transformation of the tip path plane coordinate system to have 
the x axis point along the span of the blades with the y axis 
forward along the chord and the z axis up. Obviously the 
rotation required for blade one pointing aft towards the tail 
of the aircraft is zero, and the rotation for the other blades 
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is 90 degrees from the previous blade since there are four 
blades in this rotor system. 

The translate component connected to the rotate 
component is used to represent the distance from the center of 
the rotor hub to the blade hinges, i.e. the hinge offset. 
Connected to this is a torspdm component which represents a 
torsional spring damper that is used to model the lead-lag 
damping in this degree of freedom. The lead-lag damper was 
not originally in the design for psh, but it was added since 
the rotor system model is unstable without a lead-lag damper. 
Prouty does not include one because he assumes that for the 
liner model analysis used in his book there is no effect on 
the stability of his design [Ref. 5:p. 146], but there is an 
effect since we have included the blade degrees of freedom in 
this model. The next hinge component represents the blades 
flapping hinge and degree of freedom about the respective 
axes. The chinge component is where the feathering motion 
input for the blade is input via the swashplate. 

The multiple sets of components following the 
control hinge represent sections of the blade itself. The 
first portion represents the blade spar and is modeled using 
a translate component with a pmass component which represents 
the mass of the blade at that point. Five identical blade 
segments follow, each with translate, pmass, rotate, and 
aero2d components. The end of the blade has a translate 
component which represents the distance from the center of the 
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last blade segment to the tip of the blade. A tip position 
marker, the markpos component, is used by the tpp component 
to compute flapping motion of the blades. The blade segments 
individual data fields are computed using a program file 
called the blade secrment. geom to divide the blade into a 
specified number of segments that sweep out equal areas during 
one revolution. Inputs required for this program include 
blade mass distribution, radius, the number of segments 
desired and the twist distribution of the blade. The model 
for psh uses five blade segments, with a constant mass 
distribution and the blade twist profile as shown in the 
mpl.exc and blade_twist.exc files listed in Appendix A. 

The aero2D components represent the two-dimensional 
aerodynamic characteristics of the blade airfoil and require 
tables for the lift, drag and pitching moment coefficients as 
a function of angle of attack and mach number. No data was 
available for the 0012 airfoil for changing mach number, so 
the same values are used for all mach numbers. Like the 
horizontal and vertical tail section, the data for the main 
rotor blade segments is loaded by the appropriate section of 
the psh.prolocr file. The data tables for the main rotor 
system are listed in the c mrotl.tab and c_mrot2.tab files, 
and were determined using the characteristics for the NACA 
0012 blade as listed in NACA technical note 3361 [Ref. 9], 
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c. The Inflow Group 

The inflow group represents the induced velocity 
through the main rotor system and the interference effects 
between the rotor inflow and the body of the helicopter. The 
inflow group consists of a uniformiv component which models 
the inflow velocity based on momentum theory. This component 
includes the effect of close proximity to the ground. The 
ground effect parameters are determined analytically. The 
inflow time constant for the ground effect, tau, is determined 
by the equation tau=k m /2*omega, where k m =8/(3*pi) and omega is 
the main rotor speed [Ref 10: p. 100] . The ground effect 
parameter equation used for flightlab has the first ground 
effect parameter, gefl, always equal to 0.5 and gef2 = -0.6667 
[Ref 11: p. 147]. The inflow group is shown in Figure 9. 













Figure 9 Inflow Group,PSH 
The four interfer components in the inflow group 
represent the interference between the man rotor downwash and 
the fuselage, horizontal tail, vertical tail, and tail rotor. 
The data field requirements for these components are the 
dynamic pressure ratio, downwash and sidewash effects as a 
function of angle of attack and sideslip, and the incremental 
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velocities due to the rotor as a function of main rotor wake 
skew angle, Chi, and the longitudinal flapping angle of the 
tip path plane, Alf. The values used for these components 
assume no interference or dynamic pressure ratio loss and are 
presented in the files c^wfintf.exc . c_htintf.exe . 
e_vtintf.exe . and c_trintf.exc for the fuselage, horizontal 
and vertical tails, and tail rotor respectively. 

5. The Cont Group 

The cont group includes the pilot's flight control 
components, the mechanical flight controls, and any automatic 
flight control sub-systems for the helicopter. The cont group 
has five sub groups as shown in Figure 10. 



Figure 10 Control Group,PSH 
a. Sensors Group 

The sensors group is used to provide state variable 
input tc the various control groups during execution of the 
model. The contents of the sensors group is shown in Figure 
ll. 
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The three source components at the top of this 
figure output the model's angle of attack, angle of sideslip, 
and the value of unity to the model control subgroups. The 
two euler angles are input to the solution group which is used 
to determine the values of other data fields that depend on 
these angles, e.g., the fuselage aerodynamics. The unity 
source provides a constant value of one as input to the gain 
components which represent the control axis bias in each 
control system subgroup, which is multiplied by a gain factor 
to select the reference position for that control axis. 

The ngain components in the sensor group blocks 
below the sources are connected to the body degree of freedom 
motion sensor. Each ngain acts as a demultiplexer to select 
one of the degrees of freedom from the motion sensor as its 
output. The output of each ngain component is connected to a 
gain component which converts the units of the selected output 
from radians to degrees or radians/second to degrees/second as 
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appropriate. The state variables of the body degree of 
freedom which are "sensed" by these components are phi, theta, 
psi, p, q, and r. 

The solo component connected to the gain component 
in the yaw axis rate sensor is a second order low pass filter. 
This filter is included to give the control group one degree 
of freedom which is required by the control solution method. 

The outputs from these components are set to the 
desired input location in the conficrure.exe file, which is 
listed in Appendix B. 

b. Lat, Dir r Long, Coll Control Groups 

The cont group also includes the subgroups which 
model the control system for each of the four control axes. 
Thes^ subgroups are identical except for the values assigned 
to the component data fields. The lateral control system will 
be used to explain the way the control system groups are 
modeled and is shown in Figure 12. 



Group,PSH 
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The lateral control system model begins with two 
source components which are used to represent the input of the 
pilot and the trim system. The pilot input is set to zero 
except when performing analysis or flying the model, at which 
time it gets its input value from the pilot's workstation or 
the analysis program. The trim source component is set to the 
value determined by the trim control matrix, which is 
calculated during the trim routine for the model. The output 
from these are then summed then multiplied by a gain factor 
which converts the input units from inches of control movement 
to degrees of control system movement. This output is then 
summed again with the control axis bias. The bias for each 
control axis represents the reference control position. The 
reference positions are full left lateral cyclic, full aft 
longitudinal cyclic, full left pedal, and full down 
collective. The value of the bias and the gains for each 
control axis is set in the prolog file along with all other 
object data field variables. The sign of each gain component 
is based on the convention that positive control movement is 
forward longitudinal control, right lateral control, right 
pedal control and up collective control. 

The next part of the lateral control system sums 
the previous output with a check position source. This 
component is used to check the control system model output 
during development of the control systems and may also be used 
by any augmentation control systems developed at a later time. 
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The outputs of the previously discussed components 
are then input to a limiter component, a gain component and 
finally to a sink component. The limiter is used to set the 
upper and lower limit for the control system output in terms 
of the swashplate's limits of travel. The gain component 
converts the units of the output from degrees to radians. The 
output of the source component is used by the swashplate and 
tail rotor components as the input for their respective 
control positions. 

6. The Drivetrain Group 

The drivetrain group is a generic group of components 
used to model the engine and drivetrain of the helicopter, as 
shown in Figure 13. 



Figure 13 Drivetrain Group,PSH 


The design of this group is very simplistic, but it 
functions well as an approximation for an actual system. 
There are no engine and transmission system components modeled 
yet in flightlab so all helicopter models use this drivetrain 
group to command the main rotor speed at the rotor hub. The 
source component models the accelerations of the drivetrain, 
but for this model these accelerations are constant and equal 
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to zero. The inttf components take the output of the previous 
components and integrate to determine the main rotor speed, 
omega, and azimuth position, psi. The two Bixmj components 
concatenate the acceleration, omega and psi values into a 
vector which is output by the final sink component to the main 
rotor control hinge. This component is connected to the main 
rotor hub chinge component as previously discussed. 

C. MODEL SCRIPT 

Upon completion of the model's development the contents of 
the model is saved to a file using the file menu's save 
command. All graphical model files are given the default 
extension mod, hence this file is called the psh.mod . The 
psh.mod file contains the graphical representation of the 
helicopter model as shown in the model window, however, this 
file cannot be executed by the Flightlab program. This file 
can be used to load the model back into the graphical user 
interface for modification or as the basis for creating a new 
model. 

The file menu option "generate script" is used to create 
a file that is executable by Flightlab. This command 
generates the executable file from the model file in the 
correct syntax of the scope language. This is done 
automatically by the Gscope program and requires that the user 
input only the desired name of the file. By convention all 
program scripts generated in this manner are given the 
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filename extension of .exc, so this file is called psh. exc and 
is presented in Appendix A. 
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III. MODEL ANALYSIS PROCEDURE 


A. ANALYSIS PROCEDURE OVERVIEW 

Once the model program script has been generated using 
Gscope you must then create and execute a program file which 
provides scope with instructions for solving the various 
states of the model based on given input and selected 
procedures. As a minimum this program script must load the 
model and associated data, instantiate the system components, 
create solution component structure, initialize the states of 
the model, invoke a solution method for determining the time 
history of the model states, select the desired outputs from 
the simulation, and finally execute the model simulation. The 
structure of the scope program requires a certain number of 
standard things be done in order to analyze and create a 
flight simulation for a helicopter. The following sections 
outline the method used for psh which can be modified as 
necessary for any other single main rotor helicopter model. 
Much of the information presented in the following sections is 
a summary of the detailed explanations found in the 
Flightlab/Scope Component Reference Guide [Ref. 2:p. 12-44]. 

The first section outlines the assembly procedure, i.e., 
how to load the model and component data, configure model 
parameters for analysis and reporting, create a solution 
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structure, initialize the model states, and assemble 
everything together. This procedure must be completed before 
anything else can be done with the model. 

The second section describes the trim procedure which 
explains how to determine the trim conditions in terms of 
airspeed, flight path, and control positions for the 
helicopter model based on selected initial conditions. This 
requires conducting a trim sweep over a range of airspeeds and 
must be performed whenever any significant change is made to 
the operating conditions for the model, e.g., changes in gross 
weight, altitude, main rotor speed. 

The third section outlines the analysis procedure, which 
demonstrates how to determine the response of the model to 
four basic tests: a longitudinal impulse, a lateral step, a 
lateral impulse, and a pedal doublet. Additionally this 
section describes the procedure necessary to obtain a reduced 
order linear state-space system matrix representation of the 
model and compares the output of the above tests for the 
linear and nonlinear simulations. This step is crucial since 
the linear system matrix is needed for control system design 
and analysis studies. The final part of this section also 
outlines how to determine the frequency response 
characteristics of the linear model. 


32 








B. ASSEMBLY PROCEDURE 

The assembly procedure for psh model is listed in the file 
psh.def found in Appendix B. A file with the .def extension 
is by convention a file that defines a sequence of 
instructions and other script files to execute. The psh.def 
file contains all the instructions needed to set up the psh 
model for running and must be accomplished prior to the 
execution of any other simulation script file. 

This file initially sets the path to all directories for 
the files used to assemble the model. This step is important 
since scope does not use the previously defined path for the 
unix system. This procedure also allows the use of files 
previously written for other models so they do not need to be 
copied into the current model directory. The psh sub¬ 
directory under flightlab contains all the script files used 
for the psh and can be used for new aircraft modelling and 
analysis. 

The next step in assembling the model is to load the user 
defined functions (UDF) used during the assembly procedure. 
UDF's are used to define a specific function and are 
constructed from built-in functions which are part of the 
scope language and other UDF's which have been previously 
created [Ref. l:p. 11-22]. UDF's are similar to matlab.m 
files and add to the versatility of the scope program. UDF's 
are usually given the .fun filename extension. All UDF's must 
be executed prior to calling that function in the script 
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files. The assembly procedure uses the cycle and odli 
functions. The cycle UDF defines the method used to cycle 
through the iterative solution process and the odli UDF 
performs one-dimensional linear interpolation as previously 
discussed in the modelling section. 

The third part of the assembly procedure is constructing 
the model in the correct hierarchy for the scope language. 
This is accomplished by executing files that contain these 
instructions. The psh.prolog (data) and psh.exc (model) files 
load the model and data parameters into the world level. The 
psh.epilog file then equivalences the model variable names to 
standard names for the solution and system components 
variables and also initialize the control connections and 
motion sensor gain matrix. The system.exe file creates and 
connects the system component to the model, and finally the 
solution.exe file which sets up the solution components and 
connections. 

The solution file creates the solution configuration for 
psh which includes six solution components. Each solution 
component integrates the states and propagates in time its 
associated model group. The helisolve, rotorsolve, and 
rhsolve components use the numeric method, hsolve, to compute 
the states of the heli, rotor, and the combined rotor and heli 
groups. The drivesolve and contsolve components use the 
analytic method, csolve, to compute the states of the 
drivetrain and cont groups respectively. The final solution 
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component, topaolve, uses the fully coupled numeric and 
analytic method to solve for the states of the entire model 
group. Each solution component method also sets the value of 
required data fields, an explanation of which can be found in 
the component reference manual [Ref. 2: p.31-32]. 

The final part of the assembly procedure is to initialize, 
configure and setup the model for running the simulation. 
This involves using several built-in scope functions to 
initialize and invoke the model. The init command links, 
equivalences and analyzes the data flow for all the components 
in the model. The world::setup command initialize the states 
and methods for the model. The world::reset command sets the 
initial conditions for the states and their derivatives and 
invokes the model. 

The mbc.exc file is executed in order to set up a multi- 
blade coordinate transformation of the rotor states. This 
improves the speed and accuracy of the solution of the rotor 
states during execution of the model. The confiqure.exe file 
is used to configure the model structure for the desired 
reporting and sharing of information between solution 
components by creating the cpg (compute parameter group) and 
results groups at the world level. These groups provide a 
central location for the output of the simulation to be 
stored. 
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The final part of the assembly process is the execution of 
the assemble.exe file. This file sets several of the compute 
flags (cf) to assemble and reset the various solution groups. 
The psh model is now ready for execution. 


C. TRIM PROCEDURE 

The trim procedure script, trimsweep.def . is used to 
conduct a trim sweep for the model over a user specified range 
of airspeeds. The trim sweep for psh was done from 0 to 140 
knots at sea level standard day conditions. Executing the 
trimsweep.def file and specifying the desired range of 
airspeeds and the flight path angles is all that is necessary 
to complete the trim procedure. Upon completion of the 
trimsweep program the user is given the option of displaying 
and printing the results of the trim sweep. All the files 
used for the trim sweep procedure are listed in Appendix C. 

The first file executed by the trimsweep program is the 
psh.def file. This is done to assemble the model for running 
as explained in the previous section. Once the model is 
assembled, the UDF, limitchange.fun , is loaded into the scope 
program. The limitchanae.fun is used by the trim program to 
limit the amount of control input changes to a max of 2.5 
percent of control travel during the trim sweep. This reduces 
the time needed for the trim program to converge to a 
solution. 
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The TrimSweep.exc file is executed next. This file sets 
the number of rotor revolutions used to determine the average 
value of the body accelerations and then executes the main 
trim script file, Trim.exc . 

The Trim.exc file is based on a simple algorithm designed 
to reduce the steady state translational and angular body 
accelerations, (bacc), to zero. The trim process calculates 
t' i initial bacc and then runs a trim routine that iteratively 
calculates the trim control positions changes needed to reduce 
the accelerations and then re-evaluates the body 
accelerations. The iteration cycle continues until the 
largest singular difference from the previous to the current 
body acceleration is less than the convergence limit, 0.0001, 
or the maximum number of iterations is reached. Three trim 
loop iterations are used, with a maximum of 60 iterations. 

After setting the number of rotor revolutions to use for 
determining steady state conditions, a trim matrix is computed 
at each airspeed in the trim sweep. The trim matrix is a 
diagonal matrix which is the partial of one acceleration which 
is coupled to one control as shown below: 
th ~ ud (pitch attitude couples with longitudinal accel) 


ph ~ vd (roll attitude couples with lateral accel) 
xc ~ wd (collective couples with vertical accel) 
xa ~ pd (lateral cyclic couples with roll accel) 
xb ~ qd (long. cyclic couples with pitch accel) 
xp ~ rd (pedal couples with yaw accel) 


The negative inverse of this matrix is used to determine 
the control change per unit acceleration used during the 
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iteration process. The diagonal matrix for each airspeed is 
converted into vector form and concatenated into a single trim 
matrix, trmd. This matrix is saved to Trim_Matrix. rbe at the 
end of the trim sweep for use during the analysis procedure. 

Two other scripts files are executed from within the 
Trim.exc file . Update_CW.exe and Update,exc . These files 
are used to determine the amount of control change to apply 
for each iteration and then to determine the new steady state 
body accelerations and update the control positions. Upon 
completion of the iteration for each airspeed the trim control 
positions are formed into a vector and concatenated to a 
matrix called stc. This matrix is saved to a file called 
TrimControls.rbe at the end of the trimsweep procedure. 

A graph of the bacc is presented after each iteration 
within the trim loop along with a listing of the average bacc 
and current control positions. At the end of the trimsweep 
procedure, the user is given the chance to view a plot of 
control positions versus airspeed as shown in Figure 14. 
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The values for the control positions obtained by this trim 
method for 115 knots are very close to the values obtained 
analytically by Prouty [Ref 5:pp. 527-529]. The plot also 
shows that the sample helicopter runs out of control power at 
approximately 135 knots, which would represent the maximum 
speed psh is capable of achieving with its control system 
design. 

O. ANALYSIS PROCEDURE 

The analysis.def file contains a sample method of 
analyzing the time response and frequency response of a 
helicopter at a given flight speed and condition. The 
analysis .def file sets up a method for analyzing the time 
response and frequency response of psh to control inputs. The 
time response procedure is set up to find an open loop 
nonlinear model solution based on the rhsolve solution group 
and compare that to the response of a reduced order linear 
model solution based on the topsolve solution group. The 
frequency response procedure analyzes the model based on the 
linear model solution only. 

The analysis.def file executes the psh.def to set up the 
model for running and then it adds the test directory to the 
search path for the files used to conduct the various tests. 
The next part of the script executes three additional UDF's 
required for the analysis procedure. The qsreduce.fun is used 
in the 6dof linearize.exe file which reduces the nonlinear 
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model with 37 states to a linear model with six degrees of 
freedom and eight state variables. The linear state variables 
include u, v f w, p. q. r. theta, and phi. The 
6dof_linearize.exe file can be modified to return a linear 
model with more states and degrees of freedom, such as a 10 
degree of freedom model which includes the blade motion 
degrees of freedom. The logspace.fun UDF and the freq.fun UDF 
are used to set up a logarithmically spaced input frequency 
vector and the frequency sweep test respectively. 

The linearization procedure uses the convolution integral 
method to determine the linear characteristic and control 
matrices. The coefficients of the characteristic equation, 
and the eigenvalues and vectors for the characteristic matrix 
are listed at the end of Appendix D. Upon completion of the 
linearization program, a plot of the eigenvalues for the 
linear model is provided as shown in Figure 15. 













As shown in the plot of eigenvalues, there is one pair of 


complex roots that are unstable roots. It will be shown 
later that this pair represents the longitudinal phugoid 
response mode. All of the other roots are stable, however, 
there is another complex pair that is just barely on the 
stable side of the real axis. The stable roots represent the 
longitudinal short period mode, the lateral-directional dutch 
roll mode, and the lateral spiral and roll modes. It is 
difficult to determine which root corresponds with which model 
just from the plot, so we can conduct several standard control 
input tests to determine this. 

The first test conducted by the analysis.def file is a 
longitudinal impulse test. The Long Impulse.exc file is used 
to set up a longitudinal cyclic control input with the user 
defined parameters for duration of the run, size of the input, 
and input delay time. Figure 16 shows a comparison of the 
nonlinear response and the linear response of psh in terms of 
pitch attitude and forward airspeed to a 1 inch forward cyclic 
impulse. 
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As shown in the figure above, the response of the 
nonlinear and linear models correspond well. The change in 
pitch attitude and velocity in response to the impulse input 
is oscillatory and divergent. The rate of divergence is slow, 
however, which corresponds to the location of the unstable 
roots on the eigenvalue plot. The results of this test show 
that psh is longitudinally unstable, which is the same 
conclusion reached by Ray Prouty using his analytic methods 
[Ref. 4: pp. 616-623]. 

The second test of the helicopters response is a lateral 
step input. The test script, Lat_Step.exe asks the user to 
input the duration of the test, the size of the input, the 
duration of the input, the time delay of the input, and the 
rise and fall time of the input. Figure 17 shows the response 
of psh in terms of roll attitude and roll rate to a l inch 
right cyclic step input lasting for 2.5 seconds with the 
input delay time, rise and fall time all set to 0.025 seconds. 



Figure 17 Roll Response 
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This plot shows a large difference between the nonlinear 
and linear response of the model after the first six seconds. 
The nonlinear model appears to have a much higher frequency at 
which it oscillates and may also be divergent. The linear 
response shows a much lower frequency of oscillation and it 
appears to slowly converge. This type of response is expected 
based on the location of the roots of the characteristic 
matrix as shown on the plot of eigenvalues, since we have 
already identified the phugoid mode as the unstable mode. The 
bank angle and roll rate response shown in the above plot for 
the nonlinear model includes the effects of control cross 
coupling. 

The linear response does not correspond well with the 
nonlinear response because the method used to reduce the 
nonlinear model to a linear model is based on the assumption 
that the nonlinear model response remains in the linear range, 
which the data shows in not the case. This type of response 
would not be seen when applying this type of test to a purely 
linear model, because this type of model looks at the response 
of the system with no control cross coupling. The major cause 
of the roll response shown for psh appears to be a coupling of 
pitch to roll. The pure lateral cyclic input to the right 
causes a corresponding pitch up and when the cyclic is moved 
back to the left it causes a corresponding pitch down. The 
coupling effect of roll to pitch in terms of pitch attitude 
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and pitch rate due to the lateral step input is shown in 
Figure 18. 



Lateral Step Input 

The third test of the helicopters response is a lateral 
impulse input. The test script, Latlmpulse.exc . allows the 
user to input the duration of the test, the time delay of the 
input, and the size of the input. Figure 19 shows the 
response of psh in terms of roll attitude and roll rate to a 
l inch right cyclic impulse input with an input delay time of 
0.025 seconds for a time interval of 30 seconds. 
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As shown in this plot, the response of the nonlinear and 
linear models correspond very well with each other. The 
spiral mode response of psh is considered to be stable because 
both the bank angle and roll rate are oscillatory and 
convergent. 

The fourth test of the helicopters response is a pedal 
doublet input. The doublet input is created by using two step 
inputs in opposite directions. The test script, 
PedDoublet.exc . asks the user to input the duration of the 
test, the time delay of the input, the size of the step input, 
the duration time of the step inputs, and the rise, fall and 
delay times between each step. Figure 20 shows the response 
of psh for a 20 second interval in terms of roll attitude, yaw 
attitude, and roll rate to a 1 inch pedal doublet input. Each 
step input lasts 1.5 seconds and the delay, rise, and fall 
times are all equal to 0.025. 
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This plot also shows a large difference between the linear 
and nonlinear response of the model. The nonlinear model 
again appears to oscillate at a much higher frequency than the 
linear model and it also appears to be divergent. The linear 
response shows a much lower frequency of oscillation and it 
appears to be convergent. This type of response for the 
linear model is reasonable based on the location of the linear 
model eigenvalues. The response shown in the above plot for 
the nonlinear model also includes the effects of roll to pitch 
control cross coupling. 

The final test conducted during the analysis procedure was 
to evaluate the frequency response to longitudinal sinusoidal 
inputs over a range of input frequencies. The frequency sweep 
test is conducted by using the freq.fun UDF. This function 
uses three inputs to determine the frequency response of the 
aircraft, the linear system matrix, the number of states in 
the linear model, and a vector of frequencies for the 
sinusoidal inputs. A vector of logarithmicly spaced 
frequencies from 0.1 to 100 radians/second was created using 
the logspace.fun UDF. This function was used to make creating 
a bode plot of the response possible, because one of the 
limitations of the scope program is its lack of a good log 
scale plotting feature. The system matrix was constructed by 
splitting the original linear matrix into its F, G, H, and D 
matrices, and then changing the input and output matrices to 
select only a longitudinal input and pitch attitude as the 


46 









output desired. The system matrix was then reconstructed using 
the original characteristic and control matrices and the 
modified input ar.I output matrices. 

The freq.fun UDF returns the amplitude and phase of the 
response for the given range of input frequencies. In order 
to create a bode plot the amplitude and the frequencies must 
be converted from natural logarithms to logarithms to the base 
10. This is done by dividing the natural log of both the 
amplitude and frequency by the natural log of 10 and 
multiplying the amplitude by 20 to convert to gain in 
decibels. The frequency response of psh to a longitudinal 
frequency sweep is presented in Figure 21. 



Longitudinal Frequency Sweep 

The data presented in the plot shows that the linear model 
response is typical of a second order system with a natural 
frequency of approximately 0.4 rad/sec and with a large amount 
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of damping as evidenced by the 40 decibels per decade decrease 
in gain at frequencies higher than the natural frequency. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

This thesis presents one method for modelling and 
analyzing a helicopter design using Flightlab. The example 
files and model can be modified and used as the basis for 
building models of different helicopters. The Flightlab 
program provides a very good tool for engineering design, 
analysis and simulation of helicopters using nonlinear dynamic 
modeling. The methodical procedure presented herein should 
supplement the user's manuals provided for Flightlab, and 
together they should make future modelling efforts for other 
helicopter designs and types a little easier. The analysis 
procedure shows that the time response of the helicopter to 
standard control inputs using the nonlinear modelling 
capabilities of Flightlab provides more information about the 
aircraft's flight characteristics in that it includes control 
cross coupling and is not limited to the assumption of linear 
modelling. The linear model of the helicopter which is 
extracted from the non-linear model can also be used to 
determine the frequency response to control inputs. This 
guide to using Flightlab for aircraft modelling and analysis 
provides the stepping stone for learning Flightlab and 
creating additional aircraft models for use in control system 
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analysis and additional engineering design at the Naval 
Postgraduate School. 

B. RECOMMENDATIONS FOR FURTHER WORK 

Although the procedures for modelling and analysis 
presented herein provide good results, there are several 
possible areas of improvements. This procedure uses several 
components based on basic theory, e.g., the rigid blade model 
of the rotor system uses blade element theory and the inflow 
model uses momentum theory. It is possible to develop a more 
advanced model of a rotor system using an elastic blade model 
and a better model of the inflow using the genwake theory. 
The capabilities are currently present within the Flightlab 
program. 

The helicopter model used for this work did not include 
any flight control system augmentation so the procedure for 
developing a model of such a system was not discussed. A 
procedural guide for this type of model is also needed. 

During the course of this effort, several problems were 
identified with the Flightlab program itself. Most of these 
were corrected by the developer of the program and implemented 
into this model. However, there was a recent problem 
discovered for which the solution was not included in this 
work. This problem concerns an omission of an angular 
velocity term in the calculation of the lateral acceleration 
rate for the dof6 component. A corrected dof6i component was 
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developed but it uses several different object fields from the 
dof6 component so some modification to the program files 
presented in this work would be necessary. A short comparison 
of the results between models with both types of components 
did not show much change in the aircraft response for the 
tests conducted in this work, but a large effect may occur for 
other tests and this correction should be implemented. 

The final recommendation is to develop a procedural guide 
which discusses how to modify the rigid blade element model to 
enable real time engineering flight simulation of the 
helicopter for use with the pilot's workstation. This 
procedure is necessary because the current version of the 
Flightlab program version being used does not yet take 
advantage of the full capabilities of the parallel processing 
computer systems to run in real time. This procedure would 
involve creating a map of the rotor system states based on 
azimuth, collective position, inflow, and advance ratio. This 
"rotor map" along with replacing the rotor system in the model 
with a rotor map component would create a model capable of 
real time simulation of the helicopter with the pilot's 
workstation. 

The current version of Flightlab includes the programs 
necessary to create a visual scene which uses a generic heads 
up display and a small portable box containing a three axis 
control stick and a collective and/or throttle. Once the 
Flightlab program is updated to run in parallel, this 
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procedure will no longer be necessary. Creating a flight 
simulation of the helicopter still requires the use of a 
multi-processor computer because the Flightlab process 
requires that three programs be run simultaneously, the visual 
scene program, the model program, and the program that is used 
to operate the pilot's workstation. This can only be done on 
a computer system that is capable of sharing memory and 
passing data between these programs in real time, otherwise, 
there would be an update problem between the visual system and 
the simulation input and output. 
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APPENDIX A 


iiiiiiiiitiiiiiiiiiniuiiiiiiiiiiiiiiiiiiiiiiiiiuiuiiiiiiuii 

// file: psh.prolog 
// date: 16 Jul 1992 
// 

// This script loads the data needed for psh model 

iiiiiiinuiiuiiiiiuiiiiiiiiuiiiuiiiiiiiniiiiiiuniiniii 

group data 
// Trim parameters 

trimg * l *Percentage of trim change to apply on a control update"; 

nr « 3 "Number of rotor revolutions between control updates"; 

// Useful Constants 

pi « acos(-l) "Ratio of diameter to circumferance"; 

d2r a pi/180.0 "Degrees to radians conversion factor"; 

r2d - 180.0/pi "Radians to degrees conversion factor”; 

k2f = 6076.115/3600"Knots to feet per second conversion factor”; 

f2k * 3600/6076.115"Feet per second to knots conversion factor"; 

g • 32.2,7/ "Acceleration due to gravity (fpss)"; 

gravity « [0 0 g] "Inertial gravity vector"; 

dt « 0.001 "Integration step size (sec)*; 

eps x s "Solution convergence criteria on the Q' s"; 

imax x 20 "Maximum number of convergence iterations"; 

// Inflow data 

gefl = 0.0625 

ge£2 x l.o 

dwtau x 0.01959 

agl x 90 

chimr x 0 

lam x 0 

nblades = 4 
nseg * 5 

// C G data 

Vweight x 20000 
ixx x 5000.0 
iyy x 40000.0 
izz x 35000.0 
ixy x 0.0 
ixz x 0.0 
iyz x 0.0 

// Rotor data 

lagdamper x 4000 "Lag Damping Coefficient (lbs sec /rad)"; 
imr x 0 "Longitudinal Shaft Tilt + Forward (rad)"; 
mrloc x [-0.5 0 -7.5] "Main Rotor Location (ft)"; 
rpmnom x 21.6667 "Nominal main rotor speed (r/s)”; 
rmr x 30.0 "Main rotor radius (ft)"; 

naz x 24 "Number of azimuth steps/rev"; 

dpsi x 360/naz "Change in azimuth angle/integration step (deg)"; 

dt x d2r*dpsi/rpmnom; // integration step size 


"Total vehicle weight (lbs)"; 

"Total moment of inertia about x (sl-ft2)”; 

"Total moment of inertia about y (sl-ft2)"; 
"Total moment of inertia about z (sl-ft2)"; 
"Total cross product xy (sl-ft2)"; 

"Total cross product xz (sl-ft2)"; 

"Total cross product yz (sl-ft2)"; 


"Cheeseman Bennett ground effect parameter"; 
"Cheeseman Bennett ground effect parameter”; 
"Inflow time constant (sec)"; 

"Alttitude above ground plane (ft)"; 

"Wake skew angle (rad)"; 

"Inflow velocity (nd)"; 

"Number of rotor blades"; 

"Number of blade segments"; 
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tippoa * leroa(3,nblades)‘Inertial positions of blade tips (ft)"; 
tippa ■ zeros (3,1) "Tip path plane angles aO alf blf (rad)"; 


// Load main rotor blade segment aero tables 

// 

// Rotor 

rtablel* read("c_mrotl.tab"); 
clmrl(1:25,1:1) = rtablel(:,1).. 

"Main rotor blade segment cl(alpha,mach) table (nd): high res/low 
angle"; 

cdmrl(1:25,1:1) ■ rtablel(:,2) .. 

"Main rotor blade segment cd(alpha,mach) table (nd): high res/low 
angle"; 

cmmrl(1:25,1:1) « rtablel(:,3) .. 

"Main rotor blade segment cm(alpha,mach) table (nd): high res/low 
angle"; 

clear(rtablel); 

rtable2= read("c_mrot2.tab"); 
clmr2 = rtable2(:,1) 

"Main rotor blade segment cl(alpha) table (nd): low res/high angle"; 
cdmr2 « rtable2(:,2) 

"Main rotor blade segment cd(alpha) table (nd): low res/high angle"; 
cmmr2 = rtable2(:,3) .. 

"Main rotor blade segment cm(alpha) table (nd): low res/high angle"; 
clear(rtable2); 

mrbp = 30.0*d2r .. 

"Main rotor blade segment angle of attack transition angle (rad)"; 
mraoatl = [d2r*[-30 30 2.5] 25] .. 

"Main rotor blade segment angle of attack breakpoint table for high res 
tables"; 

mraoat2 = [d2r*[-180 180 5] 73] .. 

"Main rotor blade segment angle of attack breakpoint table for low res 
tables"; 

mrmacht =[0011] 

"Main rotor blade segment Mach number breakpoint table for high res 
tables"; 

mraal = mraoatl (4) " of rows in main rotor blade segment high res 

tables"; 

mraa2 = mraoat2(4) " of rows in main rotor blade segment low res 
tables"; 

mraml = mrmacht (4) " of cols in main rotor blade segment high res 

tables"; 

exec("blade_seg_geom",l);// Compute the blade segment geometry 
// C G reference Data (DMASS) 

fscg = 296.0 "Fuselage Station of eg (in)"; 

blcg » 0.0 "Buttline Station of eg (in)"; 

wlcg = 113.0 "Waterline Station of eg (in)"; 

xcgf = 0 "Rigid body fuselage station offset (ft)"; 
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ycgf ■ 0 "Rigid body buttline station offset (ft)■; 
scgf « 0 "Rigid body waterline station offset (ft)"; 

cgloc * [xcgf ycgf scgf] "Inertial eg offset (ft)"; 


fmass • (Vweight - nblades*bw)/g"Mass of the fuselage (si)"; 
finertia * [ixx ixy ixz 
ixy iyy iys 

ixz iyz izz] "Inertia matrix of the fuselage (sl-ft2)"; 

// Aero Wing/Fuselage data (ABR03DS and INTERFER) 

fswf * 302.0 "Fuselage station of wing/fuselage (in)*; 

blwf > 0.0 "Buttline station of wing/fuselage (in)"; 

wlwf > X19.0 "Waterline station of wing/fuselage (in)"; 

fswfft * (fscg - fswf)/12.0; 
blwfft - (blcg - blwf)/12.0; 
wlwfft » (wlcg - wlwf)/12.0; 

wfloc«[fswfft blwfft wlwfft] "Aero fuselage location (ft)"; 

exec("c_wfaero.exe",1); // Fuselage Aerodynamics 

wabp x 25.0*d2r "Wing/fuselage angle of attack transition angle 
(rad)"; 

wbbp s 25.0*d2r "Wing/fuselage angle of sideslip transition angle 
(rad)"; 

alfwf • 0 "Wing/fuselage angle of attack (rad)"; 

betwf * 0 "Wing/fuselage sideslip angle (rad)"; 

alfiv x 0 "Angle of attack for the interference effects (rad)"; 

betiv a 0 "Sideslip angle for the interference effects (rad)"; 

wfiv = zeros(3,1) "Interference velocities on the wing/fuselage (fps)"; 
exec("wf_intf.exc",1); // W/F interference 

// Horizontal Stabilizer (AER02D3D and INTERFER) 

fsht x 692.0 "Fuselage station of horizontal stabilizer (in)"; 

blht a 0.0 "Buttline station of horizontal stabilizer (in)"; 

wlht = 95.0 "Waterline station of horizontal stabilizer (in)"; 

fshtft x (fscg - fsht)/12.0; 
blhtft x (blcg - blht)/12.0; 
wlhtft x (wlcg - wlht)/12.0; 

htloc x [fshtft blhtft wlhtft] "Horizontal tail location (ft)"; 

htincang x -0.052 "horizontal incidence angle + up (rad)"; 
htablel x read("c_htaill.tab"); 

clhl(1:13,1:2) x htablel(:,l) .. 

"Horizontal stabilizer cl(alpha,mach) table (nd): high res/low angle"; 
cdhl(1:13,1:2) = htablel{:,2) .. 

"Horizontal stabilizer cd(alpha,mach) table (nd): high res/low angle”; 
cmhl(1:13,1:2) * 0.0*ones(26,1) .. 

"Horizontal stabilizer cm(alpha,mach) table (nd): high res/low angle"; 
clear(htablel); 
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htable2 * read("c_htail2.tab"); 

clh2 * htable2(:, 1) 

"Horizontal stabilizer cl (alpha) 
cdh2 » htable2(:,2) .. 

"Horizontal stabilizer cd(alpha) 
cmh2 * 0.0*ones(19,1) .. 

"Horizontal stabilizer cm(alpha) 

clear(htable2); 

hbp = 30.0*d2r "Horizontal stabilizer angle of attack transition angle 
(rad)«; 

haoatl = [d2r* t-30 30 5] 13] .. 

"Horizontal stabilizer angle of attack breakpoint table for high res 
tables"; 


table (nd): low res/high angle"; 
table (nd): low res/high angle"; 
table (nd): low res/high angle"; 


haoat2 = [d2r*[-90 90 10] 19] .. 

"Horizontal angle of attack breakpoint table for low res tables"; 
hmacht = [0112] 

"Horizontal stabilizer Mach number breakpoint table for high res 
tables"; 


hnal 

= haoatl(4); 

// 

hna2 

a haoat2(4); 

// 

hnml 

a hmacht(4); 

// 


# of rows in clhl,cdhl,cmhl 

# of rows in clh2,cdh2,cmh2 

# of cols in clhl,cdhl,cmhl 


hchord = 2 "Chord length of horizontal stabilizer (ft)"; 
hlen = 9 "Span of horizontal stabilizer (ft)"; 

hdefic = 1 "Lift deficiency factor (nd)"; 

htiv a zeros(3,1) "Interference velocities on the horizontal stabilizer 
(fps)"; 

exec("ht_intf.exc",l); // H tail interference 

// Vertical Tail (AER2D3D and INTERFER) 


fsvt = 716.0; 
blvt a 0.0; 
wlvt a 149.0; 


// Fuselage Station of Vertical Tail 
// Buttline Station of Vertical Tail 

// Waterline Station of Vertical Tail 


fsvtft 
blvtft 
wlvtft 


(fscg - fsvt)/12.0; 
(blcg -blvt)/12.0; 
(wlcg -wlvt)/12.0; 


vtloc a [fsvtft blvtft wlvtft] "Vertical Tail location (ft)"; 
vtablela read("c_vtaill.teUs"); 

clvl(1:13,1:2) = vtablel(:,l) .. 

"Vertical tail cl(alpha,mach) table (nd): high res/low angle"; 
cdvl(1:13,1:2) a vtedalel{:,2) .. 

"Vertical tail cd(alpha,mach) table (nd): high res/low angle"; 
cmvl(1:13,1:2) = 0.0*ones(26,1) .. 

"Vertical tail cm(alpha,mach) table (nd): high res/low angle"; 


clear(vtablel); 


vtable2a read("c_vtail2.tab"); 
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clv2 « vtable2(:,l) 

"Vertical tail cl(alpha) table (nd): low res/high angle"; 
cdv2 > vtable2(:,2) 

"Vertical tail cd(alpha) table (nd): low res/high angle"; 
cmv2 - 0.0*ones(19,1) .. 

"Vertical tail cm(alpha) table (nd): low ree/high angle"; 
clear(vtable2); 


vbp » 30 

0*d2r; 

// 

vaoatl 

m 

[d2r*[-30 30 

5] 

vaoat2 

as 

[d2r*[-90 90 

10] 

vmacht 

as 

[0 112] ; 

// 

vnal 

m 

vaoatl(4); 

// 

vna2 

as 

vaoat2(4); 

// 

vnml 

= 

vmacht(4); 

// 

vchord 

* 

33/7.7; 

// 

vlen 

s 

7.7; 

// 

vdefic 

s 

1; 

// 

vtiv 

a 

zeros(3,1); 


exec ("vt_ 

_intf.exc",1); 

// 


Transition AoA fro clvl to clv2 
13] ;// AoA break pts for clvl,cdvl,cmvl 
19]; // AoA break pts for clv2,cdv2,cmv2 
Mach break pts for clvl,cdvl,cmvl 

# of rows in clvl,cdvl,cmvl 

# of rows in clv2,cdv2,cmv2 

# of cols in clvl,cdvl,cmvl 

Chord length of v tail 
Width of the v tail 
No lift deficiency 


No V tail interference 


// Bailer Tail Trotor (BAILEY) 


fstr = 740.0; 
bltr = 24.0; 
wltr - 185.0; 


// Fuselage Station of tail rotor 
// Buttline Station of tail rotor 

// Waterline Station of tail rotor 


fstrft = (fscg - fstr)/l2.0; 
bltrft * (blcg -bltr)/12.0; 
wltrft * (wlcg -wltr)/12.0; 


trloc = [fstrft bltrft wltrft] "tail rotor location (ft)"; 


atr 

= 

5.73; 


ttr 

= 

1.00; 


cdtr 

= 

0.0; 


dOtr 

S 

0.0087; 

// 

dltr 

= - 

0.0216; 

// 

d2tr 

s 

0.4000; 

// 

biastr 

= 

14.0; 


trblades 

= 3; 


btltr 

s 

0.92; 


bvttr 

s 

1.0; 


bvtltr 

s 

1.0; 


chordtr 

s 

l; 

// 

delttr 

s 

0.001455 

f 

omegatr= 

100; 


rtr 

= 

6.5; 


thettr 

s 

0.0; 


twsttr 

= 

-5; 


td3tr 

s 

-0.5774; 

// 

vbvttr 

as 

30*k2f; 

// 

xibtr 

s 

0.67; 


xlamtr 

= 

1.0; 



// Lift curve slope 
// Tail rotor thrust 
// Rotor head drag coefficient 
Taylor series drag coeff. 

Taylor series drag coeff. 

Taylor series drag coeff. 

// Blade pitch bias 
// Number of blades 
// Blade tip loss 
// Blockage effect parameter 
// Blockage effect parameter 
Blade chord 

// Partial of coning wrt thrust 
// Tail rotor speed 
// Tail rotor radius 
// Tail rotor collective pitch 
// Tail rotor blade twist 
Tan of delta 3 angle 
Blockage velocity parameter 
// Mass moment of inertia 
// Tail rotor inflow (inital value) 


triv = zeros(3,l); 

exec("tr_intf.exc",l);// No Tail rotor interference 
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// Atmosphere data (ATMOS) 

atmtab « read("atmo.tab”) ,-// Get ARDC62 atmosphere tables 

densityt a atmtab(:, 1) // Air density as function of altitude 
ssoundt » atmtab(:,2);// Speed of sound as f(altitude) 
natmo = prod(size(densityt));// Size of data tables 

altt - (0 240000 2000 121]';// Alitude break points 

clear(atmtab); 

wind = zeros(3,1); // No wind 

bodytr * eye(3); 
pos * [0 0 -90]; 
bodyvel * [10 0 0] * ; 


// Data for the UH60 control system 

// 

// Control system data 

mtheta = 17.25*d2r,-// main rotor pitch (rad) 
mthetad =0; // main rotor pitch (rad) 

mthetadd = 0; // main rotor pitch (rad) 


als = -1.I*d2r;// lateral cyclic (rad) 

alsd =0; // lateral cyclic (rad) 

alsdd =0; // lateral cyclic (rad) 

bis = -0.78*d2r;// long cyclic (rad) 

blsd =0; // long cyclic (rad) 

blsdd =0; // long cyclic (rad) 


phase =-4.0*d2r;// Swash plate phase angle 


// Pilot Controls, Trim signals, Fps and Sas 

Xa =0.0; 

Xb =0.0; 

Xc =0.0; 

Xp =0.0; 

Xatrm = 4.95; 

Xbtrm = 1.196; 

Xctrm = 10.817; 

Xptrrn = 1.925; 

// Control bias inputs 

Alsbias = -10.625; 

Blsbias = -10.0; 

ThObias = 1.25; 

Thrbias = 31.875; 

// Swashplate test inputs 

Alschk = 0.0; 

Blschk = 0.0; 

ThOchk = 0.0; 

Thrchk = 0.0; 
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// Swashplate limits 
Alsll - -11.0; 

Alsul * 7.0; 

Blsll - -10.0; 

Blsul * 20.0; 

ThOll . 1.25; 

ThOul - 19.0; 

Thrll - -15.0; 

Thrul - 36.875; 

// Stick to swashplate scale factor 

kxaals » 17.5/9 "Lateral pitch degrees per inch control movement"; 

kxbbls » 30/10 "Longitudinal pitch degrees per inch control 

movement"; 

kxcthO *17.75/12 "Collective pitch degrees per inch control 

movement"; 

kxpthr *-46.875/5.5 "Tail rotor pitch degrees per inch control 
movement"; 


// Pilot input stops 


xcll 

xcul 

xall 

xaul 

xbll 

xbul 

xpll 

xpul 


(ThOll-ThObias)/kxcthO 
(ThOul-ThObias)/kxcthO 
(Alsll-Alsbias)/kxaals 
(Alsul-Alsbias)/kxaals 
(Blsll-Blsbias)/kxbbls 
(Blsul-Blsbias)/kxbbls 
(Thrul-Thrbias)/kxpthr 
(Thrll-Thrbias)/kxpthr 


"Collective stick lower Btop (in)"; 
"Collective stick upper stop (in)"; 
"Lateral cyclic stick lower stop (in)"; 
"Lateral cyclic stick upper stop (in)"; 
"Long. cyclic stick lower stop (in)"; 

"Long. cyclic stick upper stop (in)"; 

"Pedal lower stop (in)"; 

"Pedal upper stop (in)"; 


parentg 

// 

// End of psh.prolog data script 
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iiiiiiiii/iiiiiiiiiniiiiiiiiiiiiiiiiiiiimiiiiiniuiiiiii 

// This is the file c_htaill.tab which contains the high resolution 
// data tables lift and drag coefficients for psh horizontal tail 

iiiiiriiiiiiiiiiiiiuiiiiiiiiniiiniiiiiiniuiiimuiiii 

26 2 


-0.9700 

0.4300 

-1.0300 

0.3700 

-1.0300 

0.3600 

-0.9300 

0.1900 

-0.7100 

0.0400 

-0.3560 

0.0220 

0.0000 

0.0100 

0.3560 

0.0220 

0.7100 

0.0400 

0.9300 

0.1900 

1.0300 

0.3600 

1.0300 

0.3700 

0.9700 

0.4300 

-0.9700 

0.4300 

-1.0300 

0.3700 

-1.0300 

0.3600 

-0.9300 

0.1900 

-0.7100 

0.0400 

-0.3560 

0.0220 

0.0000 

0.0100 

0.3560 

0.0220 

0.7100 

0.0400 

0.9300 

0.1900 

1.0300 

0.3600 

1.0300 

0.3700 

0.9700 

0.4300 


//////////////////////////////////////////////////////////// 

// This is the file c_htail2.tab which contains the low resolution 
// data tables for the horizontal tail for psh 
llllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 

19 2 


0.0000 

1.2000 

-0.2940 

1.1610 

-0.5580 

1.0500 

-0.7450 

0.8880 

-0.8470 

0.7020 

-0.8470 

0.5310 

-0.9700 

0.4300 

-1.0300 

0.3600 

-0.7100 

0.0400 

0.0000 

0.0100 

0.7100 

0.0400 

1.0300 

0.3600 

0.9700 

0.4300 

0.8470 

0.5310 

0.8470 

0.7020 

0.7450 

0.8880 

0.5580 

1.0500 

0.2940 

1.1610 

0.0000 

1.2000 
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//////////////////////////////////////////////////////////// 

// This is the file c_vtaill.tab which contains the high resolution 
// data tables for the vertical tail for psh 

llllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 


26 2 


- 1.0000 

0.3600 

- 1.0000 

0.2650 

-0.9300 

0.1740 

-0.7300 

0.1180 

-0.5000 

0.0660 

-0.2800 

0.0330 

-0.0600 

0.0180 

0.1600 

0.0210 

0.3800 

0.0440 

0.6100 

0.0920 

0.8200 

0.1620 

0.8900 

0.2480 

0.8900 

0.3550 

- 1.0000 

0.3600 

- 1.0000 

0.2650 

-0.9300 

0.1740 

-0.7300 

0.1180 

-0.5000 

0.0660 

-0.2800 

0.0330 

-0.0600 

0.0180 

0.1600 

0.0210 

0.3800 

0.0440 

0.6100 

0.0920 

0.8200 

0.1620 

0.8900 

0.2480 

0.8900 

0.3550 


111111111111111111111111111111111111111111111111111111111111 
// This is the file c_vtail2.tab which contains the low resolution 
// data tables for the vertical tail for psh 
//////////////////////////////////////////////////////////// 


19 2 


0.0000 

1.1000 

-0.1200 

1.0250 

-0.2800 

0.9650 

-0.4600 

0.8750 

-0.6600 

0.7450 

-0.8800 

0.5750 

- 1.0000 

0.3600 

-0.9300 

0.1740 

-0.5000 

0.0660 

-0.0600 

0.0180 

0.3800 

0.0440 

0.8200 

0.1620 

0.8900 

0.3550 

0.8000 

0.5800 

0.6300 

0.7500 

0.4800 

0.8750 

0.3200 

0.9650 

0.1700 

1.0200 

0.0000 

1.0800 
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UlWiW IU "I 


iiiiiiuitiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiniiiimitiiui 

// This is the file c_mrotl.tab which contains the high resolution 
// data tables for the main rotor for psh 

iiiiiiiuiiiiiiiiiiiiiiniiiiiiiniiiiiiiiiiiiiiiiiiiiiiini 

25 3 


-0.9 

0.6 0.6 


-0.75 

0.5 0.56 


-0.7 

0.41 

0.48 

-0.675 

0.35 

0.4 

-0.674 

0.25 

0.34 

-0.8 

0.155 

0.28 

-0.9 

0.1 0.12 


-1.225 

0.05 

0.0000 

-1.12 

0.045 

0.0000 

-0.85 

0.03 

0.0000 

-0.5 

0.0200 

0.0000 

-0.3 

0.02 

0.0000 

0.0 

0.0200 

0.0000 

0.3 

0.02 

0.0000 

0.5 

0.02 

0.0000 

0.8 

0.02 

0.0000 

1.0 

0.02 

0.0000 

1.25 

0.035 

0.0000 

0.985 

0.13 

0.0000 

0.9 

0.2 - .2 


0.7 

0.275 

- .4 

0.745 

0.375 

- .48 

0.8 

0.4025 

- .6 

0.85 

0.52 

- .68 

0.96 

0.602 

- .72 


//////////////////////////////////////////////////////////// 

// This is the file c_mrot2.tab which contains the low resolution 
// data tables for the main rotor for psh 
//////////////////////////////////////////////////////////// 
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0.0 0.02 0.0 


0.6 0.05 -1.0 


0.72 0.125 

1.4 

0.7 0.25 1.2 


0.67 0.3 1.28 


0.8 0.475 1.5 


0.94 0.65 

1.8 

1.025 0.8 2.0 


1.06 1.05 

2.16 

1.03 1.2 2.24 


1.0 1.395 2.32 


0.925 1.525 

2.36 

0.82 1.65 

2.376 

0.725 1.8 2.4 


0.575 1.875 

2.376 

0.425 1.95 

2.36 

0.275 2.02 

2.34 

0.1 2.05 2.32 


-0.75 2.075 

2.3 

-0.23 2.07 

2.2 

-0.4 2.025 

2.12 

-0.55 2.0 2.0 


-0.7 1.93 

1.92 

-0.804 1.83 

1.8 
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-0.975 

1.7 1.7 


-1.06 

1.575 1.6 


-1.14 

1.4 1.4 


-1.16 

1.25 1.2 


-1.15 

1.04 1.04 


-1.06 

0.875 0.8 


-0.9 

0.6 0.6 


-0.7 

0.41 0.48 


-0.674 

0.25 0.34 


-0.9 

0.1 0.12 


-1.12 

0.045 0.0 


-0.5 

0.02 0.0 


0.0 0. 

02 0.0 


0.5 0. 

02 0.0 


1.0 0. 

02 0.0 


0.985 

0.13 0.0 


0.7 0. 

275 -0.4 


0.8 0. 

4025 -0.6 


0.96 

0.602 -0.72 


1.1 0. 

8 -0.8 


1.125 

1.0 -1.0 


1.125 

1.2 -1.2 


1.11 

1.36 -1.3 


1.025 

1.5 -1.48 


0.975 

1.675 -1.6 


0.85 

1.8 -1.7 


0.71 

1.9 -1.8 


0.575 

1.98 -1.9 


0.4 2. 

03 -2.0 


0.23 

2.07 -2.096 


0.75 

2.075 -2.2 


-0.75 

2.07 -2.28 


-0.225 

2.03 -2.3 


-0.4 

1.98 -2.32 


-0.55 

1.9 -2.36 


-0.68 

1.8 -2.376 


-0.79 

1.675 -2.376 


-0.88 

1.55 -2.36 


-0.98 

1.4 -2.32 


-1.03 

1.21 -2.28 


-1.06 

1.075 -2.2 


-1.05 

0.9 -2.12 


-0.975 

0.7 -1.88 


-0.775 

0.5 -1.56 


-0.7 

0.35 -1.3 


-0.75 

0.275 -1.3 


-0.775 

0.175 -1.44 


-0.6 

0.05 -1.2 


0.0 0. 

02 0.0// 
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iiiiuiiiiiuuiiiiiiiiiiniiiiiiniiiiiuiiiiiiiuiiiiiiiu 

// This is the file c_wfaero.exe which contains the aerodynamic 
// data tables for the fuselage for psh 

llllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 

II Wing Fuselage Aero Tables Low Angle High Resolution 

// 

waoatl=[d2r*[-25 25 5] 11]; // Low angle, high resolution alpha brkpts 
waostls[d2r*[-25 25 5] 11]; // Low angle, high resolution beta brkpts 
waoat2«[d2r*[-90 90 10] 19]; // High angle, low resolution alpha 

brkpts 

waost2=[d2r*[-90 90 10] 19]; // High angle, low resolution beta 

brkpts 

waoat3*[d2r*[-90 90 180] 2];// Cross coupling angle of attack brkpts 

wnal = waoatl(4); // Number of points in low angle alpha tables 

wna2 = waoat2(4); // Number of points in high angle alpha tables 

wna3 = waoat3(4); // Number of points in cross couple alpha tables 

wnbl = waostl(4); // Number of points in low angle beta tables 

wnb2 = waost2(4); // Number of points in high angle beta tables 

// Low angle, high resolution Lift, Drag, and Pitch Moment tables as a 
// function of alpha 

clwls[-35.0; -29.0; -22.0; -15.0; -8.0; 

-1.50; 5.0; 12.0; 18.5; 25.0; 

31.5] ; 

cdwl=[30.0; 25.0; 22.5; 20.0; 18.0; 

17.5; 18.5; 21.0; 23.5; 27.0; 

32.5] ; 

cmwl=[-900.0; -770.0; -615.0; -460.0; -305.0; 

-155.0; 0.0; 145.0; 300.0; 450.0; 

590.0] ; 

// Low angle, high resolution Sideforce, Rolling amd Yawing moment 
tables as a 

// function of beta and alpha 

cywl=[87.5; 62.5; 50.0; 31.5; 14.0; 

0.0; -14.0; -31.5; -50.0; -62.5; 

-87.5] ; 

cywl = [cywl cywl ]• 

crwl=[ -95.0; -75.0; -60.0; -37.5; -18.75; 

4.0; 25.0; 43.75; 62.5; 81.25; 

100 . 0 ] ; 

crwl = [crwl crwl]; 

cnwl=[325.0; 275.0; 210.0; 150.0; 80.0; 

0.0; -80.0; -150.0; -210.0; -275.0; 

-325.0]; 

cnwl = [cnwl cnwl]; 


65 









// High angle, low resolution Lift, Drag, and Pitch Moment tables as a 
// function of alpha 

xi= (-25:5:25) ' ; // Independent vars for cmbal 

x={- 9 0:10:9 0)'; // New independent var for cmba2 
clw2 * odli(x,xi,clwl);// 1-D Linear Interpolation 
cdw2 * odli(x,xi,cdwi);// l-D Linear Interpolation 
cmw2 = odli(x,xi,cmwl);// 1-D Linear Interpolation 

// High angle, low resolution Sideforce, Rolling and Yawing moment 
tables 

// as a function of beta and alpha 

cyw2 = odli(x,xi,cywl(:,l));// 1-D Linear Interpolation 
crw2 * odli(x,xi,crwl(:,l));// 1-D Linear Interpolation 
cnw2 = odli(x,xi,cnwl(:,1));// 1-D Linear Interpolation 
cyw2 * [cyw2 cyw2]; 

crw2 = [crw2 crw2]; 

cnw2 * [cnw2 cnw2]; 

// Increments to lift, drag and pitch moment due to sideslip and angle 
of 

// attack 
clbalazeros(11,2); 
cmbalszeros(11,1) ; 
cdbl = zeros(11,1) ; 
cdb2=zeros(19,1); 

// Create a high angle low resolution table from 
// the low angle high resolution table 

xi=(-25:5:25)'; // Independent vars for cmbal,cdbl 

x=(-90:10:90)'; // New independent var for cmba2,cdb2 
cmba2 = odli(x,xi,cmbal);// 1-D Linear Interpolation 
cmbal = [cmbal cmbal] ; 
cmba2 m [cmba2 cmba2]; 

// End of c_wfaero.exe 
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//////////////////////////////////////////////////////////// 

// This is the file mpl.exc which contains the 

// data for the mass distribution along a rotor blade for psh 

// Data from a uniform blade weight distribution analysis 

iiiiiiiiniiiiiiiiiiiiiiiiiiniiiiiniiiiiiniiiiiiiiiiiiiii 


mpl»(0.000 

0.372 

0.050 

0.372 

0.100 

0.372 

0.150 

0.372 

0.200 

0.372 

0.250 

0.372 

0.300 

0.372 

0.350 

0.372 

0.400 

0.372 

0.450 

0.372 

0.500 

0.372 

0.550 

0.372 

0.600 

0.372 

0.650 

0.372 

0.700 

0.372 

0.750 

0.372 

0.800 

0.372 

0.850 

0.372 

0.900 

0.372 

0.950 

0.372 

1.000 

0.372] ; 


//////////////////////////////////////////////////////////// 

// This is the file blade_twist.exe which contains the 

// data for the twist distribution along a rotor blade for pBh 

//////////////////////////////////////////////////////////// 

// Prouty's Sample Helicopter Main Rotor Preset Blade Twist 
// From p648 Helicopter Performance, Stability and Control 


// 

Figure 

10.6 

btw= 

[ 0.000 

0.000 


0.050 

-0.500 


0.100 

- 1.000 


0.150 

-1.500 


0.200 

-2.000 


0.250 

-2.500 


0.300 

-3.000 


0.350 

-3.500 


0.400 

-4.000 


0.450 

-4.500 


0.500 

-5.000 


0.550 

-5.500 


0.600 

-6.000 


0.650 

-6.500 


0.700 

-7.000 


0.750 

-7.500 


0.800 

-8.000 


0.850 

-8.500 


0.900 

-9.000 


0.950 

-9.500 


1.000 

-10.00]; 


67 





iiiiiiiiiiiimiiiiimuiiiiiinnmiiiiniiiiuniiiiiiii 

// This is the file wf_intf.exe which contains the 
// data for the wing fuselage interference effects for psh 

IllllllllllllllllllllinillllllllllllllllUIIIIIIIIIIIIIIUI 

// 

II Dynamic Presure Loss at Wing / Fuselage 

// 

wfaostl = [d2r*[-90 90 180] 2]; 
wfaoatl = [d2r*[-90 90 180] 2]; 

wfqdyn»[ 1 1; 1 1]; // represents no dynamic pressure loss 

wfnsl • wfaostl(4); 
wfnal ■ wfaoatl(4); 

// Fuselage downwash at Fuselage 

wfaost2 - [d2r*[-90 90 180] 2] ; 
wfaoat2 = [d2r*[-90 90 180] 2] ; 

fvzwf m [ 0 0; 0 0]; // no downwash at fuselage caused by fuselage 

w£ns2 s wfaost2(4); 
wfna2 = wfaoat2(4); 

// Fuselage sidewash at Fuselage 

wfaost3 = [d2r*t-90 90 180] 2] ; 
wfaoat3 = [d2r*[-90 90 180] 2]; 

fvywf = [ 0 0; 0 0]; // no sidewash at fuselage caused by fuselage 

wfns3 * wfaost3(4); 
wfna3 * wfaoat3(4); 

// 

// Rotor Wash on Fuselage 

// 

wf chit = [d2r*[0 100 10] 11]; 
wfalft = [d2r*[-6 6 6] 3] ; 
yrfnchi = wfehit (4) ; 
wfnalf = wfalft(4); 

rvxwfi= zeros(ll,3); // interference velocity,Vx due to rotor 

rvywfi = zeros(11,3); // interference velocity,Vy due to rotor 

rvzwfi=-l*ones(11,3); // interference velocity,Vz due to rotor 

II End of wf_intf.exe 
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//////////////////////////////////////////////////////////// 

// This is the file ht_intf.exe which contains the 

// data for the horizontal tail interference effects for psh 

iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiinniiuiuiiiiiiiiiiiiini 

// 

II Dynamic Presure Loss at Horizontal Tail 

// 

htaostl - [d2r*[-90 90 180] 2]; 
htaoatl »* [d2r*[-90 90 180] 2] ; 

htqdyn-[ 1 1; 1 1]; 
htnsl ■ htaostl(4); 
htnal - htaoatl(4); 

// Fuselage downwash at Horizontal Tail 

htaost2 = [d2r*[-90 90 180] 2]; 
htaoat2 * td2r*[-90 90 180] 2]; 

fvzht=[ 0 0; 0 0]; 

htns2 = htaost2(4); 
htna2 = htaoat2(4); 

// Fuselage sidewash at Horizontal Tail 

htaost3 = [d2r*[-90 90 180] 2]; 
htaoat3 = [d2r*[-90 90 180] 2]; 

fvyht*t 0 0; 0 0] ; 

htns3 = htaost3(4); 
htna3 * htaoat3(4); 

// 

II Rotor Wash on Horizontal Tail 

// 

htchit - [d2r*[0 100 10] 11] ; 
htalft = [d2r*[-6 6 6] 3] ; 
htnchi = htchit(4); 
htnalf = htalft(4); 

rvxhti= zeros(ll,3); 

rvyhti = zeros(ll,3); 

rvzhti= zeros(ll,3); 

// End of ht_intf.exe 
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IIIIIIIIUIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIUIIIIIII 

// This is Che file tr_incf.exe which contains Che 
// data for Che tail rotor interference effects for psh 
IIIIIIIIIIIIIIIIIIIIIUIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII 
!/ Dynastic Presure Loss at Tail Rotor is saste as vertical tail 
// 

// Downwash Component of Fuselage Wash on Tail Rotor is sane as 
// horizontal tail 
// 

// Sidewash Component of Rotor Wash on Tail Rotor is same as 
// vertical tail 
// 

// Rotor wash on Tail Rotor is same as 
// horizontal tail 
// End of tr_intf 

llllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 

// This is the file vt_intf.exe which contains the 
// data for the vertical tail interference effects for psh 

llllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 

// 

// Dynamic Presure Loss at Vertical Tail 
vtaostl = [d2r*[-90 90 1801 2] ; 
vtaoatl * [d2r*[-90 90 180] 2] ; 

vtqdyn*[ 1 1; 1 1] ; 
vtnsl * vtaostl(4); 
vtnal * vtaoatl(4); 

// Fuselage downwash at Vertical Tail 

vtaost2 * [d2r*[-90 90 180] 2]; 
vtaoat2 * [d2r*[-90 90 180] 2] ; 

fvzvt=[ 0 0; 0 0] ; 

vtns2 * vtaost2(4); 
vtna2 = vtaoat2(4); 

// Fuselage sidewash at Vertical Tail 

vtaost3 = [d2r*t-90 90 180] 2] ; 
vtaoat3 = [d2r*[-90 90 180] 2]; 

fvyvt=[ 0 0; 0 0] ; 

vtns3 = vtaost3(4); 
vtna3 m vtaoat3(4); 

II 

II Rotor Wash on Vertical Tail 
vtchit * td2r*[0 100 10] 11] ; 
vtalft - [d2r*[-6 6 6] 3] ; 
vtnchi * vtchit(4); 
vtnalf = vtalft(4); 

rvxvtis zeros(ll,3); 

rvyvti = zeros(ll,3); 

rvzvti= zeros(ll,3); 

// End of vt intf.e.tc 
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IlllllllllllllllUHIIIIIIIIIIIIIIIIIItllllllUlllllUUIIII 
// This is the psh.exc file which contains the script 
// file generated from the psh.raod 

IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIU 

GOTO WORLD 
GROUP model 
GROUP Drivetrain 
SINK mrsout; 

NSUMJ mux2; 

SOURCE rra; 

INTTF rai; 

INTTF rsi; 

NSUMJ muxl; 

CONNECT mux2(1) TO mrsout(l); 

CONNECT rra(1) TO raid); 

CONNECT rai(l) TO rsi(l); 

CONNECT rsi(l) TO muxld) ; 

CONNECT rra (1) TO mux2(2); 

CONNECT rai(l) TO muxl(2); 

CONNECT muxl(1) TO mux2(1); 


PARENTG 
GROUP cont 
GROUP Dir 
SINK Thettr; 

GAIN kl; 

GAIN Thettrrad; 

LIMITER Thettrlim; 

SUMJ aj 3; 

SOURCE Thrchk; 

SOURCE Xp; 

SOURCE Xptrm; 

SUMJ aj1; 

SUMJ sj2; 

GAIN Bias; 

CONNECT Thettrrad(1) TO Thettr(1); 
CONNECT Thrchk(1) TO sj3(2); 

CONNECT aj2 d) TO aj3(l); 

CONNECT Bias(1) TO sj2(2); 

CONNECT kl(l) TO sj2(1); 

CONNECT sjl(l) TO kid); 

CONNECT Thettrlimd) TO Thettrradd); 
CONNECT sj3 (1) TO Thettrlimd); 
CONNECT Xp(l) TO ajl(l); 

CONNECT Xptrm(1) TO ajl(2); 


PARENTG 
GROUP Sensors 
SOURCE Beta; 
SOURCE Alpha; 
NGAIN r; 

NGAIN q; 

NGAIN p; 

NGAIN psi; 
NGAIN theta; 
NGAIN Phi; 
GAIN Phideg; 
GAIN Thetadeg; 
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GAIN Rdeg; 

GAIN Qdeg; 

GAIN Pdeg; 

GAIN Psideg; 

SOURCE Unity; 

SOLO S0L0283; 

CONNECT r(1) TO Rdeg(l); 

CONNECT q(l) TO Qdeg(l); 

CONNECT p(l) TO Pdeg{l); 

CONNECT psi(l) TO Psideg(l); 
CONNBCT theta(1) TO Thetadeg(l); 
CONNECT Phi(1) TO Phideg(l); 
CONNECT Rdeg(1) TO S0L0283(1); 


PARHNTG 
GROUP Coll 
SINK mTheta; 

GAIN kl; 

GAIN ThetaOrad; 

LIMITER ThetaOLim; 

SUMJ sj 3; 

SOURCE ThOchk; 

SOURCE Xc; 

SOURCE Xctrm; 

SUMJ sjl; 

SUMJ sj 2; 

GAIN Bias; 

CONNECT ThetaOrad(1) TO mTheta(1); 
CONNECT Bias(1) TO sj2(2); 

CONNECT sjl(l) TO kl(l); 

CONNECT sj2 (1) TO sj3(l); 

CONNECT ThOchk(l) TO sj3(2); 

CONNECT ThetaOLim(l) TO ThetaOrad(1); 
CONNECT sj3(1) TO ThetaOLim(1); 
CONNECT kl(1) TO sj2(l); 

CONNECT Xctrm(l) TO sjl(2); 

CONNECT Xc(l) TO sjl(l); 


PARENTG 
GROUP Lat 
SINK Als; 

GAIN Alsrad; 

LIMITER Alslim; 

SOURCE Alschk; 

SUMJ Sj 3; 

SUMJ sj2; 

GAIN kl; 

SUMJ sjl; 

SOURCE Xatrm; 

SOURCE Xa; 

GAIN Bias; 

CONNECT Alsrad(1) TOAls(l); 
CONNECT sj2(1) TO Sj3(l); 
CONNECT Alschk(1) TO Sj3(2); 
CONNECT Bias{l) TO sj2(2); 
CONNECT kl(l) TO sj2(l); 

CONNECT sjl(1) TO kl(1); 
CONNECT Alslim(1) TO Alsrad (1) ; 
CONNECT Sj3 (1) TO Alslim(1); 
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CONNECT Xatrm(l) TO sjl<2); 
CONNECT Xa(l) TO sjl(l); 


PARENTG 
GROUP Lng 
SINK Bis; 

GAIN Blsrad; 

LIMITER Blslim; 

SUHJ sj 3; 

SOURCE Blschk; 

SUMJ sj 2; 

SOURCE Xb; 

SOURCE Xbtrm; 

SUMJ sj1; 

GAIN k3; 

GAIN Bias; 

CONNECT Blsrad(l) TO Bls(l); 
CONNECT Blslim(1) TO Blsrad(l); 
CONNECT Blschk(1) TO sj3(2); 
CONNECT sj2(1) TO sj3(l); 
CONNECT k3(l) TO sj2(l); 

CONNECT Bias(1) TO sj2(2); 
CONNECT sjl(l) TO k3(1) ; 
CONNECT sj3 (1) TO Blslim(l); 
CONNECT Xb(1) TO sjl(l); 
CONNECT Xbtrm(1) TO sjl(2); 


PARENTG 

CONNECT Sensors_Unity(l) TO Lat_Bias(l); 
CONNECT Sensors_Unity(1) TO Dir_Bias(l); 
CONNECT Sensors_Unity(1) TO Lng_Bias(1); 
CONNECT Sensors_Unity(1) TO Coll_Bias(1); 


PARENTG 
GROUP heli 
ATMOSPHERE atmos; 

GROUP body 
ROTATE htxr2; 

MSENSOR dofsensor; 

TRANSLATE3 wfloc; 

TRANSLATE3 cgloc; 

AER03DS Wf; 

DMASS body; 

DOF6 bodydof; 

AER02D3D Hstab; 

AER02D3D Vstab; 

TRANSLATE3 trboom; 

TRANSLATE3 vtboom; 

ROTATE htzr; 

ROTATE htxr; 

TRANSLATE3 htboom; 

ROTATE vtyr; 

ROTATE vtxr; 

ROTATE trxr; 

BAILEY Trotor; 

CONNECT htxr2(1) TO Hstab(1); 
CONNECT htzr(1) TO htxr2(l); 
CONNECT bodydof(1) TO wfloc(1); 
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CONNECT wfloc(l) TO wf(1) ; 

CONNECT cgloc(l) TO bodyU) ; 
CONNECT bodydof(1) TO cglocU) ; 
CONNECT bodydof(l) TO dofsensor(1) 
CONNECT trxrU) TO Trotor(l); 
CONNECT trboom(l) TO trxr(l); 
CONNECT vtxr (1) TO Vstab(l) ; 
CONNECT vtyr(l) TO vtxr(1); 
CONNECT vtboom(l) TO vtyr(l ); 
CONNECT htboom(l) TO htxr(l); 
CONNECT htxr(l) TO htrr(l); 
CONNECT bodydof(1) TO trboom(1); 
CONNECT bodydof(1) TO htboom(l); 
CONNECT bodydof(1) TO vtboom(l); 


PARENTG 
GROUP rotor 
GROUP bladel 
TRANSLATE biseg6; 

TRANSLATE spar; 

PMASS sparmass; 

TRANSLATE blsegl; 

ROTATE twstl; 

PMASS massl; 

AER02D aerol; 

AER02D aero2; 

PMASS mass2; 

ROTATE twst2; 

TRANSLATE blseg2; 

AER02D aero3; 

PMASS mass3; 

ROTATE twst3; 

TRANSLATE blseg3; 

AER02D aero4; 

PMASS mass4; 

ROTATE twst4; 

TRANSLATE blseg4; 

AER02D aeroS; 

PMASS mass5; 

ROTATE twstS; 

TRANSLATE blseg5; 

MARKPOS tip; 

CHINGE feat; 

TRANSLATE hofl; 

ROTATE rotrl; 

HINGE flap; 

TORSPDM Lead; 

CONNECT Lead(l) TO flap(l); 
CONNECT feat(l) TO spar(l); 
CONNECT spar(1) TO sparmass(1); 
CONNECT blsegl(l) TO massl(1); 
CONNECT blsegl(l) TO twstl(1); 
CONNECT twstl(1) TO aerol(1); 
CONNECT twst2(1) TO aero2(l); 
CONNECT blseg2(l) TO twst2(T); 
CONNECT blseg2(1) TO mass2(l); 
CONNECT twst3(1) TO aero3(l); 
CONNECT blseg3(l) TO twst3(l); 
CONNECT blseg3(l) TO maBs3(l); 
CONNECT twst4(1) TO aero4(l); 



CONNECT blseg4(1) TO twst4(l); 
CONNKCT blseg4(1) TO mass4(1); 
CONNECT tWBtS(l) TO aero5(l); 
CONNECT blsegS<1) TO twst5(l); 
CONNECT blseg5(1) TO raassSU); 
CONNECT spar(1) TO blsegl(1); 
CONNECT blsegl(X) TO blseg2(l) 
CONNECT blseg2(1) TO blseg3(l) 
CONNECT blseg3(1) TO blseg4(X) 
CONNECT blseg4(l) TO blsegS(l) 
CONNECT blsegS(1) TO blsegS(1) 
CONNECT blsegS(1) TO tip(l); 
CONNECT flap(X) TO feat(l); 
CONNECT rotrl(l) TO hofl(l); 
CONNECT hofl(l) TO Lead(l); 


PARENTS 
TPP tppc; 

ROTATE hub; 

TRANSLATE3 mrshaft; 

GROUP blade2 
TORSPDM Lead; 

HINGE flap; 

ROTATE rotrl; 

TRANSLATE hof1; 

CHINGE feat; 

MARKPOS tip; 

TRANSLATE blseg5; 

ROTATE twstS; 

PMASS mass5; 

AER02D aeroS; 

TRANSLATE blseg4; 

ROTATE twst4; 

PMASS mass4; 

AER02D aero4; 

TRANSLATE blseg3; 

ROTATE twst3; 

PMASS mass3; 

AER02D aero3; 

TRANSLATE blseg2; 

ROTATE twst2; 

PMASS mass2; 

AER02D aero2; 

AER02D aerol; 

PMASS massl; 

ROTATE twstl; 

TRANSLATE blsegl; 

PMASS sparmass; 

TRANSLATE spar; 

TRANSLATE blseg6; 

CONNECT Lead(l) TO flap(l); 
CONNECT hof1(1) TO Lead(l); 
CONNECT rotrl(1) TO hofl(l); 
CONNECT flap(1) TO feat(l); 
CONNECT blseg6(X) TO tip{l); 
CONNECT blseg5(1) TO blseg6(1) 
CONNECT blseg4(1) TO blsegS(1) 
CONNECT blBeg3(l) TO blseg4(1) 
CONNECT blseg2(1) TO blseg3(l) 
CONNECT blsegl(1) TO blseg2(l) 
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CONNECT spar(1) TO blsegl(1); 
CONNECT blseg5(l) TO mass5(l); 
CONNECT blseg5(l) TO twst5(l); 
CONNECT twst5(l) TO aero5(l); 
CONNECT blseg4(1) TO mass4(1); 
CONNECT blseg4(1) TO twst4(l); 
CONNECT twst4(1) TO aero4(1); 
CONNECT blseg3(1) TOmass3(l); 
CONNECT blseg3(1) TO twst3(l); 
CONNECT twst3(l) TO aero3(1); 
CONNECT blseg2(1) TOmass2(l); 
CONNECT blseg2(l) TO twst2(l); 
CONNECT twst2(1) TO aero2(l); 
CONNECT twstl(l) TO aerol(1); 
CONNECT blsegl(1) TO twstl(l); 
CONNECT blsegl(l) TO masslU) ; 
CONNECT spar(l) TO sparmass(1); 
CONNECT feat(1) TO spar(l); 


PARENTG 
GROUP blade3 
TRANSLATE b)seg6; 

TRANSLATE spar; 

PMASS sparmass; 

TRANSLATE blsegl; 

ROTATE twstl; 

PMASS massl; 

AER02D aerol; 

AER02D aero2; 

PMASS mass2; 

ROTATE twst2; 

TRANSLATE blseg2; 

AER02D aero3; 

PMASS mass3; 

ROTATE twst3; 

TRANSLATE blseg3; 

AER02D aero4; 

PMASS mass4; 

ROTATE twst4; 

TRANSLATE blseg4; 

AER02D aero5; 

PMASS mass5; 

ROTATE twst5; 

TRANSLATE blseg5; 

MARKPOS tip; 

CHINGE feat; 

TRANSLATE hofl; 

ROTATE rotrl; 

HINGE flap; 

TORSPDM Lead; 

CONNECT Lead(l) TO flap(l); 
CONNECT feat(1) TO spar(l); 
CONNECT spar(l) TO sparmass(1); 
CONNECT blsegl(1) TO massl(1); 
CONNECT blsegl(1) TO twstl(l); 
CONNECT twstl(l) TO aerol(1); 
CONNECT twst2(1) TO aero2(1); 
CONNECT blseg2(l) TO twst2(1); 
CONNECT blseg2(l) TO mass2(1); 
CONNECT twst3(1) TO aero3(l); 
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CONNECT blseg3(X) TO twst3(1); 
CONNECT blseg3(l) TOmass3(l); 
CONNECT twst4(l) TO aero4(l); 
CONNECT blseg4(1) TO twst4(l); 
CONNECT blseg4(1) TO mass4(1); 
CONNECT twst5(1) TO aero5(l); 
CONNECT blseg5(l) TO twst5(l); 
CONNECT blsegS(1) TO mass5(l); 
CONNECT spar(1) TO blsegl(1); 
CONNECT blsegl(1) TOblseg2(l) 
CONNECT blseg2(1) TO blseg3(l) 
CONNECT blseg3(1) TO blseg4(1) 
CONNECT blseg4(1) TO blseg5(l) 
CONNECT blseg5(1) TO blseg6(l) 
CONNECT blseg6(l) TO tip(l); 
CONNECT flap(1) TO feat(l); 
CONNECT rotrl(l) TO hofl(l); 
CONNECT hofl(l) TO Lead(l); 


PARENTG 
GROUP blade4 
TRANSLATE blseg6; 

TRANSLATE Spar; 

PMASS sparmass; 

TRANSLATE blsegl; 

ROTATE twstl; 

PMASS massl; 

AER02D aerol; 

AER02D aero2; 

PMASS mass2; 

ROTATE twst2; 

TRANSLATE blseg2; 

AER02D aero3; 

PMASS mass3; 

ROTATE twst3; 

TRANSLATE blseg3; 

AER02D aero4; 

PMASS mass4; 

ROTATE twst4; 

TRANSLATE blseg4; 

AER02D aero5; 

PMASS mas85; 

ROTATE twst5; 

TRANSLATE blseg5; 

MARKPOS tip; 

CHINGE feat; 

TRANSLATE hofl; 

ROTATE rotrl; 

HINGE flap; 

TORSPDM Lead; 

CONNECT Lead(1) TO flap(l); 
CONNECT feat(1) TO spar(l); 
CONNECT spar(l) TO sparmass(l) 
CONNECT blsegl(1) TO massl(1); 
CONNECT blsegl(1) TO twstl(1); 
CONNECT twstl(1) TO aerol(1); 
CONNECT twst2(1) TO aero2(l); 
CONNECT blseg2(1) TO twst2<l); 
CONNECT blseg2(1) TO mass2(1); 
CONNECT twst3(1) TO aero3(l); 







CONNECT blseg3(1) TO twst3(l); 
CONNECT blseg3(l) TO mass3(1); 
CONNECT twst4(l) TO aero4(1); 
CONNECT blseg4(1) TO twst4(l); 
CONNECT blseg4(1) TO mass4(1); 
CONNECT twst5(1) TO aero5(l); 
CONNECT blsegS(1) TO tvst5{l); 
CONNECT blseg5(1) TO mass5(l); 
CONNECT spar(1) TO blsegl(l); 
CONNECT blsegl(l) TO blseg2U); 
CONNECT blseg2(1) TO blseg3(l); 
CONNECT blseg3(1) TO blseg4(l); 
CONNECT blseg4(1) TO blseg5(l); 
CONNECT blseg5(l) TO blseg6(l); 
CONNECT blsegS(1) TO tip{l); 
CONNECT flap(1) TO feat(l); 
CONNECT rotrl(l) TO hofl(l); 
CONNECT hofl(l) TO LeadU) ; 


PARENTG 

CHIN6E mrspeed; 

SWASHPLATE sp 1,world_data_nblades; 
CONNECT mrspeed(1) TO bladel_rotrl{1); 
CONNECT hub(1) TO tppc(l); 

CONNECT mrshaft(1) TO hub(l); 

CONNECT hub(l) TO mrspeed(1); 

CONNECT mrspeed(1) TO blade2_rotrl(1); 
CONNECT mrspeed(l) TO blade3_rotrl(1); 
CONNECT mrspeed(1) TO blade4_rotrl(1); 
CONNECT sp(3) TO blade3_feat(2); 
CONNECT sp(4) TO blade4_feat(2); 
CONNECT sp(2) TO blade2_feat(2); 
CONNECT sp(l) TO bladel_feat(2); 


PARENTG 

INERTIAL earth; 
GROUP inflow 
INTERFER Wfi; 
INTERFER Hti; 
INTERFER Vti; 
INTERFER Tri; 
UNIFORM!V ul; 


PARENTG 

CONNECT earth(1) TO bodybodydof(1); 

CONNECT body_bodydof(1) fo rotor_mrshaft(1); 
CONNECT rotor_hub<l) TO inflow ul(1); 


PARENTG 

CONNECT Drivetrain_mrsout(1) TO heli_rotor_mrspeed(2) ; 
CONNECT heli_body_dofsensor(2) TO cont_Sensors_Phi(1); 
CONNECT heli_body_dofsensor(2) TO cont_Sensors_theta(1); 
CONNECT heli_body_dofsensor(2) TO cont_Sensors_psi(1); 
CONNECT heli_body_dofsensor(2) TO cont_Sensors_p(1); 
CONNECT heli_body_dofsensor(2) TO cont_Sensors_q(l); 
CONNECT heli_body_dofsensor(2) TO cont_Sensors_r(1); 
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PARBNTG 

ATMOSYSTEM Atmosys; 

AEROSYSTEM Aerosys; 

CONNECT Aerosys(1) TO model_heli_inflow(1) 
CONNECT Aerosys(1) TO model_heli_rotor(1); 
CONNECT Atmosys(1) TO model(1); 


PARBNTG 
GOTO model 
GOTO Drivetrain 

mrsout_N =3; 

mux2_NIl =2; 

mux2_NI2 =1; 

mux2_NO =3 ; 

mux2_SM =[1 0 0;0 1 0;0 0 1] ; 

rra_N =1; 
rra_U =0 ; 

rai_N =1 ; 

rsi_N =1; 

muxl_NIl =1; 

mxixl_NI2 =1; 

muxl_NO =2 ; 
muxl SM = [1 0 ; 0 1] ; 


PARENTG 
GOTO cont 
GOTO Dir 

Thettr_N =1 ; 

kl_N =1; 

kl_K =&world_data_kxpthr; 

The 11 r rad_N = 1 ; 

Thettrrad_K =&world_data_d2r; ' 
Thettrlim_N =1; 

Thettrlim_LL =&world_data_Thrll; 
Thettrlim_UL =&world_data_Thrul; 

s j 3_N =1; 

s j 3_SA1 =1 ; 
s j 3_SA2 =1 ; 

Thrchk_N =1; 

Thrchk_U =&world_data_thrchk; 

Xp_N =1 ; 

Xp_U =&world_data_xp; 

Xptrm_N =1; 

Xptrm_U =&world_data_xptrm; 
s j 1_N =1 ; 
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sjlJSAl -1; 
sjl_SA2 ml; 

sj2_N -1; 
sj 2_SA1 =1; 

8 j 2_SA2 ml; 

Bias_N ml; 

Bias K =&world_data_Thrbias; 


PARENTG 
GOTO Sensors 

Beta_N ml ; 

BetaU *&world_data_betwf; 
Alpha__N mi; 

Alpha_U =&world_data_alfwf; 

r_NO =1 ; 
r_NI =6; 
r_K *[0 0 0 0 0 1]; 

q_NO =1; 
q_NI *6 ; 
q_K = [0 0 0 0 1 0]; 

p_NO =1; 
p_NI =6; 
p_K *[0 0 0 1 0 0]; 

psi_NO =1 ; 
psi~NI =6; 

psiJK *[0 0 1 0 0 0]; 

theta_NO =1; 

theta_NI =6 ; 

theta_K = [0 1 0 0 0 0]; 

Phi_NO *1; 

Phi_NI =6; 

Phi_K =[1000003; 

Phideg_N *1; 

Phideg_K =&world_data_r2d; 

Thetadeg_N =1; 

Thetadeg_K =&world_data_r2d; 

Rdeg_N =1 ; 

Rdeg_K =&world_data_r2d; 
Qdeg_N =1 ; 

Qdeg_K =&world_data_r2d; 
Pdeg_N =1 ; 

Pdeg_K =&world_data_r2d; 

Psideg_N =1 ; 

Psideg_K =&world_data_r2d; 
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Onity_N «1; 

Dnity_D *1; 

S0L0283_N -1; 

SOLO283_ZETA *0.707; 
S0L0283 OMEGA *10; 


PARENTG 
GOTO Coll 

oiTheta N =1 ; 


kl_N *1; 

kl~K =&world data_kxcthO; 


ThetaOrad_N =1; 

ThetaOrad K =&world data d2r; 


ThetaOLim_N =1; 

ThetaOLim_LL =&world_data_Th011; 
ThetaOLim UL =&world data_ThOul; 


sj3 N =1 
sj3~SAl =1 
sj3~SA2 *1 


ThOchk_N =1; 

ThOchk U =&world data thOchk 


Xc_N =1; 

Xc U =&world_data_xc; 


Xctrm_N *1 ; 

Xctrm U =&world_data xctrm; 


sjl_N =1 
s j 1_SA1 =1 
s j 1_SA2 =1 

sj2_N =1 
S j 2_SA1 =1 
s j 2_SA2 =1 


Bias_N =1; 

Bias K =&world data ThObias; 


PARENTG 
GOTO Lat 

A1S_N =1; 

Alsrad_N =1; 

Alsrad K =&world_data_d2r; 


Alslim_N =1; 

Alslim_LL =&world_data_Alsll; 

Alslim UL =&world~data_Alsul; 


Alschk N 


=l; 










Alschk D 


•iworld data alschk; 


Sj3 N 
Sj 3_SA1 
Sj 3_SA2 

8 j 2_N 
sj2_SAl 
sj2_SA2 »1 


kl_N 

kl_K 

sj 1_N 
sj 1_SA1 
sj1 SA2 


■&world data kxaals; 


Xatnn_N =1 ; 

Xatrm“u = &wor1d_data_xatrm; 
Xa_N =1; 

Xa~U =&world_data_xa; 

Bias_N *1 ; 

Bias K =&world data Alsbias; 


PARENTS 
GOTO Lng 

Bis N 


=l; 


Blsrad_N 

Blsrad_K 

Blslim_N 
Blslim_LL 
Blslim UL 


*1 ; 

s&world_data_d2 r; 

=1 ; 

=&world_data_Blsll; 
=&world data Blsul; 


sj 3_N =1 
8 j 3_SA1 *1 
sj3_SA2 =1 

Blschk_N 
Blschk O 


=l; 

=&world data blschk; 


sj2_N =1 
8 j 2_SA1 si 
s j 2_SA2 =1 


Xb_N *1; 

Xb U s&world data xb; 


Xbtrm_N =1; 

Xbtrm O =&world data xbtrm; 


sjl_N si; 
sj 1_SA1 si; 
sjl~SA2 si; 

k3 N *1; 


82 










k3 K 


tworld data kxbbls; 


Bias_N -1; 

Bias K >&world data_Blsbias; 


PARENTG 


PARBNTG 
GOTO heli 

a tmo b _XATMO »&world_datajpos; 
atmos_ALTT *&world_data_altt; 
atmos_NPOINTS =&world_data_natmo; 

atmos_RHOTABLE =&world_data_densityt; 

atmos_SOSTABLE =&world_data ssoundt; 

atmos_GASCN =8,95e+4/28.93; 
atmos_GAMMA = 1.4; 

GOTO body 

htxr2_ANGLE =&world_data_htincang; 
htxr2_AXIS =1; 

dofsensor_NO =6; 
dofsensor_K =zeros(6,18); 

wfloc_LEN =&world_data_wfloc; 

cgloc_LEN =&world_data_cgloc; 

wf_ALPBRK = &world_data_wabp; 
wf_BETBRK =&wor1d_data_wbbp; 
wf_ARGAl =&world_data~waoatl; 

wf_ARGA2 = &wor1d_data_waoat 2; 

wf_ARGA3 =&world~data_waoat3; 

wf_ARGBl =&world_data_waostl; 

wf_ARGB2 =&world_data_waost2; 

wf_NALFAl = &wor1d_data_wna1; 

wf_NALFA2 =&wo r1d_dat a_wna 2; 

wf_NALFA3 =&world_data_wna3; 

wf_NBETAl =&world_data_wnbl; 

wf_NBETA2 =&world_data_wnb2; 

wf_CLAL =&world_data_clwl; 
wf_CDAL =&world_data_cdwl; 
wf_CMAL =&world_data_cmwl; 
wf_CIAH =&world_data_clw2; 
wf_CDAH = &wor1d_data_cdw2; 
wf_CMAH =&world_data_cmw2; 
wf_CLBA =&world_data_clbal; 
wf_CDBL =&world_data_cdbl; 
wf_CMBAL =&world_data_cmbal; 

wf_CNBAL =&world_data_cnwl; 

wf_CRBAL =&world_data_crwl; 

wf_CYBAL =&world_data_cywl; 

wf_CDBH =&world_data_cdb2; 
wf_CMBAH =&wo r1d_dat a_cmb a 2; 

wf_CNBAH =&world_data_cnw2; 

wf_CRBAH =&world_data_crw2; 

wf_CYBAH =&world_data_cyw2; 
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body_MASS »&world_data_fmass; 
body_INERTIA -Aworld_data~finertia; 
body~GRAVITY »&world_data_gravity; 


Hstab_BRKPTS »&world_data_hbp; 
Hstab_CHORD =&vorld_data_hchord; 
Hatab_DEFIC »&world_data_hdefic; 
Hstab_LEN *tworld_data_hlen; 

Hat ab_HAL FA1 »&world_data_hnal; 
Hstab_NALFA2 »&world_data_hna2; 

Hat ab_NMACH »Aworld~data_hnml; 
Hstab_AOATl »&world_data_haoatl; 
Hatab_MACHTl »&world_data_hmacht; 
Hstab_AOAT2 «&world_data_haoat2; 
Hstab_CLl *&world3data_clhl; 
Hstab_CL2 «&world_data_clh2; 
HstabCDl «&world_data_cdhl; 
Hstab_CD2 =&world_data_cdh2; 
Hstab_CMX »&world_data_cmhl; 
Hstab_CM2 *&world_data_cmh2; 

VstaLb_BRKPTS =&world_data_vbp ; 
Vstab_CHORD =&world_data_vchord; 
Vstab_DEFIC »&world_data_vdefic; 
Vstab_LEN *&world_data_vlen; 

Vstab_NALFAl *&world_data_vnal; 

Vstab_NALFA2 =&world_data_vna2; 

Vs tab_NMACH =&world_data_vnml; 

Vstab_AOATl =&world_data_vaoatl; 
Vstab_MACHTl =&world_data_vmacht; 
Vstab_AOAT2 = &wo r1d_data_vaoa12; 
Vstab_CLl =&world_data_clvl; 

Vstab_CL2 =&world_data~clv2; 

V8tab_CDl =&world_data_cdvl; 

V8tab_CD2 =&world_data_cdv2; 

Vstab_CMl =&world_data_cmvl; 

Vstab_CM2 =&world_data_cmv2; 

trboom_LEN =&world_data_trloc; 

vtboom_LEN =&world_data_vtloc; 

htzr_ANGLE =-world_data_pi/2; 

htzr_AXIS =3; 

htxr_ANGLE =world_data_pi ; 
htxr_AXIS =1; 

htboom_LEN =&world_data_htloc; 

vtyr ANGLE =world_data_pi/2; 
vtyr~AXIS =2; 

vtxr_ANGLE =-world_data_pi/2; 

vtxr_AXIS =1; 

trxr_ANGLE =- (90+20)*world_data_d2r; 
trxr_AXIS *1; 

Trotor_BIAS =&world_data_biastr 


84 











TrotorJBLDS 

Trotor~BTL 

Trotor_BVT 

Trotor_BVTl 

Trotor_CHRD 

Trotor_DBLT 

Trotor_OMEG 

Trotor_R 

Trot or~THBTC 

Trotor_TWST 

Trotor_TD3 

Trotor_VBVT 

Trotor_XIB 

Trotor_XLAM 

Trotor_A 

Trotor_CD 

Trotor_DO 

Trotor_Dl 

Trotor_D2 

Trotor_DROT 

Trotor T 


»tworld_data_trblades; 
«tworld~data~btltr; 

«twor1d~data~bvt t r; 
>&«orld_data_bvtltr; 
»&world_data_chordtr; 

■twor1d_data_de1ttr; 
»tworld_data_omegatr; 
■4world_data_rtr; 
»4world_data_thfcttr; 
■&world~data_twsttr; 
*4world”data_td3 tr; 

* &world~data_vbvttr; 
«&world_data_xibtr; 

« &*ro r1d^data_x1amtr; 
■&world~data_atr; 

=&vorld”data_cdtr; 
»4world_data_d0tr; 
*&world”data_dltr; 
*&world”data_d2tr; 

=1 ; 

*&world_data_ttr; 


PARENTG 
GOTO rotor 
GOTO bladel 

blseg6_LEN *&world_data_tipdist; 
spar_LEN =world_data_s1 / 2 ; 


sparmass_MASS =&world_data_sparmass; 

sparmas s_GRAVITY =iwor1d~data_gravity; 


blsegl_LEN =&world_data_seglen(1) ; 


twstl_ANGLE =&world data_segtwst(1); 
twstl_AXIS =1 ; 


massl_MASS =4world_data_segmass(1); 
massl_GRAVITY =&wor1d_data_gravity; 


aerol BRKPTS 

aerol~CHORD 

aerolJDEFIC 

aerol_LEN 

aerol_NALFAl 

aerol_NAI»FA2 

aerol_NMACH 

aerol AQAT1 

aerol~MACHTl 

aerol_AQAT2 

aerol_CLl 

aerol_CL2 

aerol_CDl 

aerol_CD2 

aerol_CMl 

aerol CM2 


=&wor1d_data_mrbp; 
=&world_data_segcor(1) 
=&world_data_segdef(1) 
=&world_data_segwid(1) 
*&world~data_mmal ; 
=&world_data_mma2 ; 
=&world~data_mmml ; 
=&world_data_mraoatl; 
=&world_data_mrmacht; 
*&world~data_mraoat2; 
=&world_data_clmrl; 

=&wor1d~data_clmr2; 
=&world_data_cdmrl; 
=4world_data_cdmr2; 
=£tWorld_data_cmmrl ; 
s&world data cmmr2; 


aero2_BRKPTS =&world_data_mrbp; 
aero2_CH0RD =&world_data_segcor(2) 
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aero2_DEFIC *&world_data_segde£(2); 
aero2~LEN *&vorld_data~segwid(2); 
aero2_NALFAl »&world_data~mrnal ; 
aero2~NALFA2 *&world_data_mma2 ; 
aero2_NMACH «&world_data_mmml ; 
aero2_AQATl s&world_data_mraoatl; 
aero2~MACHTl =&world_data_mrmacht; 
aero22AQAT2 *&world_data_mraoat2; 
aero2_CLl *&world_data_clmrl; 
aero2_CL2 *&world_data_clmr2 ; 
aero2_CDl «&world_data_cdmrl; 
aero2__CD2 »&world_data_cdmr2 ; 
aero2~CMl =iworld_data_c3ranrl; 
aero2~CM2 »&world_data_anmr2; 

mass2_MASS »&world_data_segmass(2); 
mass2~GRAVITY *&world_data_gravity; 

twst2_ANGLE s&world_data_segtw8t(2); 
twst2_AXIS =1 ; 

blseg2_LEN =&world_data_seglen(2); 

ae ro3_BRKPTS =iwor1d_data_mrbp; 
aero3_CH0RD =4world_data_segcor(3); 
aero3_DEFIC »&world_data_segdef(3); 
aero3~LEN =&world_data_segwid(3); 
aero3_NALFAl =&world_data_mrnal ; 
aero3_NALFA2 =&world_dat a_mma 2 ; 
aero3_NMACH *&world_data_mrnml; 
aero3_A0ATl =&world_data_mraoatl; 
aero3_MACHTl =&world_data_mrtnacht ; 
aero3_AQAT2 =&world_data_mraoat2; 
aero3_CLl »&world_data_clmrl; 
aero3_CL2 =&world_data_clmr2; 
aero3_CDl =&world_data_cdmrl; 
aero3_CD2 =4world_data_cdmr2; 
aero3_CMl =&world_data_cmmrl; 
aero3_CM2 =&world_data_cmmr2; 

mass3_MASS =&world_data_segmass(3) ; 
mass3_GRAVITY =&world_data_gravity; 

tws 13_ANGLE =&world_data_segtwst(3); 
twst3_AXIS =1; 

blseg3_LEN =&world_data_seglen(3) ; 

aero4_BRKPTS =&world_data_mrbp; 
aero4_CH0RD =&world_data_segcor(4); 
aero4_DEFIC =&world_data_segdef(4); 
aero4_LEN =&world_data_segwid(4); 
aero4_NALFAl =&world_data_mmal ; 
aero4_NALFA2 =&world_data_mma2 
aero4_NMACH =&world_data_mmml ; 
aero4_A0ATl =&world_data_mraoatl ; 
aero4_MACHTl =&world_data_mrmacht; 
aero4_AQAT2 =&world_data_mraoat2; 
aero4_CLl =&world_data_clmrl; 
aero4_CL2 =&world_data_clmr2; 
aero4 CD1 =&world data cdmrl; . 
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aero4_CD2 «4world_data_cdmr2; 
aero4~CMl «&world~data_aa«nrl ; 
aero4~CM2 ■tworld_data_cnmr2; 

mass4_MASS =&world_data_segmass(4), 
mass4_GRAVITY »&world_data_gravity; 

twst4_ANGLE *&world_data_segtwst(4) ; 
twst4_AXIS *1; 

blseg4_LBN aiworld_data_seglen(4); 

aero5_BRKPTS *&world_data_mrbp; 
aero5_CH0RD B&world_data_segcor(5) ; 
aero5_DBFIC «4world_data_segde£(5) ; 
aero5_LEN =4world_data_segwid(5); 
aero5_NALFAl =&world~data_mrnal; 
aero5_NALFA2 »&world_data_mma2 ; 
aero5_NHACH =&world_data_mmml ; 
aero5_AQATl =&world_data_mraoatl; 
aero5_MACHTl *iworld_data_mrmacht; 
aero5_AOAT2 =4world_data_mraoat2; 
aero5_CLl =&world_data_clmrl; 
aero5_CL2 =4world_data_clmr2; 
aero5_CDl =&world_data_cdmrl; 
aero5_CD2 =&world_data_cdmr2 ; 
aero5_CMl =&world_data_cmmrl; 
aero5_CM2 =&world_data_cmmr2 

mass5_MASS =&world_data_segmass (5); 
mass5_GRAVITY =&world_data__gravity; 

twst5_ANGLE *&world_data_segtwst(5); 
twst5_AXIS «1 ; 

blseg5_LEN =&world_data_seglen(5); 

tip_POS =&world_data_tippos(1:3,1); 

feat_AXIS =1; 

hofl_LEN =&world_data_ho; 

rotrl_ANGLE =0; 

rotrl_AXIS =3; 

flap_AXIS =2; 

Lead_AXIS =3; 

Lead_KK =0; 

Lead CC =&world_data_lagdamper; 

Lead'ANGO =0; 

PARENTG 

tppc_NBLADES =&world_data_nblades; 
tppc_TIPPOS =&world_data_tippos; 

hub_ANGLE =world_data_pi-world_data imr; 

hub_AX!S =2; 

mrshaft LEN =&world_data_mrloc; 
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3; 


GOTO blade2 

Lead_AXIS 
LeadJCX »0; 

Lead_CC *iworld_data_lagdansper ; 
Lead_ANGO =0; 

flap_AXIS »2; 

rotrl_ANGLE xworld_data_pi/2; 
rotrl_AXIS =3; 

hofl_LEN =&wo r1d_dat a_ho; 

feat_AXIS =1; 

tip_POS *&world_data_tippos(1:3,1) ; 

blseg5_LEN =&world_data_seglen(5); 

twst5_ANGLE =&world_data segtwst(5); 
twstS_AXIS =1; 

mass5_MASS =&world_data_segmass(5); 
mass5_GRAVITY =&world_data_gravity; 

aero5_BRKPTS =&world_data_mrbp; 
aero5_CHORD =&world_data_segcor(5); 
aero5_DEFlC =&world_data~segdef(5); 
aero5_LEN =&world_data_segwid(5); 
aero5_NALFAl =&world_data_mmal ; 
aero5_NALFA2 =&world_data_mma2; 
aero5_NMACH =&world_data~mmml ; 
aero5_AOATl =&world~data~mraoatl; 
aero5_MACHTl =&world~data~mrmacht; 
aero5_AOAT2 *&world_data_mraoat2; 
aero5_CLl =&world_data_clmrl; 
aero5_CL2 =&world_data_clmr2; 
aero5_CDl =&world_data_cdmrl; 
aero5_CD2 =&world_data~cdmr2; 
aero5_CMl =&world_data_cnnnrl ; 
aero5_CM2 =4world_data_cmmr2; 

blseg4_LEN =&world_data_seglen(4); 

twst4_ANGLE =&world_data_segtwst(4); 
twst4_AXIS =1; 

mass4_MASS =&world_data_segmass(4); 
mass4_GRAVITY =&world_data_gravity; 

aero4_BRKPTS =&world_data_mrbp; 
aero4_CHORD =&world_data_segcor(4); 
aero4_DEFlC =&world_data~segdef(4) ; 
aero4_LEN =&world_data_segwid(4); 
aero4_NALFAl =&world_data_mmal ; 
aero4_NALFA2 =&wo r1d_data~mma2 ; 
aero4_NMACH =&world_data_mmml ; 
aero4_AQATl =&world_data~mraoatl; 
aero4_MACHTl =&world_data_mrmacht; 
aero4_AOAT2 =&world_data_mraoat2; 
aero4_CLl =&world_data_clmrl; 
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aero4_CL2 
aero4_CDl 
aero4_CD2 
aero4_CMl 
aero4 CM2 


=iworld_data_clmr2 
*4world_data_cdtairl 
»&world_data~cdmr2 
s&world_data_CD8nrl 
world data cmmr2 


blseg3_LEN »&world_data_seglen(3); 

twst3_ANGLE *&world_data segtwst(3); 
twst3_AXIS =1; 


mass3_MASS =&world_data_segmass(3); 
mas s 3 _GRAVITY =4world_data_gravity; 


aero3_BRKPTS 

aero3_CH0RD 

aero3_DEFIC 

aero3_LEN 

aero3_NALFAl 

aero3_NALFA2 

aero3_NMACH 

aero3_AOATl 

aero3_MACHTl 

aero3_AOAT2 

aero3_CLl 

aero3_CL2 

aero3_CDl 

aero3_CD2 

aero3_CMl 

aero3 CM2 


=&wo r1d_data_mrbp; 
»&world_data~segcor(3); 
«4world_data_segdef(3); 
=4world_data~segwid(3); 
=&world_data~mmal ; 
=&world_data_mma2 ; 
=&world_data_mmml ; 
=4world_data_mraoatl; 
=&world_data_mrmacht; 

=&wor1d_data_mraoa12 ; 
=4world_data_clmrl; 
=&world_data_clmr2; 
=4world_data_cdmrl; 
=&world_data_cdmr2; 
=*&world_data_cmmrl; 
=4world_data_cmmr2; 


blseg2_LEN 


=&world_data_seglen(2); 


tws12_ANGLE =iworld_data segtwst(2); 
twst2_AXIS =1 ; 


mass2_MASS =&world_data_segmass(2) ; 
mass2_GRAViTY = &world_data_gravity; 


aero2_BRKPTS 

aero2_CH0RD 

aero2_DEFIC 

aero2_LEN 

aero2_NALFAl 

aero2_NALFA2 

aero2_NMACH 

aero2_AOATl 

aero2_MACHTl 

aero2_AOAT2 

aero2_CLl 

aero2_CL2 

aero2_CDl 

aero2_CD2 

aero2_CMl 

aero2 CM2 


=&world_data_mrbp; 
=&world_data_segcor(2) 
=&world_data~segdef(2) 
=&world_data_segwid(2) 
=&world_data_mrnal 
=&world_data_mrna2 ; 

= £tWorld_data_mrnml ; 
=&world_data_mraoatl; 

=&wo r1d_data_mnnacht; 
=&world_data_mraoat2; 
=&world_data_clmrl; 
=&world_data_clmr2; 
=&world_data_cdmrl; 
=&world_data_cdmr2; 
=&world_data_cmmrl; 
=&world data cmmr2; 


aerol_BRKPTS 
aerol_CHORD 
aerol_DEFIC 
aerol_LEN 
aeroi nalfai 


=&world_data_mrbp; 
=&world_data_segcor(1) 
=&world_data_segdef(1) 
*&world_data_segwid{1) 
=&world_data_mrnal; 
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aerol_NALFA2 

aerol_NMACH 

aerol_AQATl 

aerol_MACHTi 

aerol_A0AT2 

aerol_CLl 

aerol_CL2 

aerol_CDl 

aerol_CD2 

aerol~CMl 

aerol CM2 


»fcworld_data_mma2; 
■&world_data_mraml ; 
»&world_data_mraoatl; 
m&world~data_mrmacht; 
■&world~data~mraoat2; 

=&world~data~clmrl; 
=&world~data_clmr2; 
=&world~data_cdmrl ; 
«&world_data_cdmr2; 
-&world~data_cmmrl ; 
»&world data cmmr2; 


massl_MASS *&world_data_segmass(1); 
massl_GRAVTTY =&world_data_gravity; 

twstl ANGLE s&world data_segtwst(1); 
twstl~AXIS =1; 

blsegl_LEN =&world_data_seglen(1) ; 


sparmass_MASS =&world_data_sparmass 

sparmass J3RAVTTY = &wor 1 d_data__gravi ty ; 

spar_LEN =world_data_sl/2; 

blseg6_LEN «&world_data_tipdist; 


PARENTG 
GOTO blade3 

blseg6_LEN 

spar_LEN 


=&wor1d_data_tipdist; 
=world_data sl/2; 


sparmass_MASS =&world_data_sparmaBS 

spannas s_GRAVITY = &wor 1d_data~gravity; 


blsegl_LEN =&world_data_seglen(l); 


twstl_ANGLE =&world_data_segtwst(1); 
twstl_AXIS =1; 


massl_MASS =&world_data_segmass(1); 
massl_GRAVTTY =&world_data_gravity; 


aero1_BRKPTS 

aerol_CHORD 

aerol_DEFIC 

aerol_LEN 

aerol_NALFAl 

aerol_NALFA2 

aero1_NMACH 

aerol_AOATl 

aerol_MACHTl 

aerol_AOAT2 

aerol_CLl 

aerol_CL2 

aerol_CDl 

aerol_CD2 

aerol_CMl 

aerol CM2 


=&world_data_mrbp; 
=&world~data_segcor(1); 
=&world~data_segdef(1); 
=&world_data_segwid(1); 
=&world__data_mmal ; 
=&world~data_mma2 ; 
=&world_data_mmml ; 
=&world_data_mraoatl; 
=&world~data_mrmacht; 
=&world~data_mraoat2 
=&world~data_clmrl; 
=&world~data_clmr2; 
=&world_data_cdmrl; 
=&world_data_cdmr2; 

=&wo r1d~data_cmmrl; 
=&world~data cmmr2; 
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aero2_BRXPTS «&world_data_mrbp; 
aero2_CH0RD ■&world_data_segcor(2); 
aero2_DBPIC »&world_data_segdef(2); 
aero2_LEN ■&world_data_segwid(2); 
aero2_NALFAl «tworld_data_mmal ; 
aero2_NALFA2 * &world_data_mma2 ; 
aero2_NMACH *&world_data_mrnml; 
aero2_A0ATl «&world_data~mraoatl; 
aero2_MACHTl *&world_data_mrmacht; 
aero2_AOAT2 =&world_data_mraoat2; 
aero2_CLl =&world_data_clmrl; 
aero2_CL2 ■&world_data~clmr2; 
aero2_CDl =&vorld_data_cdmrl; 
aero2_CD2 =&world_data_cdmr2; 
aero2_CMl «&world_data_cmmrl ; 
aero2_CM2 =&vforld_data_cmmr2; 

mass2_MASS =&world_data_segmass(2); 
mass2_GRAVlTY =&world_data_gravity; 

twst2_ANGLE =&world_data_segtwst(2); 
twst2_AXIS =1 ; 

blseg2_LEN =&world_data_seglen(2); 

aero3_BRKPTS =&world_data_mrbp; 
aero3_CH0RD =&world_data_segcor(3); 
aero3_DEFIC =&world_data_segdef(3); 
aero3_LEN =&world_data_segwid(3); 
aero3_NALFAl =&world_data_mmal ; 
aero3_NALFA2 =&world_data_mma2 ; 
aero3_NMACH =&world_data_mrnml ; 
aero3_A0ATl = &wor1d_data_mraoa11; 
aero3_MACHTl =&world~data_mrmacht; 
aero3_AOAT2 =&world_data_mraoat2; 
aero3_CLl =&world_data_clmrl; 
aero3_CL2 =&world_data_clmr2; 
aero3_CDl =&world_data_cdmrl; 
aero3_CD2 =&world_data_cdmr2; 
aero3_CMl =&world_data_cmmrl ; 
aero3_CM2 =&world_data_cmmr2; 

mass3_MASS =&world_data_segmass(3) ; 
mass3_GRAVTTY =&world_data_gravity; 

tws 13_ANGLE =&world_data_segtwst(3); 
twst3_AXIS =1; 

blseg3_LEN =&world_data_seglen(3); 

aero4_BRKPTS =&world_data_mrbp; 
aero4_CH0RD =&world_data_segcor(4); 
aero4_DEFIC =&world_data_segdef(4) ; 
aero4_LEN =&world_data_segwid(4); 
aero4_NALFAl =&world_data_mmaX; 
aero4_NALFA2 =&world_data_mrna2; 
aero4_NMACH =&world_data_mmml ; 
aero4_A0ATl =&world_data_mraoatl; 
aero4_MACHTl =&world_data_mrmacht; 
aero4_AOAT2 =&world_data_mraoat2; 
aero4_CLl =&world_data_clmrl; 
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aero4_CL2 
aero4~CDl 
aero4~CD2 
aero4~CMl 
aero4 CM2 


* 4 world_data_clmr 2 ; 
»&world~data~cdmrl; 
■&world~data_cdmr2; 
■&world~data_cai«nrl ; 
«&world data cmmr2; 


mass4 MASS -&world_data_Begmass(4); 
mass4~GRAVITY -&world_data_gravity; 


twst4 ANGLE =&world data_segtwst(4); 
twst4 AXIS -1; 


blseg4_LEN =&world_data_seglen(4); 


aero5_BRKPTS 

aero5_CH0RD 

aero5_DEFIC 

aero5~LEN 

aero5~NALFAl 

aero5~NALFA2 

ae ro5~NMACH 

aero5~AOATl 

aero5_MACHTl 

aero5_AOAT2 

aero5~CLl 

aero5~CL2 

aero5_CDl 

aero5_CD2 

aero5~CMl 

aero5 CM2 


*&world_data_mrbp; 
■&world_data_segcor(5); 
■&world”data_segdef(5); 
»&world~data~BegwiJ(5); 
*&world_data_mmal ; 
=&world_data_mma2 ; 
=&world_data_mmml ; 
=r < iworld_data_mraoatl ; 
*&world_data_mrmacht; 
*&world~data_mraoat2; 
*&world_data_clmrl; 
*tworld_data_clmr2; 
«&world_data~cdmrl; 
=&world_data_cdmr2; 
*&world_data_csmnrl ; 
*&world data cmmr2; 


mass5_MASS *&world_data_segmass(5); 
mass5 GRAVITY * &vor1d_data_gravity; 


twst5_ANGLE =&world_data_segtwst(5) ; 
twst5 AXIS =1; 


bl8eg5_LEN =&world_data_seglen(5); 


tip_POS =&world_data_tippoe(1:3,1); 


feat_AXIS =1; 

hofl LEN =&world data_ho; 


rotrl_ANGLE =world data_pi; 
rotrl AXIS =3; 


flap_AXIS =2; 

Lead_AXIS =3; 

Lead_KK =0; 

Lead_CC *=&world data lagdamper; 
Lead ANGO «0~ 


PARENTG 
GOTO blade4 

blseg6_LEN =&world_data_tipdist; 
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sparJLEN *world_data_sl/2; 

spannass_MASS *&world_data_8parmass; 

spannass_GRAVITY »fcworld_data_gravity; 

blsegl_LEN =&world_data_seglen(1); 

tws 11_ANGLE » &wor1d_data_segtwst(1) ; 
twstl_AXIS -1; 

massl_MASS »&world_data_segmass(1) ; 
massl_GRAVITY =&world_data_gravity; 

aerolJBRKPTS »&world_data_mrbp; 
aerol_CHORD ■&world_data_segcor(1) ; 
aerol_DEFIC »&world_data_segdef(l) ; 
aerol_LEN =&world_data_segwid(1); 
aerol_NALFAl »&world_data_mmal ; 
aerol_NALFA2 =&world_data_mma2 ; 
aerol_NMACH *&world_data_mmml; 
aerol_AOATl «iworld_data_mraoatX; 
aerol_MACHTl «&world_data_mrmacht; 
aerol_AQAT2 =&world_data_mraoat2; 
aeroX_CLl =&world_data_clmrl; 
aerol_CL2 =&world_data__clmr2 ; 
aerol_CDl *&world_data_cdmrl; 
aerol_CD2 =&world_data_cdmr2; 
aerol_CMl =&world_data_cmmrl; 
aerol_CM2 =&world_data_cmmr2; 

aero2_BRKPTS =&world_data_mrbp; 
aero2_CH0RD =&world_data_segcor(2); 
aero2_DEFIC =&world_data_segdef(2); 
aero2_LEN =&world_data_segwid(2); 
aero2_NALFAl =&world_dat a_mma 1 ; 
aero2_NALFA2 =&world_data_mma2 ; 
aero2_NMACH =&world_data_mrnml; 
aero2_AQATl =&world_data_mraoatl; 
aero2_MACHTl =&world_data_mrmacht; 
aero2_AOAT2 =&world_data_mraoat2; 
aero2_CLl =&world_data_clmrl; 
aero2_CL2 =&world_data_clmr2; 
aero2_CDl =&world_data_cdmrX; 
aero2_CD2 =&world_data_cdmr2; 
aero2_CMl =&world_data_cmmrl; 
aero2_CM2 =&world_data_cmmr2; 

mass2_MASS =&world_data_segmass(2); 
mass2_GRAVITY =&world_data_gravity; 

twst2_ANGLE =&world_data_segtwst(2); 
twst2_AXIS =1 ; 

bl8eg2_LEN =&world_data_seglen(2); 

aero3_BRKPTS =&world_data__mrbp; 
aero3_CH0RD =&world_data_segcor(3); 
aero3_DEFIC =&world_data_segdef(3); 
aero3_LEN =&world_data_segwid(3); 
aero3_NALFAl =&world_data_mmal ; 
aero3_NALFA2 =&world_data_mrna2; 
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aero3_NMACH 

aero3_AOATl 

a«ro3~MACHTl 

aero3_AOAT2 

aero3~CLl 

aero3~CL2 

aero3~CDl 

aero3~CD2 

aero3~CMl 

aero3 CM2 


■&world_data_mnnnl ; 

« fcwor 1 d_dat a~mraoat 1 ; 
■&wor1d_daca~mrmacht; 
«&world~data~mraoat2; 
■&world~data_clmrl; 
■Evorld~data_clmr2; 
*&world_data_cdmrl; 
=&world_data~cdmr2; 
*&vorld_data~awnrl; 

* {.world data cmmr2; 


mass3_MASS *&world_data_segmass(3); 
mass3_GRAVTTY «&world_data_gravity 


twst3_ANGLE «&world_data segtwst(3); 
tWSt3_AXIS »1; 


blseg3__LEN 


*&world_data_seglen(3); 


aero4_BRKPTS 

aero4_CHORD 

aero4_DEF!C 

aero4_LEN 

aero4_NALFAl 

aero4~NALFA2 

aero4_NMACH 

aero4_AOATl 

aero4_MACHTl 

aero4_A0AT2 

aero4_CLl 

aero4_CL2 

aero4_CDl 

aero4_CD2 

aero4_CMl 

aero4 CM2 


=&world_data_mrbp ; 
*&world_data_segcor(4); 
=&world_data~segdef(4); 
=&world_data__segwid(4) ; 
=&world_data~mmal ; 
=&world_data_mma2 ; 

*iwor 1 d_dat a_mmml ; 
=&world_data~mraoatl; 
=&world_data~mrmacht; 
=&world_data~mraoat2; 
=&world_data_clmrl; 
=&world_data_clmr2; 
=&world_data_cdmrl; 
={.world~data~cdmr2 ; 
=&world~data_cmmrl; 

=&world data~canmr2; 


mass4_MASS =&world_data_segmass (4) ; 
mass4_GRAVTTY = &wor1d_data_gravi ty 


twst4_ANGLE =&world_data_segtwst(4); 
twst4_AXIS =1; 


blseg4_LEN =&world_data_seglen(4 ); 


aero5_BRKPTS 

aero5_CH0RD 

aero5_DEFIC 

aero5 LEN 

aeroSJNALFAl 

aero5_NALFA2 

aero5_NMACH 

aero5~A0ATl 

aeroS~MACHTl 

aero5~A0AT2 

aero5~CLl 

aero5~CL2 

aero5~CDl 

aero5~CD2 

aero5~CMl 

aero5~CM2 


=&world_data_mrbp; 
=&world_data~segcor(5) 
=&world_data_segdef(5) 
=&world_data_segwid(5) 
=&world_data~mmal ; 
=&world_data~mma2 ; 
=&world_dat a~mmm 1 ; 
=&world_data_mraoatl; 
=&world_data~mrmacht; 

=&wo r1d_data_mraoa 1 2; 
=&world_data~clmrl; 

«&wor1d~data~clmr2; 
=&world~data~cdmrl; 
=&world~data_cdmr2; 
=&world~data_anmrl; 
»&world~data cmmr2; 
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mass5_MASS «&world_data_segmass(5); 
mass5_GRAVTTY »&world_data_gravity; 

twst5_ANGLE B&world data_segtwst(5); 
twst5_AXIS -1; 

blseg5_LEN *&world_data_seglen(5); 

tip_POS *&world_data_tippos(1:3,1); 

feat_AXIS «1; 

hof1_LEN «&world_data_ho; 

rotrl_ANGLE =3*world_data_pi/2; 
rotrl_AXIS *3; 

flap_AXIS =2; 

Lead_AXIS =3; 

Lead_KK =0; 

Lead_CC =&world_data_lagdampe r; 
Lead_ANG0 =0; 


PARENTG 

mrspeed_AXIS =3; 

sp_PSI =&mrspeed_u(1); 
sp_PSID =&mrspeed_u(2); 
sp_PSIDD «&mrspeed_u(3) ; 

sp_AlS =&world_data_als 
sp_AlSD =&world_data~alsd; 
sp_AlSDD *&world_data_alsdd; 

sp_BlS =&world_data_bls; 
sp_BlSD =&world_data_blsd; 
sp_BlSDD =&world_data_blsdd; 

sp_NBLADES =&world_data_nblades; 

spJPHASE =&wor1d_data_phase; 

sp_THE =&world_data_mtheta; 
sp_THED =&world_data_mthetad; 
sp_THEDD =&world_data_mthetadd; 


PARENTG 
GOTO inflow 

Wfi_ALPHA =&world_data_alfiv; 

Wfi_ARGl =&world_data_wfaostl; 

Wfi_ARG2 =&world_data_wfaoatl; 

Wfi_ARG3 =&world_data_wfaost2; 

Wfi_ARG4 =4world_data_wfaoat2 ; 

Wfi_ARG5 =&world_data_wfaost3; 

Wfi_ARG6 *&world_data_wfaoat3 ; 

Wfi_ARG7 =&world_data_wfchit ; 

Wfi_ARG8 «&world_data_wfalft; 

Wfi_AlF =&world_data_tippa(2); 
Wfi_BETA = &wor1d_data_be t i v; 

Wfi_CHI *&world_data_chimr; 
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Iff i_LAM «6world_data_lam ; 

Wf i_OMBGA ■ 4 world_data_rpmncm; 

Wfi_RMR *&world_data~rmr ;~ 
Wfi_ROWl -&world”data_wfnsl; 

Wfi_COLl »&world~data~vfnal; 

W£i_ROW2 »&world~data_wfns2; 

Wfi COL2 ■&world”data_wfna2; 

Wfi_ROW3 »&world“data~wfns3; 

Wfi_COL3 ■&world~data~wfna3; 

Wfi~ROW4 »&world”data_wfnchi; 

Wfi_COL4 «&world”data~wfnalf; 

Wfi_ROWS ■&wor 1d~data~wfnchi; 

Wfi~COL5 *&world”data_wfnalf; 

Wfi~ROW6 «&world2data~wfnchi; 

Wfi_COL6 siworld~data~wfnalf; 

Wf i_TABLEl = &wo r 1 d”da t a~ wf qdyn ; 

Wfi_TABLE2 *&world~data_fvzwf; 

Wfi~TABLE3 =&world_data~fvywf; 

Wf i_TABLE4 *&world”data_rvxwfi; 

Wf i_TABLE5 «=&world_data_rvywf i ; 

Wfi_TABLE6 *&world~data_rvzwfi; 

Wfi_TRANS =&world_data_bodytr; 

WfiVELIN = &wo r1d~databodyve1 

WfiVELOUT =&world”data_wfiv; 

Hti_ALPHA =&world_data_alfiv; 

Hti_ARGl =&world_data~htaostl 

Hti_ARG2 =&world~data_htaoatl 

Hti_ARG3 =&world_data~htaost2 

Hti_ARG4 =&world_data_htaoat2 

Hti_ARG5 =&world“data_htaost3 

Hti_ARG6 =&world~data htaoat3 

Hti_ARG7 =&world_data_htchit; 

Hti_ARG8 =&world_data_htalft; 

Hti_AlF =&world_data~tippa(2); 
Hti_BETA =&world_data_betiv; 

Hti_CHI =&world_data_chimr; 
Hti_LAM =&world_data_lam; 
Hti_OMEGA =&world_data_rpmnom; 

Hti_RMR =&world_data_rmr; 

Hti_ROWl = &world_data_htns1; 

Hti_COLl =&world_data_htnal; 

Hti_ROW2 =&world~data_htns2; 

Hti_COL2 =&world~data~htna2; 

Hti_ROW3 =&world_data_htns3; 

Hti_COL3 = &wo r1d_data_htna3; 

Hti_ROW4 =&world~data_htnchi; 

Hti_COL4 =&world~data~htnalf; 

Hti_ROW5 =&world”data~htnchi; 

Hti_COL5 = &wor1d”data~htnalf; 

Hti_ROW6 *&world_data_htnchi; 

Hti_COL6 =&world2data_htnalf; 

Hti_TABLEl =&world~data~htqdyn; 

Hti_TABLE2 =&world”data~fvzht; 

Hti_TABLE3 =&world~data~fvyht; 

Hti_TABLE4 =&world~data_rvxhti; 

Hti_TABLE5 ■&world”data~rvyhti; 

Hti_TABLE6 »&world_data~rvzhti ; 

Hti~TRANS =fitWorld~data~bodytr; 

Hti_VELIN =&world~data_bodyvel 

Hti”VELOUT *&world data htiv; 






Vti_ALPHA «4world_data_alfiv; 

Vti~ARGl «&world_data~vtaostl; 

Vti_ARG2 »iworld_data~vtaoatl; 

Vti~ARG3 » &wo r1d_data_vtaos 12; 

Vti_ARG4 « &wo r1d_data_vtaoa12; 

Vti_ARG5 a&world_data_vtaost3; 

Vti_ARG6 «=&world_data_vtaoat3 ; 

Vti_ARG7 «&world_data~htchit; 

Vti_ARG8 aiworld~data_htalft ; 

Vti_AlF =&world_data_tippa(2); 

Vt i~BETA »&world_data_betiv; 

Vti_CHI «&world_data_chimr; 
vti_IiAM »&world_data_lam; 
Vti_OMEGA aiworld_data_rpmnooi; 

Vti_RMR =&world_data_rmr; 

Vti_ROWl -&world_data_vtns1; 

Vti~COLl a&world_data_vtnal; 

Vti_R0W2 a&world_data_vtns2; 

Vti_COL2 =&world_data_vtna2; 

Vti_ROW3 «&world_data_vtns3; 

Vti_COL3 «4world_data_vtna3; 

Vti_R0W4 =&world_data_htnchi; 

Vti_COL4 =iworld_data_htnalf; 

vti_R0W5 a&world_data_htnchi; 

Vti_COL5 =&world_data_htnalf; 

Vti_ROW6 =&world_data_htnchi; 

Vti_C0L6 =&world_data_htnalf; 

Vti_TABLEl =&world_data_vtqdyn; 

Vti_TABLE2 =&world_data_fvzvt; 

Vti_TABLE3 =&world_data_fvyvt; 

Vti_TABLE4 =&world_data_rvxhti; 

Vt i_TABLE5 aiwor 1d_data_rvyhti; 

Vti_TABLE6 =&world_data_rvzhti; 

vti_TRANS =&world_data_bodytr; 

Vti_VEI<IN =&world_data_bodyvel ; 

Vti_VELOUT =&world_data_vtiv; 

Tri_ALPHA =&world_data_alfiv; 

Tri_ARGl =&world_data_vtaost1; 

Tri_ARG2 =&world_data_vtaoatl; 

Tri_ARG3 =&world_data_htaost2; 

Tri_ARG4 =&world_data_htaoat2; 

Tri_ARG5 =&world_data_vtaost3; 

Tri_ARG6 = &world_data_vtaoat3; 

Tri_ARG7 a&world_data_htchit; 

Tri_ARG8 =&world_data_htalft; 

Tri_AlF a&world_data_tippa(2); 
Tri_BETA = &wor1d_data_betiv; 

Tri_CHI a&world_data_chimr; 
Tri_LAM =&world_data_lam; 
Tri_OMEGA a&world_data_rpmnom ; 

Tri_RMR =&world_data_rmr; 

Tri_ROWl =&world_data_vtnsl; 

Tri_COLl =&world_data_vtnal; 

Tri_ROW2 =&world_data_htns2; 

Tri_COL2 =&world_data_htna2; 

Tri_R0W3 =&world_data_vtns3; 

Tri_C0L3 =&world_data_vtna3; 

Tri_ROW4 a&world_data_htnchi; 

Tri_C0L4 =&world_data_htnalf; 

Tri R0W5 a&world data htnchi; 
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Tri COLS 
Tri”ROW6 
Tri”COL6 
Tri TABLE1 
Tri”TABLE2 
Tri~TABLE3 
Tri TABLE4 
Tri~TABLE5 
Tri TABLE6 
Tri”TRANS 
Tri_VELIN 
Tri - VELOUT 


«&world_data_htnalf; 

-&wor1d”data_htnchi; 
«&world”data_htnalf; 

-iworld_data_vtqdyn; 
-iworld_data_fvsht; 
-fcworld”data_fvyvt; 
*4world”data_rvxhti; 
-&world_data_rvyhti; 
-&world”data_rvrhti; 
»&world”data_bodytr; 
*&world_data_bodyvel ; 
={.world data triv; 


ul_ALT = &world_data_agl; 
ul DWTAU =&world_data_dwtau ; 

ul“GBFl *iworld_data_gef l~ 
ul_GEF2 «&world”data_gef2; 
uljNBLADES » & wor 1 d ~dat a_nb lades ; 

ul_NSEG = &wo r1d_data_nse g~ 
ul_RPMNOM «&world_data_rpmnom; 

ul RMR =&world_data_rmr; 


PARENTG 

PARENTG 

PARENTG 
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Illlllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 

// file: psh.prolog 
//’ date: 16 Jul 1992 

// 

// This script loads the data needed for psh model 

lllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 

group data 
// Trim parameters 

trimg = l "Percentage of trim change to apply on a control update"; 

nr =3 "Number of rotor revolutions between control updates"; 

// Useful Constants 

pi = acos(-l) "Ratio of diameter to circumferance"; 

d2r = pi/180.0 "Degrees to radians conversion factor"; 

r2d = 180.0/pi "Radians to degrees conversion factor"; 

k2f = 6076.115/3600"Knots to feet per second conversion factor"; 

f2k = 3600/6076.115"Feet per second to knots conversion factor"; 

g =32.2;// "Acceleration due to gravity (fpss)"; 

gravity * [0 0 g] "Inertial gravity vector"; 

dt = 0.001 "Integration step size (sec)"; 

eps = 5 "Solution convergence criteria on the Q's"; 

imax = 20 "Maximun number of convergence iterations"; 

// Inflow data 

gefl = 0.0625 "Cheeseman Bennett ground effect parameter"; 

gef2 « 1.0 "Cheeseman Bennett ground effect parameter"; 

dwtau = 0.01959 "Inflow time constant (sec)"; 

agl = 90 "Alttitude above ground plane (ft)"; 

chimr = 0 "Wake skew angle (rad)"; 

lam = 0 "Inflow velocity (nd)"; 

nblades = 4 "Number of rotor blades"; 

nseg = 5 "Number of blade segments"; 
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// c G data 

Weight ■ 20000 
ixx * 5000.0 
iyy * 40000.0 
izz ■ 35000.0 
ixy * 0.0 
ixz a 0.0 
iyz a 0.0 

// Rotor data 

lagdamper = 4000 "Lag Damping Coefficient (lbs sec /rad)"; 
imr = 0 "Longitudinal Shaft Tilt + Forward (rad)"; 
mrloc = [-0.5 0 -7.5] "Main Rotor Location (ft)"; 
rpmnom = 21.6667 "Nominal main rotor speed (r/s)"; 
rmr = 30.0 "Main rotor radius (ft)"; 

naz a 24 "Number of azimuth steps/rev"; 

dpsi = 360/naz "Change in azimuth angle/integration step (deg)"; 

dt = d2r*dpsi/rpmnom; // Integration step size 

tippos = zeros(3,nblades)"Inertial positions of blade tips (ft)"; 
tippa = zeros (3,1) "Tip path plane angles aO alf blf (rad)"; 


"Total vehicle weight (lbs)"; 

"Total moment of inertia about x (sl-ft2)"; 

"Total moment of inertia about y (sl-ft2)*; 
"Total moment of inertia about z (sl-ft2)"; 
"Total cross product xy (sl-ft2)*; 

"Total cross product xz (sl-ft2)"; 

"Total cross product yz (sl-ft2)"; 


// Load main rotor blade segment aero tables 

// 


// Rotor 

rtablel= read("c_mrotl.tab"); 


Clmrl (1:25,1:1) = 
"Main rotor blade 
angle"; 

cdmrl(1:25,1:1) = 
"Main rotor blade 
angle"; 

cmmrl(1:25,1:1) = 
"Main rotor blade 
angle"; 


rtablel(:,1) . . 

segment cl(alpha,mach) table (nd): high res/low 
rtablel(:,2) . . 

segment cd(alpha,mach) table (nd): high res/low 
rtablel(:,3) .. 

segment cm(alpha,mach) table (nd): high res/low 


clear(rtablel); 


rtable2= read("c_mrot2 . tab") ,- 
clmr2 = rtable2(:,l) .. 

"Main rotor blade segment cl(alpha) table (nd): low res/high angle"; 
cdmr2 = rtable2(:,2) .. 

"Main rotor blade segment cd(alpha) table (nd): low res/high angle"; 
cmmr2 = rtable2(:,3) .. 

"Main rotor blade segment cm(alpha) table (nd): low res/high angle"; 
clear(rtable2); 


mrbp = 30.0*d2r .. 

"Main rotor blade segment angle of attack transition angle (rad)"; 
mraoatl = [d2r*[-30 30 2.5] 25] .. 

"Main rotor blade segment angle of attack breakpoint table for high res 
tables"; 

mraoat2 = [d2r*[-180 180 5] 73] .. 
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"Main rotor blade segment angle of attack breakpoint table for low res 
tables"; 

mrmacht =[0011] 

"Main rotor blade segment Mach number breakpoint table for high res 
tables"; 

mrnal = mraoatl(4) " of rows in main rotor blade segment high res 
tables"; 

mma2 = mraoat2(4) " of rows in main rotor blade segment low res 
tables"; 

mmml = mrmacht (4) " of cols in main rotor blade segment high res 

tables"; 

exec("blade_seg_geom",1);// Compute the blade segment geometry 
// C G reference Data (DMASS) 

fscg = 296.0 "Fuselage Station of eg (in)"; 

blcg = 0.0 "Buttline Station of eg (in)"; 

wlcg = 113.0 "Waterline Station of eg (in)"; 

xcgf = 0 "Rigid body fuselage station offset (ft)"; 

yegf = 0 "Rigid body buttline station offset (ft)"; 

zcgf = 0 "Rigid body waterline station offset (ft)"; 

cgloc = [xcgf yegf zcgf] "Inertial eg offset (ft)"; 


fmass = (Vweight - nblades*bw)/g"Mass of the fuselage (si)"; 
finertia = [ixx ixy ixz 
ixy iyy iyz 

ixz iyz izz] "Inertia matrix of the fuselage (sl-ft2)"; 

// Aero Wing/Fuselage data (AER03DS and INTERFER) 

fswf = 302.0 "Fuselage station of wing/fuselage (in)"; 

blwf = 0.0 "Buttline station of wing/fuselage (in)"; 

wlwf = 119.0 "Waterline station of wing/fuselage (in)"; 

fswfft = (fscg - fswf)/12.0; 
blwfft = (blcg - blwf)/12.0; 
wlwfft = (wlcg - wlwf)/12.0; 

wfloc=[fswfft blwfft wlwfft] "Aero fuselage location (ft)"; 

exec("c_wfaero.exc",1); // Fuselage Aerodynamics 

wabp = 25.0*d2r "Wing/fuselage angle of attack transition angle 
(rad)"; 

wbbp = 25.0*d2r "Wing/fuselage angle of sideslip transition angle 
(rad)"; 

alfwf = 0 "Wing/fuselage angle of attack (rad)"; 

betwf = 0 "Wing/fuselage sideslip angle (rad)"; 

alfiv = 0 "Angle of attack for the interference effects (rad)"; 

betiv = 0 "Sideslip angle for the interference effects (rad)"; 

wfiv = zeros(3,l) "Interference velocities on the wing/fuselage (fps)"; 
exec("wf_intf.exc",1); // W/F interference 

// Horizontal Stabilizer (AER02D3D and INTERFER) 
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fsht • 692.0 "Fuselage station of horizontal stabiliser (in)*; 

blht » 0.0 "Buttline station of horizontal stabilizer (in)"; 

wlht >95.0 "Waterline station of horizontal stabilizer (in)"; 

fshtft = (fscg - fsht)/12.0; 
blhtft = (blcg - blht)/l2.0; 
wlhtft * (wlcg - wlht)/l2.0; 

htloc = [fshtft blhtft wlhtft] "Horizontal tail location (ft)*; 

htincang * -0.052 "horizontal incidence angle + up (rad)"; 
htablel * read("c_htaill.tab"); 

clhl(1:13,1:2) = htablel(:,1) .. 

"Horizontal stabilizer cl(alpha,mach) table (nd): high res/low angle"; 
Cdhl(1:13,1:2) = htablel(:,2) .. 

"Horizontal stabilizer cd(alpha,mach) table (nd): high res/low angle"; 
cmhl(1:13,1:2) = 0.0*ones(26,1) .. 

"Horizontal stabilizer cm(alpha,mach) table (nd): high res/low angle"; 
clear(htablel); 

htable2 = read("c_htail2.tab"); 
clh2 = htable2(:,1) .. 

"Horizontal stabilizer cl (alpha) table (nd) : low res/high angle"; 
cdh2 = htable2(:,2) .. 

"Horizontal stabilizer cd(alpha) table (nd): low res/high angle"; 
cmh2 = 0.0*ones(l9,l) .. 

"Horizontal stabilizer cm(alpha) table (nd): low res/high angle"; 
clear(htable2); 

hbp = 30.0*d2r "Horizontal stabilizer angle of attack transition angle 
(rad)"; 

haoatl = [d2r*[-30 30 5] 13] .. 

"Horizontal stabilizer angle of attack breakpoint table for high res 
tables"; 

haoat2 = [d2r*[-90 90 10] 19] .. 

"Horizontal angle of attack breakpoint table for low res tables"; 
hmacht =[0112] .. 

"Horizontal stabilizer Mach number breakpoint table for high res 
tables"; 

hnal = haoatl (4); // # of rows in clhl, cdhl, cmhl 

hna2 = haoat2(4); // # of rows in clh2, cdh2, cmh2 

hnml = hmacht (4); // # of cols in clhl, cdhl, cmhl 

hchord = 2 "Chord length of horizontal stabilizer (ft)"; 
hlen = 9 "Spam of horizontal stabilizer (ft)"; 

hdefic = 1 "Lift deficiency factor (nd)"; 

htiv = zeros(3,1) "Interference velocities on the horizontal stabilizer 
(fps)"; 

exec("ht_intf.exc",1); // H tail interference 
// Vertical Tail (AER2D3D and INTERFER) 
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fsvt * 716.0; // Fuselage Station of Vertical Tail 

blvt =0.0; // Buttline Station of Vertical Tail 

wlvt = 149.0; // Waterline Station of Vertical Tail 

fsvtft = (fscg - fsvt)/12.0; 
blvtft = (blcg -blvt)/12.0; 
wlvtft = (wlcg -wlvt)/12.0; 

vtloc = [fsvtft blvtft wlvtft] "Vertical Tail location (ft)"; 
vtablel= read("c_vtaill.tab"); 

clvl(1:13,1:2) = vtablel(:,l) .. 

"Vertical tail cl(alpha,mach) table (nd): high res/low angle"; 
cdvl(1:13,1:2) = vtablel(:,2) .. 

"Vertical tail cd(alpha,mach) table (nd): high res/low angle"; 
cmvl(1:13,1:2) * 0.0*ones(26,1) .. 

"Vertical tail cm(alpha,mach) table (nd): high res/low angle"; 
clear(vtablel); 

vtable2= read("c_vtail2.tab"); 
clv2 = vtable2(:,l) .. 

"Vertical tail cl(alpha) table (nd): low res/high angle"; 
cdv2 = vtable2(:,2) 

"Vertical tail cd(alpha) table (nd): low res/high angle"; 
cmv2 = 0.0*ones(19,1) .. 

"Vertical tail cm(alpha) table (nd): low res/high angle"; 
clear(vtable2); 

vbp = 30.0*d2r; // Transition AoA fro clvl to clv2 

vaoatl = [d2r*[-30 30 5] 13);// AoA break pts for clvl,cdvl,cmvl 

vaoat2 * [d2r*[-90 90 10] 19]; // AoA break pts for clv2,cdv2, cmv2 

vmacht = [0112]; // Mach break pts for clvl, cdvl, cmv’ 

vnal = vaoatl (4); // # of rows in clvl, cdvl, cmvl 

vna2 = vaoat2(4); // # of rows in clv2, cdv2, cmv2 

vnml = vmacht (4); // # of cols in clvl, cdvl, cmvl 

vchord =33/7.7; // Chord length of v tail 

vlen =7.7; // Width of the v tail 

vdefic =1; // No lift deficiency 

vtiv = zeros(3,l); 

exec("vt_intf.exc" , 1); // No V tail interference 

// Bailer Tail Trotor (BAILEY) 

fstr = 740.0; // Fuselage Station of tail rotor 

bltr =24.0; // Buttline Station of tail rotor 

wltr = 185.0; // Waterline Station of tail rotor 

fstrft = (fscg - fstr)/12.0; 
bltrft = (blcg -bltr)/l2.0; 
wltrft = (wlcg -wltr)/l2.0; 

trloc = [fstrft bltrft wltrft] "tail rotor location (ft)"; 

atr =5.73; // Lift curve slope 

ttr =1.00; // Tail rotor thrust 

cdtr =0.0; // Rotor head drag coefficient 
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dOtr > 0.0087; // Taylor aeries drag coeff. 

dltr >-0.0216; // Taylor aeries drag coeff. 

d2tr = 0.4000; // Taylor series drag coeff. 

biastr >14.0; // Blade pitch bias 

trblades >3; // Number of blades 

btltr > 0.92; // Blade tip loss 

bvttr * 1.0; // Blockage effect parameter 

bvtltr >1.0; // Blockage effect parameter 

chordtr> 1; // Blade chord 

delttr > 0.001455; // Partial of coning wrt thrust 

omegatrs 100; // Tail rotor speed 

rtr =6.5; // Tail rotor radius 

thettr =0.0; // Tail rotor collective pitch 

twsttr = -5; // Tail rotor blade twist 

td3tr = -0.5774; // Tam of delta 3 amgle 

vbvttr = 30*k2f; // Blockage velocity parameter 

xibtr =0.67; // Mass moment of inertia 

xlamtr >1.0; // Tail rotor inflow (inital value) 

triv = zeros(3,1); 

exec("tr_intf.exc",l);// No Tail rotor interference 
// Atmosphere data (ATMOS) 

atmtab = readCatmo.tab") ;// Get ARDC62 atmosphere tables 

densityt = atmtab(:,l);// Air density as function of altitude 
ssoundt = atmtab(:,2);// Speed of sound as f(altitude) 
natmo = prod(size(densityt));// Size of data tables 
altt = [0 240000 2000 121]';// Alitude break points 

clear(atmtab); 

wind = zeros(3,1); //No wind 

bodytr a eye(3); 
pos =[0 0 -90] ; 
bodyvel = [10 0 0] ' ; 


// Data for the UH60 control system 

// 

// Control system data 

mtheta = 17.25*d2r;// main rotor pitch (rad) 
mthetad =0; // main rotor pitch (rad) 

mthetadd = 0; // main rotor pitch (rad) 

als = -1.I*d2r;// lateral cyclic (rad) 

alsd =0; // lateral cyclic (rad) 

alsdd =0; // lateral cyclic (rad) 

bis = -0.78*d2r;// long cyclic (rad) 

blsd =0; // long cyclic (rad) 

blsdd =0; // long cyclic (rad) 

phase =-4.0*d2r;// Swash plate phase angle 


// Pilot Controls, Trim Signals, Fps and Sas 

Xa =0.0; 

Xb =0.0; 
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XC -0.0; 

Xp -0.0; 

Xatrm - 4.95; 

Xbtrm - 1.196; 

Xctnn = 10.817; 

Xptrm = 1.925; 

// Control bias inputs 

Alsbias - -10.625; 

Blsbias - -10.0; 

ThObias = 1.25; 

Thrbias - 31.875; 

// Swashplate test inputs 

Alschk - 0.0; 

Blschk * 0.0; 

ThOchk = 0.0; 

Thrchk = 0.0; 

// Swashplate limits 
Alsll = -11.0; 

Alsul = 7.0; 

Blsll = -10.0; 

Blsul = 20.0; 

ThOll = 1.25; 

ThOul - 19.0; 

Thrll - -15.0; 

Thrul = 36-875; 

// Stick to swashplate scale factor 

kxaals = 17.5/9 "Lateral pitch degrees per inch control movement"; 

kxbbls = 30/10 "Longitudinal pitch degrees per inch control 

movement"; 

kxcthO -17.75/12 "Collective pitch degrees per inch control 
movement"; 

kxpthr =-46.875/5.5 "Tail rotor pitch degrees per inch control 
movement"; 


// Pilot input stops 

xcll = (ThOll-ThObias)/kxcthO 
xcul = (ThOul-ThObias)/kxcthO 
xall = (Alsll-Alsbias)/kxaals 
xaul = (Alsul-Alsbias)/kxaals 
xbll = (Blsll-Blsbias)/kxbbls 
xbul = (Blsul-Blsbias)/kxbbls 
xpll = (Thrul-Thrbias)/kxpthr 
xpul = (Thrll-Thrbias)/kxpthr 


"Collective stick lower stop (in)"; 
"Collective stick upper stop (in)"; 
"Lateral cyclic stick lower stop (in)" 
"Lateral cyclic stick upper stop (in)" 
"Long. cyclic stick lower stop (in)" 

"Long. cyclic stick upper stop (in)" 

"Pedal lower stop (in)”; 

"Pedal upper stop (in)"; 


parentg 

// 

// End of psh.prolog data script 
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APPENDIX B 


iiiiiiiiiiiiiiiiiiiiiiiiiiiniiiiiiiiiiiiuiiiuuuiiiiiiiiiininiiiii 

// file: psh.def // 

// Definition file for psh with 4 rigid blades and // 

// five segments and the primary flight control system // 

iiiiiiiiiiiiiiiiiiiiiiiiHiiiiiiuiiiiiiiiiniiiiiiiiiiiiiiniiiiiiniii 

II Set path to directories where files are located 

path("analysis/") 

path("functions/") 

path("data-tables/"); 

path("operations/") 

path("system/") 

path("solution/") 

path("rotors/") 

// load the functions used 

exec("cycle.fun",1); // Load the cycle function 

exec("odli.fun",1); // Load the Id linearinzation function 

// Model construction: 

goto world; // Go to the top group of the scope program 

exec ("psh.prolog", l) // Load the data file for the psh 

exec("psh.exc",l) // Load the model file for the psh 

exec("psh.epilog",l) // equivalence data for solution and system 

components 

exec("system.exc",1); // Create and connect the system 

component 

exec("solution.exc",l); // Setup solution 

init; // Links, equivalences, and analyzes data-flow 

exec("mbc.exc",1); // Setup a MBC transformation 

exec("configure.exc",1) // Configure for reporting and display 

world_data_omega = 650/30 // Set main rotor speed 

init; // Relink since we equivalenced fields in 

configure.exc 

world:'.setup // One time internal initialization and method setup 

world::reset // Set i.c. for states/derivatives and invokes model 

exec ("assemble.exc"),- // setup the model for running 


//// end of psh.def //////// 





iiiiiiiuiiiiiiiiiiiiiiiuniiiiiiiiiiiiiiiiniiiiimiiiiiiiiiiiiiiiiii 

// file: psh.epilog // 

// Epilog file for psh rigid blade simulation // 

iiiiiiiiiiiuiiuiiiiiiiiiiiiiiiihiiiiiiiiiiiiuiiiiiiiiiihiiiiiiiiiii 


world_data_alfwf 
world_data_betwf 
world_data_wfiv 
world~data~triv 
wor1d_data~vtiv 
world data htiv 


&world_model_he1i_body_wf_alpha; 
iworld_model_he1i”body~wf~beta; 
iworld_model_he1i”body~wf_dwash; 
&world_model~heli”body~trotor_dwash; 
&world_model_he1i”body_vstab_dwash; 
iwor1d_model_he1i~body~hstab_dwash; 


world_data_bodytr 
world_data_bodyve1 
world~data_pos 
world_data_pos « [0 


« &world_model_heli_body_bodydof_cfrt; 

• &world_model_heli_body_bodydof~cfrl(4:6) 

• &world_model~he1i_body_bodydof_xi(1:3); 

0 -90]; 


world_data_chimr ■ iworld_model_heli_inf1ow_ul_chi; 

world_data_lam > &world_model_he1i_inf1ow_ul_x; 


world_data_tippa 


iwor1d_model_he1i_rotor_tppc_tippp; 


world_model_drivetrain_rai_xic * world_data_rpmnom; 
world_model_heli_rotor_mrspeed_u(2)=world_data_rpmnom; 

world_model_heli_body_bodydof_xiic(l:3) * world_data_pos; 

world_model_heli_body_bodydof_xicf=ones(6,1); 
world_model_heli_body_bodydof_xcf-ones(6,1); 
world_model_heli_body_bodydof_xdcf=ones(6,1); 

world_model_heli_inflow_ul_xic = 0.05; 
world_model_heli_inflow_ul~x * 0.05; 

// Control connection 


init(model_cont); 

world_data_alfiv 
world_data_betiv 

world_data_bls 
world_data_als 
world_data_mtheta 
world data thettr 


&world_model_cont_sensors_alpha_y; 
iworld_model_cont_sensors_beta_y; 

&world_model_cont_lng_Bls_y; 

&world_model_cont_lat_Als_y; 
&world_model_cont_coll_mtheta_y; 
&world_model_cont_dir_Thettr_y; 


// Motion sensor gain matrix 

world_model_heli_body_dofsensor_k(1,10)*1 
world_model_heli_body_dofsensor_k(2, 11)=1 
world_model_heli_body_dofsensor_k(3,12)=1 
world~model_heli_body_dofsensor~k(4,13)»1 
world_model~heli_body_dofsensor_k(5,14)=1 
world_model_heli_body_dofsensor_k(6,15)*1 

// End of uh60.epilog 
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//////////////////////////////////////////////////////////////////////// 

// file: 8yetem.exe // 

// This script creates and connects the system component // 

iniinininiiuuiiiiiiiiiiiiiiiininiuiiiiiiiiiiiiiiniimiiiiiiu 

goto world 

system sys 
connect sys to world 

init ; 

// End of system.exe script 

///////////////////////////////////////////////////////////////// 

// file: solution.exe // 

// This file creates the solution configuration for psh // 

// date: 19 July 1993 // 

llllllllllllllllllllllllllimilllIHUIIIIIIIIIIIIIIIIIIIIIIIIII 

pushg(world); 


HSOLVE helisolve; 

helisolve_dt world_data_dt; 

helisolve_delta = l/2**14; 

helisolve_epsilon=& world_data_eps; 
helisolve_maxiter=& world_data_imax; 
helisolve_reassemble= 1; 

CSOLVE drivesolve; 

drivesolve_dt world_data_dt; 

drivesolve_delta= 1/2**14; 
drivesolve_epsilon=& world_data_eps; 
drivesolve_maxiter=& world_data_imax; 
drivesolve_reassemble= 0; 
drivesolve_ycf = 0; 

CSOLVE contsolve 

contsolve_dt =& world_data_dt; 

contsolve_delta = 1/2**14 ; 

contsolve~epsi1on=& world_data_eps; 
conteolve~maxiter=& world_data_imax; 
contsolve_reassemble= 0; 
contsolve__ycf = 0 ; 

HSOLVE rotorsolve; 

rotorsolve_dt =& world_data_dt; 
rotorsolve_delta= 1/2**14; 
rotorsolve_epsilon=& world_data_eps; 
rotorsolve_maxiter«& world_data_imax; 
rotorsolve_reassetnble= 1; 

HSOLVE rhsolve; 

rhsolve_dt »& world_data_dt; 

rhsolve_delta = 1/2**14; 

rhsolve_epsilon =& world_data_eps; 
rhsolve_maxiter =& world”data_imax; 
rhsolve_reassemble= 1; 
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PSOLVE topsolve; 

topsolve_dt =4 world_data_dt; 
topsolve_delta * ~l/2**14; 
topsolve_epsilon=& world_data_eps; 
topsolve_maxiter*4 world_data_imax; 
topsolve_reassemble= 0; 
topsolve_ycf » 0 ; 

// Low-level solution groups: 

connect contsolve to model_cont; 

connect contsolve to model~heli_body_dofsensor; 

connect drivesolve to model_drivetrain; 

connect helisolve to model_heli; 

connect rotorsolve to model_heli_rotor; 
connect rotorsolve to model_heli_inflow_ul; 

// Mid-level: RHSOLVE combines ROTOR and HELI groups 
connect rhsolve to helisolve; 
connect rhsolve to rotorsolve; 

// Top-level: 

connect topsolve to drivesolve; 
connect topsolve to rhsolve; 
connect topsolve to contsolve; 

// Default architecture: 

// Iterate between RH, DRIVE, and CONT solution groups: 
topsolve_cf = 1; 
rhsolve_cf * 0; 

// These must always be zero, as they are bottom-level groups: 

drivesolve_cf = 0 ; 

contsolve_cf = 0; 

helisolve_cf = 0; 

rotorsolve_cf = 0 ; 


popg 

exec("load_operation8.exe",1) 

// End of solution.exe 

/////////////////////////////////////////////////////////////// 

// file: load_operations.exe // 

// This file loads the simulation opns that propagate the model// 
// date: 22 July 1993 // 

/////////////////////////////////////////////////////////////// 

// Load the simulation operations 

exec("step.op",l) 
exec("steprotor.op",1) 
exec("stepnc.op",1) 

// End of load_operations.exe 
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//////////////////////////////////////////////// 

// File: step.op 
// Author: Joe Bnglish 

// Created: 3 May 1993 

///////////////////////////////////////////////// 

operation step 

// 

// STEP operation for helicopter model: 

// Propagate forward one time step and solve for the new states, 

// using CONTSOLVE, DRIVESOLVE, and RHSOLVE. 

// 

// Integrate everything: 

[contsolve::integ, rhsolve::integ, drivesolve::integ],- 

// Drivetrain and control system outputs: 

[drivesolve::geny, contsolve::geny]; 

// propagate the control system again bo the nonlinear components 
converge: 

world_contsolve::genq; 
world_contsolve::geny; 

// Rotor / helicopter: 

[rhsolve:;genq, rhsolve::geny] 

// convergence iterations on the helicopter: 
repeat 2, 
rhsolve::solve; 
rhsolve::genq; 
end 

end 

// End of step.op 

/////////////////////////////////////////// 

// File: steprotor.op 
// Author: Joe English 

// Created: 3 May 1993 
//////////////////////////////////////////// 
operation steprotor 
// 

// STEPROTOR operation for helicopter 

// Runs the rotor alone (just drivetrain and rotorsolve) 

// Iterates on ROTORSOLVE twice 
// Integrate everything: 

[drivesolve::integ, rotorsolve::integ]; 

[drivesolve::geny]; 

[helisolve::genq]; 

// Rotor / helicopter: 

[rotorsolve::genq, rotorsolve::geny] 

// two convergence iterations on the rotor is enough 
repeat 2, 

rotorsolve::solve; 
rotorsolve::genq; 
end 

end 

// End of steprotor.op 
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/////////////////////////////// 

// Pile: stepnc.op 
// Author: Joe English 

// Created: 3 May 1993 

////////////////////////////// 

operation stepnc 

// 

// STEPNC operation for helicopter model: 

// Propagate everything but the control system 
// using DRIVESOLVE and RHSOLVE. 

// 

// Integrate everything: 

[rhsolve::integ, drivesolve::integ]; 

// Drivetrain outputs: 

[drivesolve::geny, contsolve::geny]; 

// Rotor / helicopter: 

[rhsolve::genq, rhsolve::geny] 

// convergence iterations on the helicopter: 
repeat 2, 
rhsolve::solve; 
rhsolve::genq; 
end 

end 


// End of stepnc.op 

////////////////////////////////////////////////////////// 

// file: mbc.exc 
// date: 21 July 1993 

// This script sets up a multi blade coordinate 
// transformation of the rotor states to improve the speed and 
// accuracy of a rotating blade states for the 
// psh rigid blade element model 
///////////////////////////////////////////////////////// 
pushg(world_model_heli_rotor); 

describe("MBC: Rigid blade element") 

nblades =& world_data_nblades; 
ncomp =2; // flapping, lead-lag 

ndof = 1; 


mbcxfrm mbc 1,nblades*ncomp; 


j=i; 

for 


end 


mbc_nblades 
mbc_ncomp = 

mbc_ndof = 

mbcjpsi = 

mbc_omega = 

mbc_mbcf = 

i=l:nblades 
connect mbc(j) 
connect mbc(j) 


&nblades; 
ncomp; 
ndof ; 

&mrspeed_u(1); 

&mrspeed_u(2); 

0; // leave off initially 


to blade$i_flap; j=j+l; 
to blade$i_lead; j=j+l; 


popg; 

// End of mbc.exc 
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iiiiiiiiiiiiiiihiiiiiihiiiiiuiiiiiiiiiiihuiiihi 

// file: configure.exe 
// date: 2 September 1993 

// This file configures the parameter group by selectively 
// addressing relevant model and data information 

iiiiiiiiiiiiiiiiiiiiniiiiiiiiiiiiiiniiiiiiiiiiiiui 

pushg(world); group results; popg 
pushg(world); group cpg; 


g « &world_data_g; 
d2r * &world~data_d2r; 
r2d » iwor1d_data_r2 d; 
k2f » &world_data~k2f; 
f2k = &world_data_f2k; 
pi » &world_data_pi; 
naz = world_data_naz; 

// Rigid body data 

fmass = &world_data_fmass; 

finertia = &world_data_finertia; 

rbpfrl = &world_model_heli_body_body_pfrl(1:3); 

rblf = &world_model~heli_body_body_pfo; 

trb = &wor1d_model_he1i_body_body_pfrt; 

// Aircraft Center of Gravity 

pos = &world_model_heli_body_bodydof_xi(1:3); 
ean = &world_model_heli_body_bodydof_xi(4:6); 
ivel = &world_model_heli_body_bodydof_xid(l:?); 
eand = &world_model_heli_body_bodydof_xid(4:6); 
bve1 * iworld_model_he1i_body_bodydof~xd(1:3); 
brat * &world_model_heli_body_bodydof_xd(4:6); 
bacc a &world_model_heli_body_bodydof_xdd; 
bata a &world_model_heli_body_bodydof_xdd(l:3); 
bara a &world_model_heli_body_bodydof_xdd(4:6); 
tran a &world_model_heli_body_bodydof_cfrt; 

ixx = &world_data_ixx; 
iyy = &world_data_iyy; 
izz a &world_data_izz; 
ixz a &world_data_ixz; 
ixy a &world_data_ixy; 
iyz a &world_data_iyz; 

// Atmospheric data 

tamb a &world_model_heli_atmos_tamb; 
rho a &world_model_heli_atmos_rho; 
rhoO a world_data_densityt(1) ; 

// Aircraft station data 

mrloc a &world_data_mrloc; 
cgloc a &world_data3cgloc; 
wfloc a &world_data^wfloc; 
htloc a &wor1d_data”h11oc; 
trloc a &world_data_trloc; 
vtloc a &world data vtloc; 
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mrpfrl 

wfpfrl 

htpfrl 

vtpfrl 

trpfrl 

cgpfrl 

rbpfrl 


iwo r1d_model_he1i_rotor_hub_cf r1(1:3); 
&world_model_heli_body_wf jpf rl (1:3) ; 
&world_model~heli_body_hstab_pfrl(1:3); 
&world_model_heli_bodyjvstab_pfrl(1:3); 
&world_model_heli_body_trotor_pfrl(1:3); 
&world_model_he1i_body~bodydof_cf r1(1:3); 
&world_model_heli_body_bodyjpfrl(1:3); 


// Data for inflow and tip path plane data 

tippa * &world_data_tippa; 
chimr = &world~data_chimr; 
dwash = &world_model_heli_inflow_ul_x; 
mu = &world_model_heli_inflow_ul_mu; 

eta = &world_model_heli_inflow_ul_ct; 

gndef = &world__model_heli_inf low_ul_gndef ; 
ahubf = &world2model_heli_inflow_ul~faa(1:3); 
ahubm = 4world_model_heli_inflow_ul_faa(4:6); 

psimr = &world__model_heli_rotor_mrspeed_u (1) 
omega = &world_model_heli_rotor_mrspeed_u(2) 
omegad = &world_model_heli_rotor_mrspeed_u(3) 

hublf = &world__model_heli_rotor_hub_cfo; 
thubf = &world_model_heli_rotor_hub_cfo(1:3); 
thubm = &world_model_heli_rotor_hvib_cfo(4 :6) ; 
hbtr = &world_model_heli_rotor_hub_cfrt; 


// Fuselage aerodyncunic data 


wf cv 

wf iv 

wf tv 

wfalfa 

wfbeta 

wfqdyn 

edwf 

cywf 

clwf 

erwf 

cmwf 

enwf 

wf If 

wf for 

wfmom 

wf if 

wfifor 

wfimom 

twf 


&world_model_heli_body_wf_pfrl(4:6); 
&wor1d_data_wfiv; 

&wor1d_model_he1i_body_wf_ve1; 
&world_model_heli_body_wf_alpha; 

&world_model_he1i_body_wf_beta; 

&wor1d_model_he1i_body_wf_qdyn; 
&world_model_he1i_body_wf_cd; 

&worid_model_he1i_body_wf_cy; 
&world_model_heli_body_wf_cl; 

&wo r1d_mode1_he1i_body_wf_cr; 
&world_model_heli_body_wf_cm; 
&world_model_heli_body_wf_cn; 
&world_model_heli_body_wf jpfo; 
&world_model_heli_body_wf_pfo(1:3); 
&world_model_he1i_body_wf_pfo(4:6); 
&wor1d_model_he1i_body_wf_faero; 
&world_model_heli_body~wf_faero(1:3); 
&world_model_heli_body_wf_faero(4:6); 
StWorld_model_heli_body_wf_pf rt ; 
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// Horizontal tail aerodynamic data 


htcv 

htiv 

httv 

htalfa 

htbeta 

htqdyn 

cdht 

clht 

cmht 

htlf 

htfor 

htmom 

htif 

htifor 

htimom 

tht 


& wor 1d_mode1_he1i_body_hstab_p fr1(4:6); 
&world~data_htiv; _ 
&world~modeT_heliJbody_hstab_vel; 
&world~model_he1i~body_hstab~alpha; 
iworld~mode l_he 1 i~body_hs 18Lb_beta ; 

&wor1d~mode1_he1i~body_hstab_qdyn; 

&wo r1d_mode1_he1i_body_hs t ab~ cd; 
iwo r1d_mode1~he1i_body~hstab_c1; 

&world_model~he1i_body~hstab_cm; 
fcworld_model_heli_body_hstab_pfo; 
&world_model_heli_body_hBtab_pfo(1:3); 
Siworld_model_heli_body_hstab_pf o (4:6); 
&world_model_heli_body_hstab_faero; 
&world_model_heli_body_hstab_faero(1:3); 
&world_model_heli_body_hstab_faero(4:6); 
&world_model_heli_body_hstab_pfrt; 


// Vertical tail aerodynamic data 


vtcv * 
vtiv s 
vttv ■ 
vtalfa = 
vtbeta = 
vtqdyn * 
cdvt = 
clvt = 
cmvt = 
vtlf 

vtfor >= 
vtmom = 
vtif * 
vtifor = 
vtimom » 
tvt 


&world_model_heli_body_vstab_pfrl(4:6); 
&world_data_vtiv; 

&world_model_he1i_body_vstab_ve1; 
&world_model_heli_body_vstab_alpha; 

&wor1d_model_he1i_body_vstab_beta; 

&wor1d_model_he1i_body_vstab_qdyn; 

&world_model_he1i_body_vstab_cd; 
&world_model_heli_body_vstab_cl; 
Stworld_model_heli_body_vstab_cm; 
&world_model_he1i_body_vstabjpfo; 
&world_model_heli_body_vstab_pfo(1:3); 
&world_model_heli_body_vstab_pfo(4:6); 
&wor1d_model_he1i_body_vstab_faero; 
&world_model_heli_body_vstab_faero(1:3); 
&world_model_heli_body_vstab_faero(4:6); 
&world_model~heli_body_vstcd)_p£rt; 


// Tail rotor data 


trthr = 
trtorq = 
trcv = 
triv = 
trtv = 
trlam = 
tromega= 
trrad = 
trlf 

trfor * 
trmom = 
trif 

trifor = 
trimooi = 
trtr = 


&world_model_heli_body_trotor_t; 
&world_model_heli_body_trotor_q; 
&world_model_heli_body_trotorjpfrl(4:6); 
&world_data_triv; 

&wor1d_mode1_he1i_body_t rotor_ve1; 

&wor1d_model_he1i_body_trotor_xl am ; 
&world_model_he1i_body_trotor_omeg; 
Stworld_model_heli_body_trotor_r ; 
Stworld_model_heli_body_trotor_pf o ; 
&world_model_heli_body_trotor_pfo(1:3); 
&world_model_heli_body_trotor_pfo(4:6); 
&wor1d_model_he1i_body_trotor_faero; 
&world_model_heli_body_trotor_faero(1:3); 
&world_model_heli_body_trotor_faero(4:6); 
&world_mode1_he1i_body_trotor_pfrt; 
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It Rotor data which includes inflow, tip path plane and hub 
for i»l:world_data_nblades, 

flap$i * tworld_model_heli_rotor_blade$i_flap_x; 
flapd$i m &world_model~heli_rotor_blade$i_flap_xd; 
flapdd$i * &world_model_heli_rotor_blade$i_flap_xdd; 
lag$i « &world_model_heli_rotor_blade$i_lead_x; 
lagd$i « 4world”model~heli_rotor_blade$i_lead~xd; 
lagdd$i « &worldjmodel_heli~rotor_blade$i lead_xdd; 
feat$i * tworld_model_helT_rotor_blade$T_feat_u(l); 
featd$i * 4world_model~heli_rotor_blade$i~feat~u(2); 
featdd$i = &world_model_heli_rotor_blade$i_feat~u(3); 

end 
i * l ; 

for j =o 1 :world_data_nseg, 

alfmr$i$j * &worTd_model_heli_rotor_blade$i_aero$j_alpha; 

mnmr$i$j = &world_model_heli_rotor_blade$i~aero$j_mach; 

clmr$i$j = &world_model~heli_rotor~blade$i~aero$j 3 cl; 

cdmr$i$j = &world_model_he 1 i~rotor_blade$i_aero$j _cd; 

cmmr$i$j = &wor 1 djmodel_he 1 i_rotor~blade$i_aero$ j _cm; 

qdyn$i$j * iwor 1 d_model_he 1 i_rotor_blade$i_aero$j_qdyn 

tv$i$j * &world_model_heli_rotor_blade$i_aero$j~vel; 

iv$i$j * &world_model~he 1 i~rotor_blade$i_aero$j_dwash; 

cv$i$j = &world_model_heli_rotor_blade$i_aero$j_pfrl(4:6); 

cfor$i$j = &world_model_heli_rotor_blade$i_aero$j jpfo(l:3); 

cmom$i$j s &world~model_heli_rotor_blade$i~aero$j_pfo(4:6); 

ifor$i$j * &world_model_heli~rotor_blade$i_aero$j_faero( 1 :3) 

imom$i$j = &world_model_heli_rotor~blade$i_aero$j_faero(4:6) 

end 

// Control system outputs 
xatrm= &world_data_xatrm; 
xbtrm= &world_data_xbtrm; 
xctrm= &world_data_xctrm; 
xptzm= &world_data_xptrm,- 

bis = &world_model_cont_lng_bls_y; 
als a &world_model_cont_lat_als_y; 
mth * &world_model_cont_coll_mtheta_y; 
tth = &wor1d_model_cont_dir_the 11r_y; 
iht = &world_data_htincang; 
xa = &world_data_xa; 

xb = £world_data_xb; 

xc = &world_data_xc; 

xp = &world_data_xp; 

lower = zeros( 6 , 1 ); 
upper = zeros( 6 , 1 ); 
lower( 1 ) =- 10 *d 2 r; 
lower(2) = -5*d2r; 
lower(3) = world_data_xcll; 
lower(4) = world_data_xall; 
lower(5) = world_data_xbl 1 ; 
lower(6) = world_data_xpll; 

upper( 1 ) = 10 *d 2 r; 
upper(2) = 5*d2r; 

upper(3) s world_data_xcul; 
upper(4) 3 world_data_xaul; 
upper(5) 3 world_data_xbul; 
upper(6) 3 world~data~xpul; 

popg 

// End of configure.exe 













/////////////////////////////////////////////////// 

// fils: assemble.exe 
// Assemble the helicopter model: 

// date: 21 July 1993 

iiimiiiiiiiiiiiiiituuiiuuiiiiiiiimuiuiiii 

topsolve_cf » 1; 

rhsolve_cf * 0; 
rotorsolve_cf ■ 0; 
helisolve_cf * 0; 
drivesolve_cf ■ 0; 
contsolve_cf « 0; 

drivesolve_reaBsemble « 0; 
drivesolveT:assemble; 
drivesolve::reset; 

contsolve_reassemble = 0; 
contsolve::assemble; 
contsolve::reset; 

// Run to steady state: 

for i*l:36; contsolve::step, contsolve::iter; end 
contsolve::saveic; 

exec("mbc_setup.exc",1); 

rotorsolve_reassemble = 0; 
rotorsolve::reset; 
rotorsolve::assemble; 

rhsolve_reassemble = 0 ; 

rhsolveT:reset; 
rhsolve::assemble; 

for i»l:10, rhsolve::genq; rhsolve::solve; end 
rhsolve::assemble; 

rhsolve_niter«0; rhsolve::iter, rhsolve_niter 
// End of assemble.exe 

/////////////////////////////////////////////////// 

// file: mbc_setup.exe 

// This file sets up the MBC transformation for the rigid rotor 
// date: 22 July 1993 

/////////////////////////////////////////////////// 

mbeset(model_heli_rotor_mbc, 1); 

// End of mbc_setup.exe 
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APFKMDXX C 


iiniiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiuiuiiiiuiiiiiu 

// file: trimsweep.def // 

// Definition file to conduct trim sweep for psh with rigid blades // 

IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIUIIIIIIIIIIIIIIII 

II execute the psh.def file to assemble the model for running 
exec{"psh.def *,1) 

// load additional functions used 

exec("limitchange.fun" , 1); // function to limit control change size 
exec("plotutils.exc",1) // plotting utilities 

exec("Trim_Sweep.exe") // Conduct trim sweep for the range of airspeeds 
desired 

// display trim control plots if desired 

showplt = Enter("Enter a 1 to plot control pen vs airspeed or 0 to 
exit:"); 
if showplt == 1 
plot("clear") 

plot ("title a Control Position vs Airspeed,xlabel * Veq(kts)") 
plot("ylabel ^Collective (inches)") 
plot(cpg_stc(l,:)',cpg_stc(4,:) ') 

pli * Enter("Enter 1 to print plot 0 to continue:*); 

if pli a= l; plot("print");elseif pli a= 0;disp("plot not printed");end 
plotCylabel aLongitudinal Cyclic (inches) ") 
plot(cpg_stc(l,:)',cpg_stc(6,:)') 

pli a Enter("Enter 1 to print plot 0 to continue:"); 

if pli a= i; plot("print");elseif pli aa 0;disp("plot not printed");end 
plotCylabel =Lateral Cyclic (inches)") 
plot(cpg_stc(l,:)',cpg_stc(5,:)') 

pli a Enter("Enter 1 to print plot 0 to continue:"); 

if pli aa l; plot("print");elseif pli a= O;disp("plot not printed");end 
plotCylabel aPedal (inches)") 
plot(cpg_stc(1,:)',cpg_stc(7,:)') 

pli a Enter("Enter 1 to print plot 0 to continue:"); 

if pli a= l; plot ("print") ,-elseif pli •* O;disp("plot not printed") ;end 
elseif showplt a= 0 

dispCTrim Control Matrix is saved in file called Trim_Controls.rbe") 

dispCTrim Matrix is saved in file called Trim_Matrix.rbe") 

end 

dispCend of trimsweep.def file") 

/////// end of trimsweep.def//////// 
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///////////////////////////////////////////////////////////////// 

// file: limitchange.fun 

// Thi* file limits the change of control during an update 
// date: 27 Auguat 1993 

////////////////////////////////////////////////////////////// 

function y * limitchange(x(n),pet,range(n)) 

// 

// 

// Y . LIMXTCHANGE(X,PCT,RANGE) 

// 

// Inputs: 

//. 

// X - Input vector 

// PCT - Allowable percentage change (0-1) 

// RANGE - Total range of X (must be same size as X) 

// 

// Outputs: 

//. 

// 

// Y - X limited to the allowable percentage of the range 

// 

// This function limits the amount the signal can change to a 
// specified percentage of the total range 

// 

function sign; 

y * x; 
for i«l:n 

if abs(x(i)) > pct*range(i) 

y(i) - sign(pct*range(i),x(i)); 
disp("Limiting Change on "$i) 

end 

end 

end 

//////////end of limitchange.fun ///////////// 
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IllllillllllllllllllllllUIIIIIIIIIIIIIIIIIIIIIIIIIIIIII 

// file: Trim_Sweep.exe 
// date: 31 August 1993 

// This file performs a trim sweep over a range of 
// airspeeds in increment specified by user 

IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIUIIIIIIIIIIIIII 

pushg(world_cpg) 

nprev * world data_naz; //# of azimuth steps/revolution 
nrevss » 4; //I of revolutions to reach steady state 

nrevavg » 1; //# of revolutions for obtaining average 

converg ■ 0.0001; 

bias ■ zeros(6,1); 

Trimg ■ ones(bias) "Trim Gain (nd)"; 

//exec("Compute_Trim_Matrix.exe",1) 

//disp("Trim Matrix") 

//disp( (trm) ) 
stc * [] ; 
output ( []) 
output(bacc) 
popg 

exec("Set_Airspeed.exe",1) // Set airpseed for initial setup 
exec("Steady_State.exc",1) // Run model to steady state 
exec("Steady_State.exc",1) 
exec("Steady_State.exc", 1) 
exec("Steady_State.exc",1) 
exec("Steady_State.exc",1) 
exec("Steady_State.exc",1) 
exec("Steady_State.exe",1) 
exec("Steady_State.exc",1) 

pushg(world_cpg) 
vstart • veq; 

vinc m Enter("Enter increment velocity :"); 
vend * Enter("Enter ending velocity :"); 

tc = [] ; 
trmd= [] ; 
output ( []) 
output(bata,bara) 


for ijk= vstart:vine:vend 
veq = ijk; 

trimg = ones(trimg); 

exec("Compute_Trim_Matrix.exc",1) 

exec("Save_Trim_Matrix.exc",1) 

disp((trm) ) 

exec ("Trim.exc" , 1) 

tc=(veq ean(2) ean(l) xctrm xatrm xbtrm xptrm mth als bis tth gamvr 
gamhr]'; 

exec("Save_Controls.exc",1) 
end 

save("Trim_Controls.rbe",stc) 
save("Trim_Matrix.rbe",trmd) 

popg 

// End of Trim_Sweep.exe //// 
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///////////////////////////////////////////////////// 

// file: Compute_Trim_Matrix.exe 

// This script computes a diagonal matrix which is 
// the partial of one acceleration which is coupled to 
// one control. 

// date: 27 July 1993 

iiiiiiiiiiiiiiniiiiuiiiiiiiiiinniniiiiiiiiiiiiui 

// This script computes a diagonal matrix which is 
// the partial of one acceleration which is coupled to 
// one control. 

// 

II th - ud (pitch attitude couples with longitudinal accel) 


// ph - vd (roll attitude couples with lateral accel) 
II xc - wd (collective couples with vertical accel) 
// xa - pd (lateral cyclic couples with roll accel) 
// xb - qd (long. cyclic couples with pitch accel) 
// xp ~ rd (pedal couples with yaw accel) 


// 

d(ud) /d(th) 

0 

0 

0 

0 

0 

// 

0 d (vd) 

/d(ph) 

0 

0 

0 

0 

// 

0 

0 

d(wd)/d(xc) 

0 

0 

0 

// 

0 

0 

0 

d(pd) /d(xa) 

0 

0 

// 

0 

0 

0 

0 

d (qd) /d(xb) 

0 

// 

II 

0 

0 

0 

0 

0 

d(rd) /d(xp) 

pushg(world_cpg) 







// Steady_State.exe scripts computes the average values 
// of the output if outputs have been selected. The results 
// are returned in aO. 


output ([] ) 
output(bacc) 

exec("Steady_State.exc",l) 
exec("Steady_State.exc",1) 
exec("Steady_State.exe",1) 
aOO = aO; 

xO = ean(2); 

trm = zeros(6,6) ; 


i = l; 

ean(2) = ean(2) + 0.01; 
exec("Steady_State.exc", 1) 
trm(i,i)= (aO(i) - aOO(i))/0.01; 
ean(2) * ean(2) - 0.01; 

i » 2; 

ean(l) * ean(l) + 0.01; 
exec("Steady_State.exc",1) 
trm(i,i)« (a0(i) - aOO(i))/0.01; 
ean(l) - ean(l) - 0.01; 

i * 3 ; 

xctrm * xctrm + 0.01; 
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exec("Steady_State.exc",1) 
trm(i,i)» (aO(i) - aOO(i))/0.01; 
xctrm ■ xctrm - 0.01; 


i - 4; 

xatrm « xatrm + 0.01; 
exec<*Steady_State.exc",1) 
trm(i,i)» (aO(i) - aOO(i))/0.01; 
xatrm » xatrm - 0.01; 

i - 5; 

xbtrm > xbtrm + 0.01; 
exec("Steady_State.exc", 1) 
trm(i,i)a (a0(i) - aOO(i))/0.01; 
xbtrm a xbtrm - 0.01; 

i - 6; 

xptrm a xptrm + 0.01; 
exec("Steady_State.exc",1) 
trm(i,i)a (aO(i) - aOO(i))/0.01; 
xptrm = xptrm - 0.01; 

trm a inv(trm); 
trm a - trm; 

save("Trim_Matrix.rbe",trm) 
popg 

// End of Compute_Trim_Matrix.exe 

////////////////////////////////////////////////////// 

// file: Trim.exc 

// This script applies an extremely simplistic trim algorithm 
// to drive the accelerations to zero. 

// date: 3 Setember 1993 

////////////////////////////////////////////////////// 

// This script applies an extremely simplistic trim algorithm to 
// drive the accelerations to zero, 
pushg(world_cpg) 

// Set body velocities from prescribed flight condition 
exec("Conpute_Body_Axis_Velocity.exe",1) 
world_rhsolve7:assemble 
world_topsolve::iter; 

// Run to steady state 
output ([]) 

output(world_cpg_bacc) 
dispCRun to steady state"); 
exec("Steady_State.exc",1); 

//aaaaaaaaaaaaaaaaaaaaaTRXM LOOP j****************************** 

// 

itercntaO; 

trimfailaO; 

disp("Perform gradient trim"); 
while norm(aO-bias) > converg && iterent < 20 
iterent * iterent + 1; 
exec("Update_CW.exc",1); 

disp("Iteration count is "$itercnt$" at "$veq$" knots"); 
disp("Average accelerations are"); 
disp((a0')) 
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if (itercnt«»20) 
trimfail-1; 

end 

end 

// 

//if crim fails, reduce gradient step by 50% and try once more 

// 

if (triufail»l) 

world cpg trimg*.5*world cpg trimg; 

// ~ 

//*********** **********<TRIM LOOP 2******************************** 

// 

itercnt-0; 

trimfail»0; 

disp{"Perform gradient trim"); 
while norm(aO-bias) > converg && itercnt < 20 
itercnt * itercnt + 1; 
exec("Update_CW.exc",1); 

disp("Iteration~count is "$itercnt$" at "$veq$" knots"); 

disp("Average accelerations are"); 
disp((aO')) 

if (itercnt**20) 
trimfailsl; 

end 

end 

//********** . **************...******.* 

end 

// 

//if trim fails, reduce gradient step by 50% and try once more 

// 

if (trimfail==l> 

world cpg trimg=.5*world cpg trimg; 

// ~ - - 

//*********************TR IM LOOP 3******************************** 

// 

itercnt=0 ,- 
trimfail=0 ; 

disp("Perform gradient trim"); 
while norm(aO-bias) > converg && itercnt < 20 
itercnt = itercnt + 1; 
exec("Opdate_CW.exc",1); 

disp("Iteration^count is "$itercnt$" at "$veq$" knots"); 

disp("Average accelerations are"); 
disp((aO')) 

if (itercnt=*20) 
trimfail=l; 

end 

end 

//************************.****************************** 

end 


p°pg 

// End of Trim.exc 

/////////////////////////////////////////////////////////////// 

// file: Update_CW.exc 

// This file update controls and modified trim weights 
// date: 27 August 1993 

//////////////////////////////////////////////////////////////// 

aOO * aO; // Save the last average accels 
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exec("Update.exc",1) // Update controls 

for i»l:6 // Zf the sign of the accels changed 

sal = aO(i)/abs(aO(i));// Then decrease that trim gain by 
saO aOO (i) /abs(aOO (i)) ;// by 10 percent. Ijower limit is 0.5 
if sal <> saO 

trimg(i) a max(0.5,0.9*trimg(i)) ; 
end 

end 

// End of Update_OJ.exc 

///////////////////////////////////////////////////////////////// 

// file: Ujpdate.exe 

// This file update the controls 

// date: 27 July 1993 

iiiiiiiiiiiiiiiiiiiiiiiiimiiiiiiiiiiiiiiiiiiiiiniiiiiiiini 

pushg(world_cpg) 

deltac » trm*(trimg.*a0); 

deltac a limitchange(deltac,0.02,abs(upper-lower)); 

ean(2) = min(upper(1),max(lower(1),ean(2) + deltac(1))); 
ean(l) a min(upper(2),max(lower(2),ean(1) + deltac(2))); 
xctrm = min(upper(3),max(lower(3),xctrm + deltac(3))); 
xatrm a min(upper(4),max(lower(4),xatrm + deltac(4))); 
xbtrm a miniupper(5),max(lower(5),xbtrm + deltac(5))); 
xptrm a min(upper(6),max(lower(6),xptrm + deltac(6))),- 

dispCNew control setting (theta phi xc xa xb xp)") 
disp i [ [ean (2) eaui(l)]*r2d xctrm xatrm xbtrm xptrm]) 
execi"Compute_Body_Axis_Velocity.exc",1) 

p°pg 

exec("Steady_State.exe",1) 

// end of Update.exe 

/////////////////////////////////////////////////////////// 

// file: Compute_Body_Axis_Velocity 

// This file uses the equivalent airspeed and the flight path 
// angle sines and cosines to set the aircraft body axis velocity 
// date: 4 Nov 1992 

/////////////////////////////////////////////////////////// 

pushg(world_cpg) 

world_rhsolve::genq; 

vt a veq*sqrt(rhoO/rho)*k2f; 

cgv a cos(gamvr); sgv a sin(gamvr); 

cgh a cos igamhr) ,• sgh a sin(gamhr) ; 

vi a vt*(cgv*cgh cgv*sgh -sgv]'; 

tib(1:3,1:3) a tran; 

bvel a tib*vi; 

world_rhsolve::genq 

popg 

// End of Cc«npute_Body_AxiB_Velocity 
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fiitiiiiniiiiiiiiiiiinntiiiiiuiiiiiiiiiiiiiiiiiiiiiiuit 

// file: Save_Trim_Matrix.exe 
// date: 2 September 1933 

// Saves the trim matrix for each airspeed as a vector 

iiiniiiiiiiiiiimiiiiiiiiiuiiiiuiiiiiiumiiiiiiiu/iiii 

pushg(world_cpg) 

trmda[trmd (veq diag(trm) 

describe("Saved Trim Matrix: veq d(ud)/d(th) d(vd)/d(ph) d(wd)/d(xc) 
d(pd)/d(xa) d(qd)/d(xb) d(rd)/d(xp)■); 

popg 

// Bnd of Save_Trim_Matrix.exe 

//////////////////////////////////////////////////////////// 

// file: Save_Controls.exe 
// date: 27 July 1993 

// Include the trim controls which you wish to save 
// in this list 

lllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 

pushg(world_cpg) 

stc*[stc [veq ean(2) ean(l) xctrm xatrm xbtrm xptrm mth als bis tth 
gamvr gamhr]']; 

describe("Saved Trim Controls: veq ean(2) ean(l) xctrm xatrm xbtrm xptrm 
mth als bis tth gamvr gamhr"); 

popg 

// Bnd of Save_Controls.exe 

llllllllllllllllllllllllllllllllllllllllllllllllllllll 

II file: Set_Airspeed.exe 

// User specifies airspeed in knots 

// vertical and horizontal flight path angle in degrees 
// date: 3 August 1993 

iiiiiiiiniiiiiiiiiiiiiiiiiuiiiiiiiiiiiiiiiiiniiniii 

pushg(world_cpg) 

veq = Enter("Enter equivalent airspeed (knots) : "); 

gamvr = Enter("Vertical flight path angle(deg) + down: ")*d2r; 
gamhr = Enter("Horizontal flight path angle(deg) + r : ")*d2r; 
exec("Compute_Body_Axis Velocity.exe",1) 
popg 

// End of Set_Airspeed.exe 

IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIUIIII 
II file: Steady_State.exe 

// This script runs the model to steady state 
// date: 27 July 1993 

////////////////////////////////////////////////////////// 

pushg(world_cpg) 
y = nrun(nprev*nrevss) ; 
if y <> [J 

plot("xlabel=Time steps for nrevss,ylabelsaverage body accels") 

// plot(nprev*nrevss*world_data_dt,y) 
plot(y) 

sk = ((nrevss-nrevavg)*nprev+l):(nrevss*nprev); 
aO = sum(y(sk,:))'/(nprev*nrevavg); 
end 
popg 

// End of Steady_State.exe 
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APPENDIX D 


IlllllllllltlllllllllllllllllllllllllllUIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII 

// file: analysis.def // 

// Definition file to analyze the psh with rigid blades // 

// using the linearization routine and sample tests // 

IIIIIIIIIIIIIIIIIIIIIIUIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII 

II execute the psh.def file to assemble the model for running 
exec("psh.def",1) 

// add path to directory where test files are located 
path("tests/") 

// load additional functions 

exec("qsreduce.fun",1); // linear reduction function 
exec("plotutils.exc",1) // plotting utilities 

exec("logspace.fun.fun",1); // function for creating frequency range 

exec("freq.fun",1); // function for the freq sweep test 

exec("Restore_Trim.exe") // Restores the model to trim at selected 
airspeed 

exec("Steady_State.exc",1) 
exec("Steady_State.exc",1) 
exec("Steady_State.exc",1) 
exec("Steady_State.exc",1) 

// Run the model and determine the non linear characteristic matrix 

frh=nlmodel(rhsolve); // open loop system matrix based on rhsolve 

solution 

// To compare the nonlinear results with a reduced 
// order linear model, you need to create a reduced 
// order linear model. 

cpg_havembc=l; //flag to indicate use of multi-blade coord system 
exec("6dof_linearize.exc",1); // returns linear matrix for 6dof model 
// Where 

// full f matrix in f 
// full g matrix in g 
// reduced f matrix in fq 
// reduced g matrix in gq 

// change ndof to 10, 11, 15 in 6dof _linearize.exe to get other 
// reduced models 

// Setup matrices to compare linear and nonlinear model response 
// Select the states (u,v,w,p,q,r,phi,theta) as outputs 
// by setting the hq matrix to identity and the dq matrix 
// to zero; 

hq = eye (fq) ; 
dq = zeros(gq); 

// Construct a linear system matrix 
sq = [fq gq 
hq dq]; 
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// Discretize the linear system matrix (convolution integral) 
sqd * disc(sq,8,topsolve_dt); 


// Make a plot of the eigenvalues for the linear system matrix sq 

plot("clear") 
eignuxUeig(sq(l:8,l:8)) ; 

plot("titleaBigenvalues Linear Model,xlab =Real Axis,ylab=Imag Axis") 

plot("mark=l,line=0"); 

plot(real(eigmod),imag(eigmod)); 

pli * Enter("Enter 1 to print plot 0 to continue:"); 

if pli == l; plot("print");elseif pli == 0;disp("plot not printed");end 
plot("clear") 

disp(”The reduced linear characteristic matrix") ;disp( (fq)) ,-pause 
disp("The reduced linear control matrix") ,-disp ((gq)) ,-pause 

{V,D}=eig(fq) ,- //eigenvectors and eigenvalues of linear matrix 
ceqn=poly(fq); //characteristic equation of linear matrix 

disp("The characteristic equation coefficients of the linear matrix fq") 

disp((ceqn)) 

pause 

disp("The eigenvectors of the linear characteristic matrix") 

disp ((V)) 

pause 

disp("The eigenvalues of the linear characteristic matrix") 

disp((diag(D))) 

pause 

// conduct longitudinal axis tests 

exec("IiOng_Impulse.exc",1) //use 60 sec, 0.025 sec delay, 1 inch 

// Propagate the linear model, the input vector must 
// have the same dimensionality as the gq matrix 
// which is [xb xa xp xc]' 

yllng = filp(sqd,[ulng 0*ulng 0*ulng 0*ulng]); 

// Get the bias of the nonlinear run 
xO = ylng(1,1:8); 

// Add the nonlinear bias to the linear model results 
// and plot the results 
for i=l:8 

yllng(:,i) = yllng(:,i) + x0(i); 
end 

// Plot the response of the linear and nonlinear models 
plot("clear") 

plot("title=Phugoid Response") 

plot("upper,xlab=time (sec),ylab=theta (deg)") 

plot(tv,ylng(:,8).*world_data_r2d,tv,yllng(:,8).*world_data_r2d) 
pli = Enter("Enter 1 to print~plot 0 to continue:"); 

if pli == 1; plot ("print") ;elseif pli == 0;disp("plot not printed") ,-end 
plot("lower,title= ,ylab=velocity (fps)") 
plot(tv,ylng(:,1),tv,yllng(:,1)) 

pli = Enter("Enter 1 to print plot 0 to continue:"); 

if pli == 1; plot("print");elseif pli == 0;disp("plot not printed");end 
// end of longitudinal inpulse test 
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// conduct lateral axis tests 

exec("Lat_Step.exc",l) //use 15 sec,.025 delay,1.0 inch,2.5 sec, .025 
rise/fall 

// Propagate the linear model and 
// Get the bias of the nonlinear run 

// Add the nonlinear bias to the linear model results 
// and plot the results 

yllat « filp(sqd,[0*ulat ulat 0*ulat 0*ulat]); 
for isl:8 

yllat(:,i) - yllat(:,i) + x0(i); 
end 

plot("clear") 

plot("title=Roll Response") 

plot("upper,xlabstime (sec),ylab=Phi (deg)") 

plot(tv,ylat(:,7).*world_data_r2d,tv,yllat(:,7).*world_data_r2d) 
pli ■ Enter("Enter 1 to print plot 0 to continue:"); 

if pli l; plot("print");elseif pli » 0;disp("plot not printed");end 
plot("lower,title= ,ylab»Roll Rate (deg/sec)") 

plot(tv,ylat(:,4).*world_data_r2d,tv,yllat(:,4).*world_data_r2d) 
pli = Enter("Enter 1 to print plot 0 to continue:"); 

if pli == 1; plot("print");elseif pli ** 0;disp("plot not printed");end 

exec("Lat_Impulse.exe",1) //use 30 sec, 0.025 sec delay, 1 inch right 
stick 

// Propagate the linear model and 
// Get the bias of the nonlinear run 

// Add the nonlinear bias to the linear model results 
// and plot the results 

yllat2 = filp(sqd,[0*ulat2 ulat2 0*ulat2 0*ulat2]); 
xO = ylat2(1,1:8) ; 
for i=l:8 

yllat2(:,i) = yllat2(:,i) + x0(i); 
end 

plot("clear") 

plot("title=Spiral Response") 

plot("upper,xlab=time (sec),ylab=Phi (deg)") 

plot(tv,ylat2(:,7).*world_data_r2d,tv,yllat2(:,7).*world_data_r2d) 
pli = Enter("Enter 1 to print plot 0 to continue:"); 

if pli == 1; plot("print");elseif pli == 0;disp("plot not printed");end 
plot("lower,title= ,ylab=Roll Rate (deg/sec)") 

plot(tv,ylat2(:,4).*world_data_r2d,tv,yllat2(:,4).*world_data_r2d) 
pli a Enter("Enter 1 to print plot 0 to continue:"); 

if pli == 1; plot("print");elseif pli == 0;disp("plot not printed");end 
//conduct vertical axis tests 

exec("Ped_Doublet.exe",1) // use 20 sec,0.025,1 inch,1.5,.025,.025,.025 

// Propagate the linear model and 
// Get the bias of the nonlinear run 

// Add the nonlinear bias to the linear model results 
// and plot the results 

ylped = filp(sqd,[0*uped 0*uped uped 0*uped]); 
xO = yped(1,1:8); 










for i*l:8 

ylped(:,i) » ylped(:,i) + x0(i); 

end 

plot("clear") 

plot("title*Dutch-roll Response") 

plot("upper,xlabxtime (sec),ylab=Phi (deg)") 

plot(tv,yped(:,7).*world_data_r2d,tv,ylped(:,7).*world_data_r2d) 
pli * Enter("Enter 1 to print plot 0 to continue:"); 

if pli ■■ i; plot("print");elseif pli » O;disp("plot not printed”);end 
plot("lower,left,titles ,xlab=time (sec),ylab*Roll Rate (cteg/sec)") 
plot(tv,yped(:,4).*world_data_r2d,tv,ylped(:,4).*world_data_r2d) 
pli = Enter("Enter 1 to print~plot 0 to continue:"); 

if pli == 1; plot ("print") ,-elseif pli == O;disp("plot not printed") ;end 
plot("lower,right,titles ,xlab=time (sec),ylab=Yaw Rate (deg/sec)") 
plot(tv,yped(:,6).*world_data_r2d,tv,ylped(:,6).*world_data_r2d) 
pli * Enter("Enter 1 to print~plot 0 to continue:"); 

if pli *x 1; plot ("print") ,-elseif pli ■■ 0,-disp ("plot not printed") ;end 
// conduct freq response test 

// simplify system matrix to have only the input column desired for the 
test 

// you want and the only the ouputs for that test 

// example frequency response of q/long stick and theta/long stick 
[f,g,h,d]=split(sq,8); 

g=g(:,1); // Pick up only the column associated with long stick 

// Output index for pitch rate is 5 (u v w p q r phi theta) and 8 for 
// theta, so 
sk = [8] 
h = h(sk,:) ; 

[no,ns]=size(h); 

[ns,ni]ssize(g); 
d=zeros(no,ni); 
s=[f g; 
h d] ; 

w=logspace(0.1,100,200) ; 

[amp,phase]=freq(s,8,w); 

// The scope plot package does not handle log plots very well, so to get 
// equivalent system, convert the data to a log base 10 scale. The 
// x axis now represent powers of 10. 

sc=l/log(10); //conversion from In to log base 10 
scg=20/log(10); //conversion from In to db 

// plot pitch attitude response to longitudinal freq sweep 
plot("clear") 

plot("upper,title=frequency response theta vs delta long,xlab= 

,ylab=Gain(db)") 

plot(sc*log(w),scg*log(amp)) 

pli * Enter("Enter 1 to print plot 0 to continue:"); 

if pli == 1; plot ("print") ,-elseif pli == 0;disp("plot not printed") ,-end 
plot("lower,title= ,xlab=Frequency ((rad/sec) in power of 10),ylabsPhase 
(deg)") 

plot(sc*log(w),phase) 

pli = Enter("Enter 1 to print plot 0 to continue:"); 

if pli == 1; plot ("print") ,-elseif pli == 0,-disp ("plot not printed") ,-end 
disp("end of analysis.def file") 

//end of analysis.def 
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iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiniiniiiiiiiiiiiiiiiiiiiiiiiiniiii 

// file: Restore_Trim.exe 
// date: 27 July 1993 

// This script restores the trim matrix and trim control matrix 
// and runs the model to steady state for the selected a/s 

llinillllllllllllllllllllllllllllllllllllllllllllllllllllllll 

pushg(cpg) 

load("Trim_Controls.rbe") 
load("Trim_Matrix.rbe") 
lr,c]-size(cpg_stc); 

test=diag(ones(c));for i = 1:c;test(i)*test(i)*i;;end 
disp(([test cpg_stc(l,:)'])) 

Ti « Enter("Enter the row number of the airspeed to analyze:"); 

world_cpg_veq * cpg_stc(1,ti); 

world_cpg_ean(2) * cpg_stc(2,ti); 

world_cpg~ean(l) = cpg_stc(3, t i); 

world_cpg_xctrm = cpg_stc(4,ti); 

world_r <g”xatrm = cpg_stc(5,ti); 

world_ g_xbtrm = cpg_stc(6,ti); 

world_qpg~xptrm = cpg_stc(7,ti); 

world_cpg~geunvr = cpg_stc (12, ti) ; 

world_cpg_gamhr * cpg_stc(13,ti); 

trm*diag(cpg_trmd(2:7,ti)'); 
disp("Trim Matrix="); disp((trm)) 
tc=cpg_stc(:,ti)'; 

disp ("Trim Control Matrix=") ,-disp ((tc)) 

nprev = world_data_naz; 

nrevss = 4; 

nrevavg = 1; 

converg = 0.0001; 

bias = zeros(6,1); 

Trimg = ones(bias) "Trim Gain (nd)"; 

3tc = [] ; 

output ( (] ) 
output(bacc) 

veq = cpg_tc(l,l); 

disp("Equivalent Arispeed =");disp((veq)) 

gamvr = Enter("Enter Vertical flight path angle(deg) + down: ")*d2r; 
gamhr = Enter("Enter Horizontal flight path angle(deg) + r : ")*d2r; 
exec("Compute_Body_Axis_Velocity.exe",1) 

p°pg 

output ( [] ) 

output(world_cpg_bacc) 
exec("Steady_State.exe",l) 
exec("Steady_State.exc",1) 
exec("Steady_State.exc" ,1) 
exec("Steady_State.exc",1) 
exec("Steady_State.exc", 1) 
exec("Steady_State.exc",1) 
exec("Steady_State.exe", l) 
exe c("Steady_State.exc",l) 
disp("Average accelerations") 
disp((world_cpg_aO)) 

// End of Restore Trim.exc 
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///////////////////////////////////////////////////////////////// 

// file: Long_Impulse.exe 

// This file Set • up a longitudinal inpulse 

// date: 3 August 1993 

/////////////////////////////////////////////////////////// 

// This file sets up a Longitudinal Impulse 

pushg(world_cpg) 

exec("Select_Outputs.exc", 1); 
input {[]) 
input(xb); 
popg 

tend >■ Enter ("Enter run duration (sec) : ”); 

tdly « Enter("Enter delay time of impulse (sec): "); 
ampl a Enter("Enter inpulse anplitude pos+ (inches) : "); 

ulng a ipulse(world_topsolve_dt,tend,tdly,anpl) "Impulse input signal 
(in)"; 

[r,c]asize(ulng); 

tv a world_topsolve_dt*(0:r-l) ' "Time Vector (sec)"; 
plot("clear") 

plot("xlabaTime (sec) f ylab=Long. Input Signal") 

plot(tv,ulng) 

world::saveic 

exec("Unfreeze.exc",1) 

ylng a nrun(ulng) "Nonlinear model response to longitudinal impulse"; 
exec("Freeze.exc",1) 
plot("clear") 

plot("xlabel = Time (sec),ylabel a Theta (deg) ") 
plot(tv,ylng(:,8)*world_data_r2d) 

//pli a Enter("Enter 1 to print plot 0 to continue:"); 

//if pli aa i; plot("print");elseif pli «« 0;disp("plot not 
printed") ,-end 

plot("ylabel * Velocity (fps) ") 
plot(tv,ylng{:,1)) 

//pli * Enter("Enter 1 to print plot 0 to continue:"); 

//if pli =a l; plot("print");elseif pli aa 0;disp("plot not 
printed") ,-end 
world::restoreic; 
pushg(world_cpg) 
output ([]) 
output(bacc) 
popg 

plot("clear") 

exec("Steady_State.exc",l) 
dispC’end of Long_Inrpusle.exe") 

// End of Long_Impulse.exe 
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/////////////////////////////////////////////////////// 

// file: Lat_Step.exe 

// This file sets up a Lateral Step 

// date: 3 August 1993 

iiiiiiiiiiiiiiiiiihiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii 

// This file sets up a Lateral Step 

pushg(world_cpg) 

exec("Select_Outputs.exc", 1); 
input ([] ) 
input(xa); 
popg 


tend * Enter("Enter run duration (sec) : "); 
tdly * Enter("Enter delay time of step (sec): ") ; 
ampl a Enter(”Enter step amplitude (inches) : "); 
delta* Enter("Enter step duration (sec) : ") ; 
trise* Enter("Enter rise time (sec) : ") ; 
tfall= Enter("Enter fall time (sec) : "); 


ulat * istep(world_topsolve dt,tend,tdly,delta,triae,tfall,ampl); 
describe(ulat,"Step input signal (in)"); 

[r, c] *size (ulat) ,- 

tv * world_topsolve_dt*(0:r-l)' "Time Vector (sec)"; 
plot("clear") 

plot("xlab=Time (sec),ylab=Input Signal") 
plot(tv,ulat) 
world::saveic 
exec("Unfreeze.exe" ,1) 

ylat = nrun(ulat) "Nonlinear model response to lateral step"; 
exec("Freeze.exc",1) 
plot("clear") 

//plot("title a NL Roll Response a "$veq$" kts") 
plot("xlabel = Time (sec),ylabel * Phi (deg) ") 
plot(tv,ylat(:,7)*world_data_r2d) 

//pli * Enter("Enter 1 to print plot 0 to continue:"); 

//if pli == 1; plot ("print") ,-elseif pli == 0;disp("plot not 
printed");end 

plot("ylabel * Roll Rate (deg/sec)") 
plot(tv,ylat(:,4)*world_data_r2d) 

//pli = Enter("Enter 1 to print plot 0 to continue:"); 

//if pli == 1; plot ("print") ,-elseif pli == 0;disp("plot not 
printed");end 
world::restoreic; 
pushg(world_cpg) 
output ( [] ) 
output(bacc) 
popg 

plot("clear") 

exec("Steady_Stete.exe", 1) 
disp("end of Lateral_Step.exe") 

// End of Lat_Step.exe 
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/////////////////////////////////////////////////////// 

// file: Lat_Xapulse.exe 

// This file sets up a Lateral Impulse 

// date: 3 August 1993 

/////////////////////////////////////////////////////// 

// This file sets up a Lateral Impulse 
pushg(world_cpg) 

exec(■ Select_Outputs.exc",1); 
input ([]) 
input (xa) ; 
popg 

tend « Enter("Enter run duration (sec) : "); 

tdly * Enter("Enter delay time of impulse (sec): "); 
ampl « Enter ("Enter inpulse amplitude right* (inches) : ") ; 

ulat2 ■ ipulse(world_topsolve_dt,tend,tdly,ampl) "impulse input signal 
(in)"; 

[r,c] =»size (ulat2) ; 

tv » world_topsolve_dt*(0:r-l) ' "Time Vector (sec)"; 
plot("clear") 

plot("xlab*Time (sec),ylab*Input Signal") 
plot(tv,ulat2) 

world: .-saveic 

exec("Onfreeze.exe", 1) 

ylat2 * nrun(ulat2) "Nonlinear model response to longitudinal impulse"; 
exec("Freeze.exc",1) 
plot("clear") 

//plot("title * NL Spiral Response ft "$veq$" Jets") 
plot("xlabel * Time (sec).ylabel • Phi (deg) ") 
plot(tv,ylat2(:,7)*world_data_r2d) 

//pli = Enter("Enter 1 to print plot 0 to continue:"); 

//if pli 1; plot("print");elseif pli 0;disp(”plot not 
printed") ,-end 

plot("ylabel * Roll Rate (deg/sec)") 
plot(tv,ylat2(:, 4)*world_data_r2d) 

//pli = Enter ("Enter 1 to print plot 0 to continue-."); 

//if pli =* 1; plot("print");elseif pli *■ O;disp("plot not 
printed");end 
world: .-restoreic; 
pushg(world_qpg) 
output ([]) 
output(bacc) 

pqpg 

plot ("clear") 

exec("Steady_State.exc",1) 
exec("Steady_State.exc", 1) 
disp("end of Lateral_Impulse.exe") 

// End of Lat_lmpulse.exe 
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iiiiiiiiniiiiiiiiiiiiiiiiiiiiiniiiiiiiiiiiiiiiiuiiii 

// file: Ped_Doublet.exe 

// This file sets up a Pedal Doublet 

// date: 3 August 1993 

IIIIIIIIIIIIIIIUIIIIIIIIIHIIIUIIIIIIIIIIIIUIIIIIUI 

// This file sets up a Pedal Doublet 

pushg(world_cpg) 

exec("Select_Outputs.exc", 1); 
input(U) 
input(xp); 
popg 


tend - Enter("Enter run duration (sec) : *) 
tdly > Enter("Enter delay time of doublet (sec): ■) 
ampl a Enter("Enter step amplitude (inches) : ") 
delta- Enter("Enter step duration (sec) : ") 
trisea Enter("Enter rise time (sec) : ") 
tfall* Enter("Enter fall time (sec) : ") 
tnxt a Enter("Enter time between steps (sec) : ") 


uped a 

idoublet(world_topsolve dt,tend,tdly,delta,trise,tfall,tnxt.ampl); 
describe(uped,"Doublet Input signal (in)”); 

[r,c]asize(uped); 

tv a world_topsolve_dt*(0:r-l)' "Time Vector (sec)"; 
plot("clear") 

plot("xlabaTime (sec),ylabainput Signal") 

plot(tv,uped) 

world::saveic 

exec("Unfreeze.exc",1) 

yped a nrun(uped) "Nonlinear model response to pedal doublet"; 
exec("Freeze.exc" , 1) 
plot("clear") 

plot("xlabel a Time (sec) ");plot("ylabel « Phi (deg) ") 
plot(tv,yped(:,7)*world_data_r2d) 

//pli a Enter("Enter 1 to print plot 0 to continue:*); 

//if pli aa l; plot("print");elseif pli aa 0;disp("plot not 
printed") ,-end 

plot("ylabel a Roll Rate (deg/sec)”) 
plot(tv,yped(:,4)*world_data r2d) 

//pli a Enter("Enter 1 to print plot 0 to continue:"); 

//if pli aa l; plot("print");elseif pli aa 0;diBp("plot not 
printed");end 

plot("ylabel a Yaw Rate (deg/sec)") 
plot(tv,yped(:,6)*world_data_r2d) 

//pli a Enter("Enter 1 to print plot 0 to continue:"); 

//if pli aa 1 ; plot{"print");elseif pli aa 0;disp("plot not 
printed");end 
world::re8toreic; 
pushg(world_cpg) 
output ([] ) 
output(bacc) 
popg 

plot("clear") 

exec("Steady_State.exc", 1) 
exec("Steady_State.exc",1) 
disp("end of~Pedal_Doublet.exe") 

// End of Ped Doublet.exe 
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/////////////////////////////////////////////////////// 

// file: Select Output*.exc 

// This file selects the outputs for the tests 
// date: 3 August 1993 

/////////////////////////////////////////////////////// 

output (U) 

output(bvel,brat,ean) 
for i»l:vorld_data_nblades 
output(flap$i) "" 
end 

for i»l:world_data_nblades 
output(lag$i) 

end 

describe("Outputs 1-3: Body axis translational velocities u,v,w 
(ft/sec)"); 

describe("Outputs 4-6: Body axis rotational velocities p,q,r 
(rad/sec)"); 

describe("Outputs 7-9: Euler angle roll, pitch yaw (rad)"); 
describe("Outputs 10-13: Flap response of blades 1,2,3,4 (rad)"); 
describe("Outputs 14-17: Lag response of blades 1,2,3,4 (rad)"); 

// End of Select_Outputs.exe 

//////////////////////////////////////////////////////////// 

// File: 6dof_linearize 

// This file generates a 6dof linear model for pah: 

// date: 27 July 1993 

llllllllllllllllllllllllllllllllinillllllllllllllllllllll 

[f,g]=LinearModel(topsolve); 

(ns,ns]«size(f); 

// Fetch state and input indices of relevant components: 
ix6*StateIndices(topsolve,world_model_heli_body_bodydof) "Body state 
indices"; 

ix6r»ix6([13:18 4 5]) "Reduced body state indices"; 
iup*[ InputIndices(topsolve,vorld_model_cont_lng_xb); 

Inputlndices(topsolve,world~model_cont”lat~xa); 

Inputlndices(topsolve,world_model_cont_dir_xp); 

Inputlndices(topsolve,wor1d~model_cont~coll_xc)] "Input indices" 
nup*prod(size(iup)); 
savg=[f g(:,iup)]; 

// average linear model if required 

nrun(3); 

[f,g]=LinearModel(topsolve); 
s*[f g(:,iup)] ; 
savg»savg+s; 
nrun(3); 

[f,g]=LinearModel(topsolve) ; 
s*[f g(:,iup)]; 
savg* (savg-fs) /3 ; 
f*savg(1:ns,1:ns); 
gpssavg(1:ns,ns+l:ns+nup); 

//Set up number of DOFs for reduced order model 
ndof*6; 

//Perform quasistatic reduction 
exec("qsreduce.exc", l); 

// end of 6dof_linearize.exe 
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Illlllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 

// File: qsreduce.exe 

// Author: Ron DuVal and Joe English 

// Description: Computes quasistatically reduced linear model of 

helicopter. 

// See-Also; linearize.exe 

// Notes: Most of the inputs to this script are produced by 

linearize.exe 

llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 


// 

II 

II 

II 

// 

II 

// 

II 

// 


INPUTS: 

NDOF 


6 

10 

11 

15 

0 


# degrees of freedom in reduced model 
Must be set to: 

DOF6 body states only 

DOF6 body states + flapping 

DOF6 body states + flapping + inflow 

DOF6 body states + flapping + lead-lag + inflow 

Do not do any reduction 


// (computed in linearize.exe:) 

// F(ns,ns) -- full state matrix 

// GP(ns,4) -- full control matrix (inputs are trim controls) 

// H(14,ns) 

// DP(14,4) -- output matrices for body velocities & accelerations 

(???) 

II 

II OUTPUTS: 

// FQ(n f n) 

// GQ (n, 4) -- state and control matrix of reduced model if ndof 

<> 0 

// HQ(???,n) 

// DQ(???,4) -- steady state output matrices for reduced model 

//Conqpute indices for quasistaic reduction 

ixSsStatelndices(topsolve,world_model_heli_body_bodydof) "Body state 
indices"; 

ix6r=ix6([13:18 4 5]) "??? body state indices"; 

idrtrn=StateIndices(topsolve,drivesolve) "Drivetrain state indices"; 
if (world_cpg_havembc == 1), 

irotor=StateIndices(topsolve,model_heli_rotor_mbc); 
iflap=irotor([1:4 9:12]) "Flapping state indices"; // %%% will 
change 

ilag=irotor([5:8 13:16]) "Lead-lag state indices"; // %%% will 
change 
else 

irotor= [] ; 
if lap= [] ; 
ilag* [] ; 

end 

iinflowsStatelndices(topsolve,model_heli_inflow_ul) "Inflow state 
indices"; 

// Determine dynamic states based on NDOF: 

if ndof == 6, 

idyn=[ix6r] ; 
elseif ndof == 10, 

idyn=[ix6r;iflap]; 
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els^if ndof 11, 

idyn-[ix6r;iflap;iinflow]; 
elseif ndof 15, 

idyn-[ix6r;iflap;ilag;iinflow]; 
elseif ndof ■■ 0, 

// ??? 

else 

error("NDOF » "$ndof$"; must be 6,10,11 or 15"); 

end 

// Always eliminate drivetrain and ??? dof6 states 
ielim« [idrtm;ix$ ([1:3 6:12])]; 

// Compute reduced model: 
if(ndof <> 0), 

[fq,gq,hq,dq]«qsreduce(f,gp,idyn,ielim, 0); 

end; 

// end of qsreduce.exe 




//////////// results of analysis.def ///////////////// 

Trim Matrix^ 


0.0312 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


0.0000 

-0.0312 

0.0000 

0.0000 

0.0000 

0.0000 


0.0000 

0.0000 

0.1276 

0.0000 

0.0000 

0.0000 


0.0000 

0.0000 

0.0000 

-0.5613 

0.0000 

0.0000 


0.0000 

0.0000 

0.0000 

0.0000 

1.8720 

0.0000 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0.5594 


The reduced 
-0.0300 

linear characteristic 
0.0091 -0.0591 

matrix 

-1.2212 

1.5360 

-0.2055 

0.0000 

-32.0124 

0.0019 

0.0307 

-0.0095 

-0.7664 

-0.0879 

0.4489 

32.0036 

-0.0810 

-0.0456 

-0.0067 

-0.7201 

-2.1623 

-1.2583 

2.7010 

0.7461 

3 .4700 

0.0015 

-0.0462 

0.0134 

-6.6040 

-1.1185 

1.2524 

0.0000 

0.0000 

0.0036 

0.0100 

-0.0025 

0.2859 

-1.6059 

-0.3555 

0.0000 

0.0000 

-0.0057 

0.0149 

-0.0193 

0.0024 

0.0888 

-0.5162 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0025 

-0.1084 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.9997 

0.0233 

0.0000 

0.0000 


The reduced linear control matrix 


1.8661 

-0.1586 

0.4224 

-0.8579 

0.3245 

0.2813 

-2.4906 

-0.2769 

6.7358 

-0.0642 

1.0253 

-8.1918 

0.1939 

1.8587 

-2.1082 

0.1150 

-0.5083 

0.0599 

0.6609 

0.1455 

0.1236 

-0.0548 

1.7156 

0.0738 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


The characteristic equation coefficients of the linear matrix fq 


1.0000 

+ 

O.OOOOi 

9.4455 

+ 

O.OOOOi 

21.4968 

+ 

O.OOOOi 

18.5608 

+ 

O.OOOOi 

9.8123 

+ 

O.OOOOi 

4.5450 

- 

O.OOOOi 

1.4522 

- 

O.OOOOi 

0.2193 

+ 

O.OOOOi 

0-0759 

+ 

O.OOOOi 
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The eigenvectors of the linear characteristic matrix 
columns 1 thru 4 


-0.1651 

+ 

0.0627i 

0.1248 

- 

0.9841i 

-1.2950 

- 

0.0458i 

0.4863 

- 

1.0047i 

-0.5704 

♦ 

0.2168i 

0.0238 

- 

0.1875i 

-0.6219 

+ 

0.0176i 

0.4434 

- 

1.4703i 

-0.2499 

■f 

0.0950i 

-0.0241 

+ 

0.1903i 

-0.3135 

+ 

0.9051i 

-0.1009 

+ 

0.0154i 

-0.6684 

+ 

0.2540i 

0.0023 

- 

O.OlSli 

-0.0164 

+ 

0.0062i 

-0.0052 

♦ 

0.0069 i 

0.0396 

- 

O.OlSli 

-0.0104 

+ 

0.0820i 

0.0140 

- 

0.0073i 

0.0036 

- 

0.0059i 

0.0001 

- 

O.OOOOi 

0.0007 

- 

0.0053i 

-0.0849 

+ 

0.0114i 

-0.0101 

- 

0.0243i 

0.1017 

- 

0.0387i 

-0.0013 

+ 

O.OlOli 

0.0130 

- 

0.0039i 

0.0206 

+ 

0.0087i 

-0.0060 + 
columns ! 

0.0023i 

5 thru 8 

0.0060 

“ 

0.0476i 

-0.0213 

+ 

0.0048i 

-0.0141 

• 

0.0071i 

3.3901 

+ 

6.8954i 

9.2088 

- 

0.1581i 

4.6426 


24.9243i 

21.3797 

-! 

53.2807i 

-2.2410 

- 

0.6718i 

4.4177 

+ 

0.2052i 

11.6962 


32.8607i 

-16.3188 

«• 6.2655i 

-0.4637 

- 

0.5540i 

2.1109 

+ 

6.4724i 

1.6378 

♦ 

1.6409i 

-3.1432 

♦ 

4.3881i 

-0.0040 

- 

0.0015i 

0.1155 

+ 

0.0459i 

0.0036 

♦ 

0.1960i 

-0.0288 

♦ 

0.0132i 

0.0155 

+ 

0.0141i 

-0.0984 

- 

0.0538i 

0.0132 

- 

0.1568i 

0.1073 

- 

0.1140i 

-0.0889 

- 

0.0234i 

0.6018 

+ 

0.0922i 

0.5115 

- 

0.3090i 

-0.6494 

+ 

0.2244i 

0. OOf-5 

- 

0.0198i 

-0.0921 

- 

0.0294i 

-0.4941 

- 

0.1196i 

0.0596 

+ 

0.1438i 

0.0544 

- 

0.0416i 

0.1509 

+ 

0.0365i 

0.3539 

+ 

0.0597i 

0.4285 

+ 

0.2791i 


The eigenvalues of the linear characteristic matrix 


-6.5720 - 
-1.7182 - 
-0.6057 + 
-0.0069 + 
0.0349 + 
-0.6057 - 
-0.0069 - 
0.0349 - 


O.OOOOi 
O.OOOOi 
0.1953i 
0.4620i 
0.2765i 
0.1953i 
0.4620i 
0.2765i 


// end of analysis 


results 


// 
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